Newer
Older
1001
1002
1003
1004
1005
1006
1007
1008
1009
1010
1011
1012
1013
1014
1015
1016
1017
1018
1019
1020
1021
1022
1023
1024
1025
1026
1027
1028
1029
1030
1031
1032
1033
1034
1035
1036
1037
1038
1039
1040
1041
1042
1043
1044
1045
1046
1047
1048
1049
1050
1051
1052
1053
1054
1055
1056
1057
1058
1059
1060
1061
1062
1063
1064
1065
1066
1067
1068
1069
1070
1071
1072
1073
1074
1075
1076
1077
1078
1079
1080
1081
1082
1083
1084
1085
1086
1087
1088
1089
1090
1091
1092
1093
1094
1095
1096
1097
1098
1099
1100
1101
1102
1103
1104
1105
1106
1107
1108
1109
1110
1111
1112
1113
1114
1115
1116
1117
1118
1119
1120
1121
1122
1123
1124
1125
1126
1127
1128
1129
1130
1131
1132
1133
1134
1135
1136
1137
1138
1139
1140
1141
1142
1143
1144
1145
1146
1147
1148
1149
1150
1151
1152
1153
1154
1155
1156
1157
1158
1159
1160
1161
1162
1163
1164
1165
1166
1167
1168
1169
1170
1171
1172
1173
1174
1175
1176
1177
1178
1179
1180
1181
1182
1183
1184
1185
1186
1187
1188
1189
1190
1191
1192
1193
1194
1195
1196
1197
1198
1199
1200
1201
1202
1203
1204
1205
1206
1207
1208
1209
1210
1211
1212
1213
1214
1215
1216
1217
1218
1219
1220
1221
1222
1223
1224
1225
1226
1227
1228
1229
1230
1231
1232
1233
1234
1235
1236
1237
1238
1239
1240
1241
1242
1243
1244
1245
1246
1247
1248
1249
1250
1251
1252
1253
1254
1255
1256
1257
1258
1259
1260
1261
1262
1263
1264
1265
1266
1267
1268
1269
1270
1271
1272
1273
1274
1275
1276
1277
1278
1279
1280
1281
1282
1283
1284
1285
1286
1287
1288
1289
1290
1291
1292
1293
1294
1295
1296
1297
1298
1299
1300
1301
1302
1303
1304
1305
1306
1307
1308
1309
1310
1311
1312
1313
1314
1315
1316
1317
1318
1319
1320
1321
1322
1323
1324
1325
1326
1327
1328
1329
1330
1331
1332
1333
1334
1335
1336
1337
1338
1339
1340
1341
1342
1343
1344
1345
1346
1347
1348
1349
1350
1351
1352
1353
1354
1355
1356
1357
1358
1359
1360
1361
1362
1363
1364
1365
1366
1367
1368
1369
1370
1371
1372
1373
1374
1375
1376
1377
1378
1379
1380
1381
1382
1383
1384
1385
1386
1387
1388
1389
1390
1391
1392
1393
1394
1395
1396
1397
1398
1399
1400
1401
1402
1403
1404
1405
1406
1407
1408
1409
1410
1411
1412
1413
1414
1415
1416
1417
1418
1419
1420
1421
1422
1423
1424
1425
1426
1427
1428
1429
1430
1431
1432
1433
1434
1435
1436
1437
1438
1439
1440
1441
1442
1443
1444
1445
1446
1447
1448
1449
1450
1451
1452
1453
ALLOCATE(XLBXUM(4,NJU,NKU))
ALLOCATE(XLBXVM(2,NJU,NKU))
ALLOCATE(XLBXWM(2,NJU,NKU))
ALLOCATE(XLBXTHM(2,NJU,NKU))
END IF
!
IF ( NRR > 0 ) THEN
IF ( LHORELAX_RV .OR. LHORELAX_RC .OR. LHORELAX_RR .OR. LHORELAX_RI &
.OR. LHORELAX_RS .OR. LHORELAX_RG .OR. LHORELAX_RH &
) THEN
NSIZELBXR_ll=2* NRIMX+2
ALLOCATE(XLBXRM(2*NRIMX+2,NJU,NKU,NRR))
ELSE
NSIZELBXR_ll=2
ALLOCATE(XLBXRM(2,NJU,NKU,NRR))
ENDIF
ELSE
NSIZELBXR_ll=0
ALLOCATE(XLBXRM(0,0,0,0))
END IF
!
IF ( NSV > 0 ) THEN
IF ( ANY( LHORELAX_SV(:)) ) THEN
NSIZELBXSV_ll=2* NRIMX+2
ALLOCATE(XLBXSVM(2*NRIMX+2,NJU,NKU,NSV))
ELSE
NSIZELBXSV_ll=2
ALLOCATE(XLBXSVM(2,NJU,NKU,NSV))
END IF
ELSE
NSIZELBXSV_ll=0
ALLOCATE(XLBXSVM(0,0,0,0))
END IF
!
ELSE ! 3D case
!
CALL GET_SIZEX_LB(CLUOUT,NIMAX_ll,NJMAX_ll,NRIMX, &
IISIZEXF,IJSIZEXF,IISIZEXFU,IJSIZEXFU, &
IISIZEX4,IJSIZEX4,IISIZEX2,IJSIZEX2)
CALL GET_SIZEY_LB(CLUOUT,NIMAX_ll,NJMAX_ll,NRIMY, &
IISIZEYF,IJSIZEYF,IISIZEYFV,IJSIZEYFV, &
IISIZEY4,IJSIZEY4,IISIZEY2,IJSIZEY2)
!
IF ( LHORELAX_UVWTH ) THEN
NSIZELBX_ll=2*NRIMX+2
NSIZELBXU_ll=2*NRIMX+2
NSIZELBY_ll=2*NRIMY+2
NSIZELBYV_ll=2*NRIMY+2
ALLOCATE(XLBXUM(IISIZEXFU,IJSIZEXFU,NKU))
ALLOCATE(XLBYUM(IISIZEYF,IJSIZEYF,NKU))
ALLOCATE(XLBXVM(IISIZEXF,IJSIZEXF,NKU))
ALLOCATE(XLBYVM(IISIZEYFV,IJSIZEYFV,NKU))
ALLOCATE(XLBXWM(IISIZEXF,IJSIZEXF,NKU))
ALLOCATE(XLBYWM(IISIZEYF,IJSIZEYF,NKU))
ALLOCATE(XLBXTHM(IISIZEXF,IJSIZEXF,NKU))
ALLOCATE(XLBYTHM(IISIZEYF,IJSIZEYF,NKU))
ELSE
NSIZELBX_ll=2
NSIZELBXU_ll=4
NSIZELBY_ll=2
NSIZELBYV_ll=4
ALLOCATE(XLBXUM(IISIZEX4,IJSIZEX4,NKU))
ALLOCATE(XLBYUM(IISIZEY2,IJSIZEY2,NKU))
ALLOCATE(XLBXVM(IISIZEX2,IJSIZEX2,NKU))
ALLOCATE(XLBYVM(IISIZEY4,IJSIZEY4,NKU))
ALLOCATE(XLBXWM(IISIZEX2,IJSIZEX2,NKU))
ALLOCATE(XLBYWM(IISIZEY2,IJSIZEY2,NKU))
ALLOCATE(XLBXTHM(IISIZEX2,IJSIZEX2,NKU))
ALLOCATE(XLBYTHM(IISIZEY2,IJSIZEY2,NKU))
END IF
!
IF ( NRR > 0 ) THEN
IF ( LHORELAX_RV .OR. LHORELAX_RC .OR. LHORELAX_RR .OR. LHORELAX_RI &
.OR. LHORELAX_RS .OR. LHORELAX_RG .OR. LHORELAX_RH &
) THEN
NSIZELBXR_ll=2*NRIMX+2
NSIZELBYR_ll=2*NRIMY+2
ALLOCATE(XLBXRM(IISIZEXF,IJSIZEXF,NKU,NRR))
ALLOCATE(XLBYRM(IISIZEYF,IJSIZEYF,NKU,NRR))
ELSE
NSIZELBXR_ll=2
NSIZELBYR_ll=2
ALLOCATE(XLBXRM(IISIZEX2,IJSIZEX2,NKU,NRR))
ALLOCATE(XLBYRM(IISIZEY2,IJSIZEY2,NKU,NRR))
ENDIF
ELSE
NSIZELBXR_ll=0
NSIZELBYR_ll=0
ALLOCATE(XLBXRM(0,0,0,0))
ALLOCATE(XLBYRM(0,0,0,0))
END IF
!
IF ( NSV > 0 ) THEN
IF ( ANY( LHORELAX_SV(:)) ) THEN
NSIZELBXSV_ll=2*NRIMX+2
NSIZELBYSV_ll=2*NRIMY+2
ALLOCATE(XLBXSVM(IISIZEXF,IJSIZEXF,NKU,NSV))
ALLOCATE(XLBYSVM(IISIZEYF,IJSIZEYF,NKU,NSV))
ELSE
NSIZELBXSV_ll=2
NSIZELBYSV_ll=2
ALLOCATE(XLBXSVM(IISIZEX2,IJSIZEX2,NKU,NSV))
ALLOCATE(XLBYSVM(IISIZEY2,IJSIZEY2,NKU,NSV))
END IF
ELSE
NSIZELBXSV_ll=0
NSIZELBYSV_ll=0
ALLOCATE(XLBXSVM(0,0,0,0))
ALLOCATE(XLBYSVM(0,0,0,0))
END IF
END IF
!
!
!-------------------------------------------------------------------------------
!
!* 5. INITIALIZE ALL THE MODEL VARIABLES
! ----------------------------------
!
!
!* 5.1 Grid variables and RS localization:
!
!* 5.1.1 Horizontal Spatial grid :
!
IF( LEN_TRIM(CPGD_FILE) /= 0 ) THEN
!--------------------------------------------------------
! the MESONH horizontal grid will be read in the PGD_FILE
!--------------------------------------------------------
CALL READ_HGRID(1,CPGD_FILE,YPGD_NAME,YPGD_DAD_NAME,YPGD_TYPE)
! control the cartesian option
IF( LCARTESIAN ) THEN
WRITE(NLUOUT,FMT=*) 'PREP_IDEAL_CASE : IN GENERAL, THE USE OF A PGD_FILE &
& IMPLIES THAT YOU MUST TAKE INTO ACCOUNT THE EARTH SPHERICITY'
WRITE(NLUOUT,FMT=*) 'NEVERTHELESS, LCARTESIAN HAS BEEN KEPT TO TRUE'
END IF
!
!* use of the externalized surface
!
CSURF = "EXTE"
!
! determine whether the model is flat or no
!
IF( ABS( MAXVAL(XZS(NIB:NIU-JPHEXT,NJB:NJU-JPHEXT)) ) < 1.E-10 ) THEN
LFLAT=.TRUE.
ELSE
LFLAT=.FALSE.
END IF
!
ELSE
!------------------------------------------------------------------------
! the MESONH horizontal grid is built from the PRE_IDEA1.nam informations
!------------------------------------------------------------------------
!
ALLOCATE(XXHAT(NIU),XYHAT(NJU))
!
! define the grid localization at the earth surface by the central point
! coordinates
!
IF (XLONCEN/=XUNDEF .OR. XLATCEN/=XUNDEF) THEN
IF (XLONCEN/=XUNDEF .AND. XLATCEN/=XUNDEF) THEN
!
! it should be noted that XLATCEN and XLONCEN refer to a vertical
! vorticity point and (XLATORI, XLONORI) refer to the mass point of
! conformal coordinates (0,0). This is to allow the centering of the model in
! a non-cyclic configuration regarding to XLATCEN or XLONCEN.
!
ALLOCATE(ZXHAT_ll(NIMAX_ll+2*JPHEXT),ZYHAT_ll(NJMAX_ll+2*JPHEXT))
ZXHAT_ll=0.
ZYHAT_ll=0.
CALL SM_LATLON(XLATCEN,XLONCEN, &
-XDELTAX*(NIMAX_ll/2-0.5+JPHEXT), &
-XDELTAY*(NJMAX_ll/2-0.5+JPHEXT), &
XLATORI,XLONORI)
DEALLOCATE(ZXHAT_ll,ZYHAT_ll)
!
WRITE(NLUOUT,FMT=*) 'PREP_IDEAL_CASE : XLATORI=' , XLATORI, &
' XLONORI= ', XLONORI
ELSE
WRITE(NLUOUT,FMT=*) 'PREP_IDEAL_CASE : LATITUDE AND LONGITUDE OF THE CENTER &
& POINT MUST BE INITIALIZED ALL TOGETHER OR NOT'
WRITE(NLUOUT,FMT=*) '-> JOB ABORTED'
!callabortstop
CALL CLOSE_ll(CLUOUT,IOSTAT=IRESP)
CALL ABORT
STOP
END IF
END IF
!
IF (NPROC > 1) THEN
CALL GET_DIM_EXT_ll('B',IXDIM,IYDIM)
IBEG = IXOR-JPHEXT-1
IEND = IBEG+IXDIM-1
XXHAT(:) = (/ (FLOAT(JLOOP)*XDELTAX, JLOOP=IBEG,IEND) /)
IBEG = IYOR-JPHEXT-1
IEND = IBEG+IYDIM-1
XYHAT(:) = (/ (FLOAT(JLOOP)*XDELTAY, JLOOP=IBEG,IEND) /)
!
ELSE
XXHAT(:) = (/ (FLOAT(JLOOP-NIB)*XDELTAX, JLOOP=1,NIU) /)
XYHAT(:) = (/ (FLOAT(JLOOP-NJB)*XDELTAY, JLOOP=1,NJU) /)
END IF
END IF
!
!* 5.1.2 Orography and Gal-Chen Sommerville transformation :
!
IF ( LEN_TRIM(CPGD_FILE) == 0 .OR. .NOT. LREAD_ZS) THEN
SELECT CASE(CZS) ! 'FLAT' or 'SINE' or 'BELL'
CASE('FLAT')
LFLAT = .TRUE.
IF (XHMAX==XUNDEF) THEN
XZS(:,:) = 0.
ELSE
XZS(:,:) = XHMAX
END IF
CASE('SINE') ! sinus-shaped orography
IF (XHMAX==XUNDEF) XHMAX=300.
LFLAT =.FALSE.
XZS(:,:) = XHMAX & ! three-dimensional case
*SPREAD((/((SIN((XPI/(NIMAX_ll+2*JPHEXT-1))*JLOOP)**2)**NEXPX,JLOOP=IXOR-1,IXOR+NIU-2)/),2,NJU) &
*SPREAD((/((SIN((XPI/(NJMAX_ll+2*JPHEXT-1))*JLOOP)**2)**NEXPY,JLOOP=IYOR-1,IYOR+NJU-2)/),1,NIU)
IF(L1D) THEN ! one-dimensional case
XZS(:,:) = XHMAX
END IF
CASE('BELL') ! bell-shaped orography
IF (XHMAX==XUNDEF) XHMAX=300.
LFLAT = .FALSE.
IF(.NOT.L2D) THEN ! three-dimensional case
XZS(:,:) = XHMAX / ( 1. &
+ ( (SPREAD(XXHAT(1:NIU),2,NJU) - FLOAT(NIZS) * XDELTAX) /XAX ) **2 &
+ ( (SPREAD(XYHAT(1:NJU),1,NIU) - FLOAT(NJZS) * XDELTAY) /XAY ) **2 ) **1.5
ELSE ! two-dimensional case
XZS(:,:) = XHMAX / ( 1. &
+ ( (SPREAD(XXHAT(1:NIU),2,NJU) - FLOAT(NIZS) * XDELTAX) /XAX ) **2 )
ENDIF
IF(L1D) THEN ! one-dimensional case
XZS(:,:) = XHMAX
END IF
CASE('COSI') ! (1+cosine)**4 shape
IF (XHMAX==XUNDEF) XHMAX=800.
LFLAT = .FALSE.
IF(L2D) THEN ! two-dimensional case
DO JILOOP = 1, NIU
ZDIST = XXHAT(JILOOP)-FLOAT(NIZS)*XDELTAX
IF( ABS(ZDIST)<(4.0*XAX) ) THEN
XZS(JILOOP,:) = (XHMAX/16.0)*( 1.0 + COS((XPI*ZDIST)/(4.0*XAX)) )**4
ELSE
XZS(JILOOP,:) = 0.0
ENDIF
END DO
ENDIF
CASE('SCHA') ! exp(-(x/a)**2)*cosine(pi*x/lambda)**2 shape
IF (XHMAX==XUNDEF) XHMAX=800.
LFLAT = .FALSE.
IF(L2D) THEN ! two-dimensional case
DO JILOOP = 1, NIU
ZDIST = XXHAT(JILOOP)-FLOAT(NIZS)*XDELTAX
IF( ABS(ZDIST)<(4.0*XAX) ) THEN
XZS(JILOOP,:) = XHMAX*EXP(-(ZDIST/XAY)**2)*COS((XPI*ZDIST)/XAX)**2
ELSE
XZS(JILOOP,:) = 0.0
ENDIF
END DO
ENDIF
CASE('DATA') ! discretized orography
LFLAT =.FALSE.
WRITE(NLUOUT,FMT=*) 'CZS="DATA", ATTEMPT TO READ ARRAY &
&XZS(NIB:NIU-JPHEXT:1,NJU-JPHEXT:NJB:-1) &
&starting from the first index'
CALL POSKEY(NLUPRE,NLUOUT,'ZSDATA')
DO JJLOOP = NJMAX_ll+2*JPHEXT-1,JPHEXT+1,-1 ! input like a map prior the sounding
READ(NLUPRE,FMT=*) ZZS_ll
IF ( ( JJLOOP <= ( NJU-JPHEXT + IYOR-1 ) ) .AND. ( JJLOOP >= ( NJB + IYOR-1 ) ) ) THEN
IJ = JJLOOP - ( IYOR-1 )
XZS(NIB:NIU-JPHEXT,IJ) = ZZS_ll(IXOR:IXOR + NIU-JPHEXT - NIB )
END IF
END DO
!
CASE DEFAULT ! undefined shape of orography
WRITE(NLUOUT,FMT=*) 'PREP_IDEAL_CASE: ERRONEOUS TERRAIN TYPE'
WRITE(NLUOUT,FMT=*) '-> JOB ABORTED'
!callabortstop
CALL CLOSE_ll(CLUOUT,IOSTAT=IRESP)
CALL ABORT
STOP
END SELECT
!
CALL ADD2DFIELD_ll(TZ_FIELDS_ll, XZS)
CALL UPDATE_HALO_ll(TZ_FIELDS_ll,IINFO_ll)
CALL CLEANLIST_ll(TZ_FIELDS_ll)
!
END IF
!
!IF( ( LEN_TRIM(CPGD_FILE) /= 0 ) .AND. .NOT.LFLAT .AND. &
! ((CLBCX(1) /= "OPEN" ) .OR. &
! (CLBCX(2) /= "OPEN" ) .OR. (CLBCY(1) /= "OPEN" ) .OR. &
! (CLBCY(2) /= "OPEN" )) ) THEN
! WRITE(NLUOUT,FMT=*) 'STOP:WITH A PGD FILE YOU CANNOT BE IN CYCLIC LBC'
! !callabortstop
! CALL CLOSE_ll(CLUOUT,IOSTAT=IRESP)
! CALL ABORT
! STOP
!END IF
!
IF (LWEST_ll()) THEN
DO JILOOP = 1,JPHEXT
XZS(JILOOP,:) = XZS(NIB,:)
END DO
END IF
IF (LEAST_ll()) THEN
DO JILOOP = NIU-JPHEXT+1,NIU
XZS(JILOOP,:)=XZS(NIU-JPHEXT,:)
END DO
END IF
IF (LSOUTH_ll()) THEN
DO JJLOOP = 1,JPHEXT
XZS(:,JJLOOP)=XZS(:,NJB)
END DO
END IF
IF (LNORTH_ll()) THEN
DO JJLOOP =NJU-JPHEXT+1,NJU
XZS(:,JJLOOP)=XZS(:,NJU-JPHEXT)
END DO
END IF
!
IF ( LEN_TRIM(CPGD_FILE) == 0 .OR. .NOT. LREAD_ZS) THEN
IF (LSLEVE) THEN
CALL ZSMT_PIC(NSLEVE,XSMOOTH_ZS)
ELSE
XZSMT(:,:) = 0.
END IF
END IF
!
IF (LCARTESIAN) THEN
CALL SM_GRIDCART(CLUOUT,XXHAT,XYHAT,XZHAT,XZS,LSLEVE,XLEN1,XLEN2,XZSMT,XDXHAT,XDYHAT,XZZ,XJ)
XMAP=1.
ELSE
CALL SM_GRIDPROJ(CLUOUT,XXHAT,XYHAT,XZHAT,XZS,LSLEVE,XLEN1,XLEN2,XZSMT,XLATORI,XLONORI, &
XMAP,XLAT,XLON,XDXHAT,XDYHAT,XZZ,XJ)
END IF
!* 5.4.1 metrics coefficients and update halos:
!
CALL METRICS(XMAP,XDXHAT,XDYHAT,XZZ,XDXX,XDYY,XDZX,XDZY,XDZZ)
!
CALL UPDATE_METRICS(CLBCX,CLBCY,XDXX,XDYY,XDZX,XDZY,XDZZ)
!
!* 5.1.3 Compute the localization in index space of the vertical profile
! in CSTN and RSOU cases :
!
IF (CTYPELOC =='LATLON' ) THEN
IF (.NOT.LCARTESIAN) THEN ! compute (x,y) if
CALL SM_XYHAT(XLATORI,XLONORI, & ! the localization
XLATLOC,XLONLOC,XXHATLOC,XYHATLOC) ! is given in latitude
ELSE ! and longitude
WRITE(NLUOUT,FMT=*) 'CTYPELOC CANNOT BE LATLON IN CARTESIAN GEOMETRY'
WRITE(NLUOUT,FMT=*) '-> JOB ABORTED'
!callabortstop
CALL CLOSE_ll(CLUOUT,IOSTAT=IRESP)
CALL ABORT
STOP
END IF
END IF
!
ALLOCATE(ZXHAT_ll(NIMAX_ll+ 2 * JPHEXT),ZYHAT_ll(NJMAX_ll+2 * JPHEXT))
CALL GATHERALL_FIELD_ll('XX',XXHAT,ZXHAT_ll,NRESP) !//
CALL GATHERALL_FIELD_ll('YY',XYHAT,ZYHAT_ll,NRESP) !//
IF (CTYPELOC /= 'IJGRID') THEN
NILOC = MINLOC(ABS(XXHATLOC-ZXHAT_ll(:)))
NJLOC = MINLOC(ABS(XYHATLOC-ZYHAT_ll(:)))
END IF
!
IF ( NILOC(1) == 1 ) NILOC = 2
IF ( NJLOC(1) == 1 ) NJLOC = 2
!
IF ( L1D .AND. (NILOC(1) /= NIB .OR. NJLOC(1) /= NJB) ) THEN
NILOC = NIB
NJLOC = NJB
WRITE(NLUOUT,FMT=*) 'FOR 1D CONFIGURATION, THE RS INFORMATIONS ARE TAKEN AT &
& I=NIB AND J=NJB (CENTRAL VERTICAL)'
END IF
!
IF ( L2D .AND. ( NJLOC(1) /= NJB) ) THEN
NJLOC = NJB
WRITE(NLUOUT,FMT=*) 'FOR 2D CONFIGURATION, THE RS INFORMATIONS ARE TAKEN AT &
& J=NJB (CENTRAL PLANE)'
END IF
!
!* 5.2 Prognostic variables (not multiplied by rhoJ) : u,v,w,theta,r
! and 1D anelastic reference state
!
IF(LPV_PERT .AND. .NOT.(LGEOSBAL)) THEN
WRITE(NLUOUT,FMT=*) 'FOR PV INVERSION, LGEOSBAL HAS TO BE TRUE'
!callabortstop
CALL CLOSE_ll(CLUOUT,IOSTAT=IRESP)
CALL ABORT
STOP
ENDIF
!
IF(LPV_PERT .AND. NPROC>1) THEN
WRITE(NLUOUT,FMT=*) 'PREP_IDEAL_CASE : THE USE OF A PV INVERSION HAS TO BE &
& PERFORMED WITH MONOPROCESSOR MODE'
WRITE(NLUOUT,FMT=*) '-> JOB ABORTED'
!callabortstop
CALL CLOSE_ll(CLUOUT,IOSTAT=IRESP)
CALL ABORT
STOP
ENDIF
!
!* 5.2.1 Use a Radiosounding : CIDEAL='RSOU''
!
IF (CIDEAL == 'RSOU') THEN
WRITE(NLUOUT,FMT=*) 'CIDEAL="RSOU", attempt to read DATE'
CALL POSKEY(NLUPRE,NLUOUT,'RSOU')
READ(NLUPRE,FMT=*) NYEAR,NMONTH,NDAY,XTIME
TDTCUR = DATE_TIME(DATE(NYEAR,NMONTH,NDAY),XTIME)
TDTEXP = TDTCUR
TDTSEG = TDTCUR
TDTMOD = TDTCUR
READ(NLUPRE,*) YKIND
BACKSPACE(NLUPRE) ! because YKIND read again in set_rsou
WRITE(NLUOUT,FMT=*) 'CIDEAL="RSOU", ATTEMPT TO PROCESS THE SOUNDING DATA'
IF (LGEOSBAL) THEN
CALL SET_RSOU(CEXPRE,CFUNU,CFUNV,NILOC(1),NJLOC(1),LBOUSS,LPV_PERT,&
LRMV_BL,XJ,LSHIFT,XCORIOZ)
ELSE
CALL SET_RSOU(CEXPRE,CFUNU,CFUNV,NILOC(1),NJLOC(1),LBOUSS,LPV_PERT,&
LRMV_BL,XJ,LSHIFT)
END IF
!
!* 5.2.2 N=cste and U(z) : CIDEAL='CSTN'
!
ELSE IF (CIDEAL == 'CSTN') THEN
WRITE(NLUOUT,FMT=*) 'CIDEAL="CSTN", attempt to read DATE'
CALL POSKEY(NLUPRE,NLUOUT,'CSTN')
READ(NLUPRE,FMT=*) NYEAR,NMONTH,NDAY,XTIME
TDTCUR = DATE_TIME(DATE(NYEAR,NMONTH,NDAY),XTIME)
TDTEXP = TDTCUR
TDTSEG = TDTCUR
TDTMOD = TDTCUR
WRITE(NLUOUT,FMT=*) 'CIDEAL="CSTN", ATTEMPT TO PROCESS THE SOUNDING DATA'
IF (LGEOSBAL) THEN
CALL SET_CSTN(CEXPRE,CFUNU,CFUNV,NILOC(1),NJLOC(1),LBOUSS,LPV_PERT,&
LRMV_BL,XJ,LSHIFT,XCORIOZ)
ELSE
CALL SET_CSTN(CEXPRE,CFUNU,CFUNV,NILOC(1),NJLOC(1),LBOUSS,LPV_PERT,&
LRMV_BL,XJ,LSHIFT)
END IF
!
END IF
!
!* 5.3 Forcing variables
!
IF (LFORCING) THEN
WRITE(NLUOUT,FMT=*) 'FORCING IS ENABLED, ATTEMPT TO SET FORCING FIELDS'
Gaelle Tanguy
committed
CALL POSKEY(NLUPRE,NLUOUT,'ZFRC ','PFRC')
1455
1456
1457
1458
1459
1460
1461
1462
1463
1464
1465
1466
1467
1468
1469
1470
1471
1472
1473
1474
1475
1476
1477
1478
1479
1480
1481
1482
1483
1484
1485
1486
1487
1488
1489
1490
1491
1492
1493
1494
1495
1496
1497
1498
1499
1500
1501
1502
CALL SET_FRC(CEXPRE)
END IF
!
!! ---------------------------------------------------------------------
! Modif PP ADV FRC
! 5.4.2 initialize profiles for adv forcings
IF (L2D_ADV_FRC) THEN
WRITE(NLUOUT,FMT=*) 'L2D_ADV_FRC IS SET TO TRUE'
WRITE(NLUOUT,FMT=*) 'ADVECTING FORCING USED IS USER MADE, NOT STANDARD ONE '
WRITE(NLUOUT,FMT=*) 'IT IS FOR 2D IDEALIZED WAM STUDY ONLY '
CALL POSKEY(NLUPRE,NLUOUT,'ZFRC_ADV')
CALL SET_ADVFRC(CEXPRE)
ENDIF
IF (L2D_REL_FRC) THEN
WRITE(NLUOUT,FMT=*) 'L2D_REL_FRC IS SET TO TRUE'
WRITE(NLUOUT,FMT=*) 'RELAXATION FORCING USED IS USER MADE, NOT STANDARD ONE '
WRITE(NLUOUT,FMT=*) 'IT IS FOR 2D IDEALIZED WAM STUDY ONLY '
CALL POSKEY(NLUPRE,NLUOUT,'ZFRC_REL')
CALL SET_RELFRC(CEXPRE)
ENDIF
!* 5.4 3D Reference state variables :
!
!
!* 5.4.1 metrics coefficients and update halos:
!
CALL METRICS(XMAP,XDXHAT,XDYHAT,XZZ,XDXX,XDYY,XDZX,XDZY,XDZZ)
!
CALL UPDATE_METRICS(CLBCX,CLBCY,XDXX,XDYY,XDZX,XDZY,XDZZ)
!
!* 5.4.2 3D reference state :
!
CALL SET_REF(0,'NIL',CLUOUT, &
XZZ,XZHAT,XJ,XDXX,XDYY,CLBCX,CLBCY, &
XREFMASS,XMASS_O_PHI0,XLINMASS, &
XRHODREF,XTHVREF,XRVREF,XEXNREF,XRHODJ )
!
!
!* 5.5.1 Absolute pressure :
!
!
!* 5.5.2 Total mass of dry air Md computation :
!
CALL TOTAL_DMASS(CLUOUT,XJ,XRHODREF,XDRYMASST)
!
!
!* 5.6 Complete prognostic variables (multipliy by rhoJ) at time t :
!
! U grid : gridpoint 2
IF (LWEST_ll()) XUT(1,:,:) = 2.*XUT(2,:,:) - XUT(3,:,:)
IF (LSOUTH_ll()) XVT(:,1,:) = 2.*XVT(:,2,:) - XVT(:,3,:)
!
!
!* 5.7 Larger scale fields initialization :
!
XLSUM(:,:,:) = XUT(:,:,:) ! these fields do not satisfy the
XLSVM(:,:,:) = XVT(:,:,:) ! lower boundary condition but are
XLSWM(:,:,:) = XWT(:,:,:) ! in equilibrium
XLSTHM(:,:,:)= XTHT(:,:,:)
XLSRVM(:,:,:)= XRT(:,:,:,1)
!
! enforce the vertical homogeneity under the ground and above the top of
! the model for the LS fields
!
XLSUM(:,:,NKB-1)=XLSUM(:,:,NKB)
XLSUM(:,:,NKU)=XLSUM(:,:,NKU-1)
XLSVM(:,:,NKB-1)=XLSVM(:,:,NKB)
XLSVM(:,:,NKU)=XLSVM(:,:,NKU-1)
XLSWM(:,:,NKB-1)=XLSWM(:,:,NKB)
XLSWM(:,:,NKU)=XLSWM(:,:,NKU-1)
XLSTHM(:,:,NKB-1)=XLSTHM(:,:,NKB)
XLSTHM(:,:,NKU)=XLSTHM(:,:,NKU-1)
IF ( NRR > 0 ) THEN
XLSRVM(:,:,NKB-1)=XLSRVM(:,:,NKB)
XLSRVM(:,:,NKU)=XLSRVM(:,:,NKU-1)
END IF
!
ILBX=SIZE(XLBXUM,1)
ILBY=SIZE(XLBYUM,2)
IF(LWEST_ll() .AND. .NOT. L1D) THEN
XLBXUM(1:NRIMX+1, :,:) = XUT(2:NRIMX+2, :,:)
XLBXVM(1:NRIMX+1, :,:) = XVT(1:NRIMX+1, :,:)
XLBXWM(1:NRIMX+1, :,:) = XWT(1:NRIMX+1, :,:)
XLBXTHM(1:NRIMX+1, :,:) = XTHT(1:NRIMX+1, :,:)
XLBXRM(1:NRIMX+1, :,:,:) = XRT(1:NRIMX+1, :,:,:)
ENDIF
IF(LEAST_ll() .AND. .NOT. L1D) THEN
XLBXUM(ILBX-NRIMX:ILBX,:,:) = XUT(NIU-NRIMX:NIU, :,:)
XLBXVM(ILBX-NRIMX:ILBX,:,:) = XVT(NIU-NRIMX:NIU, :,:)
XLBXWM(ILBX-NRIMX:ILBX,:,:) = XWT(NIU-NRIMX:NIU, :,:)
XLBXTHM(ILBX-NRIMX:ILBX,:,:) = XTHT(NIU-NRIMX:NIU, :,:)
XLBXRM(ILBX-NRIMX:ILBX,:,:,:) = XRT(NIU-NRIMX:NIU, :,:,:)
ENDIF
IF(LSOUTH_ll() .AND. .NOT. L1D .AND. .NOT. L2D) THEN
XLBYUM(:,1:NRIMY+1, :) = XUT(:,1:NRIMY+1, :)
XLBYVM(:,1:NRIMY+1, :) = XVT(:,2:NRIMY+2, :)
XLBYWM(:,1:NRIMY+1, :) = XWT(:,1:NRIMY+1, :)
XLBYTHM(:,1:NRIMY+1, :) = XTHT(:,1:NRIMY+1, :)
XLBYRM(:,1:NRIMY+1, :,:) = XRT(:,1:NRIMY+1, :,:)
ENDIF
IF(LNORTH_ll().AND. .NOT. L1D .AND. .NOT. L2D) THEN
XLBYUM(:,ILBY-NRIMY:ILBY,:) = XUT(:,NJU-NRIMY:NJU, :)
XLBYVM(:,ILBY-NRIMY:ILBY,:) = XVT(:,NJU-NRIMY:NJU, :)
XLBYWM(:,ILBY-NRIMY:ILBY,:) = XWT(:,NJU-NRIMY:NJU, :)
XLBYTHM(:,ILBY-NRIMY:ILBY,:) = XTHT(:,NJU-NRIMY:NJU, :)
XLBYRM(:,ILBY-NRIMY:ILBY,:,:) = XRT(:,NJU-NRIMY:NJU, :,:)
ENDIF
DO JSV = 1, NSV
IF(LWEST_ll() .AND. .NOT. L1D) &
XLBXSVM(1:NRIMX+1, :,:,JSV) = XSVT(1:NRIMX+1, :,:,JSV)
XLBXSVM(ILBX-NRIMX:ILBX,:,:,JSV) = XSVT(NIU-NRIMX:NIU, :,:,JSV)
IF(LSOUTH_ll() .AND. .NOT. L1D .AND. .NOT. L2D) &
XLBYSVM(:,1:NRIMY+1, :,JSV) = XSVT(:,1:NRIMY+1, :,JSV)
IF(LNORTH_ll() .AND. .NOT. L1D .AND. .NOT. L2D) &
XLBYSVM(:,ILBY-NRIMY:ILBY,:,JSV) = XSVT(:,NJU-NRIMY:NJU, :,JSV)
END DO
!
!
!* 5.8 Add a perturbation to a basic state :
!
IF(LPERTURB) CALL SET_PERTURB(CEXPRE)
!
!
!* 5.9 Anelastic correction and pressure:
!
IF ( .NOT. L1D ) CALL PRESSURE_IN_PREP(XDXX,XDYY,XDZX,XDZY,XDZZ)
!
!
!* 5.10 Compute THETA, vapor and cloud mixing ratio
!
IF (CIDEAL == 'RSOU') THEN
ALLOCATE(ZEXN(NIU,NJU,NKU))
ALLOCATE(ZT(NIU,NJU,NKU))
ALLOCATE(ZTHL(NIU,NJU,NKU))
ALLOCATE(ZRT(NIU,NJU,NKU))
ALLOCATE(ZCPH(NIU,NJU,NKU))
ALLOCATE(ZLVOCPEXN(NIU,NJU,NKU))
ALLOCATE(ZLSOCPEXN(NIU,NJU,NKU))
ALLOCATE(ZFRAC_ICE(NIU,NJU,NKU))
ALLOCATE(ZRSATW(NIU,NJU,NKU))
ALLOCATE(ZRSATI(NIU,NJU,NKU))
ZRT=XRT(:,:,:,1)+XRT(:,:,:,2)+XRT(:,:,:,4)
ZEXN=(XPABST/XP00) ** (XRD/XCPD)
ZT=XTHT*(XPABST/XP00)**(XRD/XCPD)
ZCPH=XCPD+ XCPV * XRT(:,:,:,1)+ XCL *XRT(:,:,:,2) + XCI * XRT(:,:,:,4)
ZLVOCPEXN = (XLVTT + (XCPV-XCL) * (ZT-XTT))/(ZCPH*ZEXN)
ZLSOCPEXN = (XLSTT + (XCPV-XCI) * (ZT-XTT))/(ZCPH*ZEXN)
ZTHL=XTHT-ZLVOCPEXN*XRT(:,:,:,2)-ZLSOCPEXN*XRT(:,:,:,4)
DEALLOCATE(ZEXN)
DEALLOCATE(ZT)
DEALLOCATE(ZCPH)
DEALLOCATE(ZLVOCPEXN)
DEALLOCATE(ZLSOCPEXN)
CALL TH_R_FROM_THL_RT_3D('T',ZFRAC_ICE,XPABST,ZTHL,ZRT,XTHT,XRT(:,:,:,1), &
XRT(:,:,:,2),XRT(:,:,:,4),ZRSATW, ZRSATI)
DEALLOCATE(ZTHL)
DEALLOCATE(ZRT)
! Coherence test
IF ((.NOT. LUSERI) ) THEN
WRITE(NLUOUT,FMT=*) "*********************************"
WRITE(NLUOUT,FMT=*) 'WARNING'
WRITE(NLUOUT,FMT=*) 'YOU HAVE LUSERI=FALSE '
WRITE(NLUOUT,FMT=*) ' BUT WITH YOUR RADIOSOUNDING Ri/=0'
WRITE(NLUOUT,FMT=*) MINVAL(XRT(:,:,:,4)),MAXVAL(XRT(:,:,:,4))
WRITE(NLUOUT,FMT=*) "*********************************"
ENDIF
ENDIF
IF ((.NOT. LUSERC)) THEN
WRITE(NLUOUT,FMT=*) "*********************************"
WRITE(NLUOUT,FMT=*) 'WARNING'
WRITE(NLUOUT,FMT=*) 'YOU HAVE LUSERC=FALSE '
WRITE(NLUOUT,FMT=*) 'BUT WITH YOUR RADIOSOUNDING RC/=0'
WRITE(NLUOUT,FMT=*) MINVAL(XRT(:,:,:,2)),MAXVAL(XRT(:,:,:,2))
WRITE(NLUOUT,FMT=*) "*********************************"
ENDIF
ENDIF
! on remet les bonnes valeurs pour NRR
IF(CCLOUD=='NONE') NRR=1
IF(CCLOUD=='REVE') NRR=2
END IF
!
!-------------------------------------------------------------------------------
!
!* 6. INITIALIZE SCALAR VARIABLES FOR CHEMISTRY
! -----------------------------------------
!
! before calling chemistry
CCONF = 'START'
1651
1652
1653
1654
1655
1656
1657
1658
1659
1660
1661
1662
1663
1664
1665
1666
1667
1668
1669
1670
1671
1672
1673
1674
1675
1676
1677
1678
1679
1680
1681
1682
1683
1684
1685
1686
1687
1688
1689
1690
1691
1692
1693
1694
1695
1696
1697
1698
1699
1700
1701
1702
1703
1704
1705
1706
1707
1708
1709
1710
1711
1712
1713
1714
1715
1716
1717
1718
1719
1720
1721
1722
1723
1724
1725
1726
1727
1728
1729
1730
1731
1732
1733
1734
CALL CLOSE_ll(CEXPRE,IOSTAT=NRESP) ! Close the EXPRE file
!
IF ( LCH_INIT_FIELD ) CALL CH_INIT_FIELD_n(1, NLUOUT, NVERB)
!
!-------------------------------------------------------------------------------
!
!* 7. WRITE THE FMFILE
! ----------------
!
CALL SECOND_MNH2(ZTIME1)
!
NNPRAR = 22 + 2*(NRR+NSV) & ! 22 = number of grid variables + reference
+ 8 + 17 ! state variables + dimension variables
! 2*(8+NRR+NSV) + 1 = number of prognostic
! variables at time t and t-dt
NTYPE=1
CDESFM=ADJUSTL(ADJUSTR(CINIFILE)//'.des')
!
CALL FMOPEN_ll(CINIFILE,'WRITE',CLUOUT,NNPRAR,NTYPE,NVERB,NNINAR,NRESP)
!
CALL WRITE_DESFM_n(1,CDESFM,CLUOUT)
!
#ifdef MNH_NCWRIT
NC_WRITE = LNETCDF
CALL WRITE_LFIFM_n(CINIFILE,' ') ! There is no DAD model for PREP_IDEAL_CASE
IF ( LNETCDF ) THEN
DEF_NC=.FALSE.
CALL WRITE_LFIFM_n(CINIFILE,' ') ! There is no DAD model for PREP_IDEAL_CASE
DEF_NC=.TRUE.
END IF
#else
CALL WRITE_LFIFM_n(CINIFILE,' ') ! There is no DAD model for PREP_IDEAL_CASE
#endif
!
CALL SECOND_MNH2(ZTIME2)
!
XT_STORE = XT_STORE + ZTIME2 - ZTIME1
!
!-------------------------------------------------------------------------------
!
!* 8. EXTERNALIZED SURFACE
! --------------------
!
!
IF (CSURF =='EXTE') THEN
IF (LEN_TRIM(CINIFILEPGD)==0) THEN
IF (LEN_TRIM(CPGD_FILE)/=0) THEN
CINIFILEPGD=CPGD_FILE
ELSE
WRITE(NLUOUT,FMT=*) 'STOP : CINIFILEPGD needed in NAM_LUNITn'
!callabortstop
CALL CLOSE_ll(CLUOUT,IOSTAT=IRESP)
CALL ABORT
STOP
ENDIF
ENDIF
CALL ALLOC_SURFEX(1)
CALL READ_ALL_NAMELISTS('MESONH','PRE',.FALSE.)
! Switch to model 1 surface variables
CALL GOTO_SURFEX(1,.TRUE.)
!* definition of physiographic fields
! computed ...
IF (LEN_TRIM(CPGD_FILE)==0 .OR. .NOT. LREAD_GROUND_PARAM) THEN
CPGDFILE = CINIFILE
CALL PGD_GRID_SURF_ATM('MESONH',CINIFILE,'MESONH',.TRUE.)
CALL SPLIT_GRID('MESONH')
CALL PGD_SURF_ATM ('MESONH',CINIFILE,'MESONH',.TRUE.)
CPGDFILE = CINIFILEPGD
ELSE
! ... or read from file.
CPGDFILE = CPGD_FILE
CALL INIT_PGD_SURF_ATM('MESONH','PGD', &
' ',' ',&
TDTCUR%TDATE%YEAR, TDTCUR%TDATE%MONTH, &
TDTCUR%TDATE%DAY, TDTCUR%TIME )
!
END IF
!
!* forces orography from atmospheric file
IF (.NOT. LREAD_ZS) CALL MNHPUT_ZS_n
!
! on ecrit un nouveau fichier PGD que s'il n'existe pas
IF (LEN_TRIM(CPGD_FILE)==0 .OR. .NOT. LREAD_GROUND_PARAM) THEN
!* writing of physiographic fields in the file
CSTORAGE_TYPE='PG'
CALL FMOPEN_ll(CINIFILEPGD,'WRITE',CLUOUT,NNPRAR,NTYPE,NVERB,NNINAR,NRESP)
CALL FMWRIT(CINIFILEPGD,'PROGRAM ',CLUOUT,'--',CPROGRAM,0,1,' ',NRESP)
CALL FMWRIT(CINIFILEPGD,'SURF ',CLUOUT,'--','EXTE',0,1,' ',NRESP)
CALL FMWRIT(CINIFILEPGD,'L1D ',CLUOUT,'--',L1D,0,1,' ',NRESP)
CALL FMWRIT(CINIFILEPGD,'L2D ',CLUOUT,'--',L2D,0,1,' ',NRESP)
CALL FMWRIT(CINIFILEPGD,'PACK ',CLUOUT,'--',LPACK,0,1,' ',NRESP)
CALL WRITE_HGRID(1,CINIFILEPGD,' ')
NC_FILE='sf1'
NC_WRITE=LNETCDF
CALL WRITE_PGD_SURF_ATM_n('MESONH')
IF ( LNETCDF ) THEN
DEF_NC=.FALSE.
CALL WRITE_PGD_SURF_ATM_n('MESONH')
DEF_NC=.TRUE.
NC_WRITE = .FALSE.
END IF
#else
CALL FMWRIT(CINIFILEPGD,'PROGRAM ',CLUOUT,'--',CPROGRAM,0,1,' ',NRESP)
CALL FMWRIT(CINIFILEPGD,'SURF ',CLUOUT,'--','EXTE',0,1,' ',NRESP)
CALL FMWRIT(CINIFILEPGD,'L1D ',CLUOUT,'--',L1D,0,1,' ',NRESP)
CALL FMWRIT(CINIFILEPGD,'L2D ',CLUOUT,'--',L2D,0,1,' ',NRESP)
CALL FMWRIT(CINIFILEPGD,'PACK ',CLUOUT,'--',LPACK,0,1,' ',NRESP)
CALL WRITE_HGRID(1,CINIFILEPGD,' ')
CALL WRITE_PGD_SURF_ATM_n('MESONH')
#endif
CSTORAGE_TYPE='TT'
ENDIF
1766
1767
1768
1769
1770
1771
1772
1773
1774
1775
1776
1777
1778
1779
1780
1781
1782
1783
1784
1785
1786
1787
1788
1789
1790
1791
1792
1793
1794
1795
1796
1797
1798
1799
1800
1801
1802
1803
1804
1805
1806
1807
1808
1809
1810
1811
1812
1813
1814
1815
1816
1817
1818
1819
1820
1821
1822
1823
1824
1825
1826
1827
1828
1829
1830
1831
1832
1833
1834
1835
1836
1837
1838
1839
1840
1841
1842
1843
1844
1845
1846
1847
1848
1849
1850
1851
1852
1853
1854
1855
1856
1857
1858
1859
1860
1861
1862
1863
1864
!
!* deallocation of physiographic fields
CALL DEALLOC_SURF_ATM_n
!
!* rereading of physiographic fields and definition of prognostic fields
!* writing of all surface fields
COUTFMFILE = CINIFILE
CALL PREP_SURF_MNH(' ',' ')
CALL DEALLOC_SURFEX
ELSE
CSURF = "NONE"
END IF
!
!-------------------------------------------------------------------------------
!
!* 9. CLOSES THE FILE
! ---------------
!
IF (CSURF =='EXTE' .AND. (LEN_TRIM(CPGD_FILE)==0 .OR. .NOT. LREAD_GROUND_PARAM)) THEN
CALL FMCLOS_ll(CINIFILEPGD,'KEEP',CLUOUT,NRESP)
ENDIF
CALL FMCLOS_ll(CINIFILE,'KEEP',CLUOUT,NRESP)
IF( LEN_TRIM(CPGD_FILE) /= 0 ) THEN
CALL FMCLOS_ll(CPGD_FILE,'KEEP',CLUOUT,NRESP)
ENDIF
!
!
!-------------------------------------------------------------------------------
!
!* 10. PRINTS ON OUTPUT-LISTING
! ------------------------
!
IF (NVERB >= 5) THEN
WRITE(NLUOUT,FMT=*) 'PREP_IDEAL_CASE: LCARTESIAN,CIDEAL,CZS=', &
LCARTESIAN,CIDEAL,CZS
WRITE(NLUOUT,FMT=*) 'PREP_IDEAL_CASE: LUSERV=',LUSERV
WRITE(NLUOUT,FMT=*) 'PREP_IDEAL_CASE: XLON0,XLAT0,XBETA,XRPK,XLONORI,XLATORI=', &
XLON0,XLAT0,XBETA,XRPK,XLONORI,XLATORI
WRITE(NLUOUT,FMT=*) 'PREP_IDEAL_CASE: XDELTAX,XDELTAY=',XDELTAX,XDELTAY
WRITE(NLUOUT,FMT=*) 'PREP_IDEAL_CASE: NVERB=',NVERB
IF(LCARTESIAN) THEN
WRITE(NLUOUT,FMT=*) 'PREP_IDEAL_CASE: No map projection used.'
ELSE
IF (XRPK == 1.) THEN
WRITE(NLUOUT,FMT=*) 'PREP_IDEAL_CASE: Polar stereo used.'
ELSE IF (XRPK == 0.) THEN
WRITE(NLUOUT,FMT=*) 'PREP_IDEAL_CASE: Mercator used.'
ELSE
WRITE(NLUOUT,FMT=*) 'PREP_IDEAL_CASE: Lambert used, cone factor=',XRPK
END IF
END IF
END IF
!
IF (NVERB >= 5) THEN
WRITE(NLUOUT,FMT=*) 'PREP_IDEAL_CASE: IIB, IJB, IKB=',NIB,NJB,NKB
WRITE(NLUOUT,FMT=*) 'PREP_IDEAL_CASE: IIU, IJU, IKU=',NIU,NJU,NKU
END IF
!
!
!* 28.1 print statistics!
!
!
CALL SECOND_MNH2(ZTIME2)
XT_START=XT_START+ZTIME2-ZEND
!
! Set File Timing OUTPUT
!
CALL SET_ILUOUT_TIMING(NLUOUT)
!
! Compute global time
!
CALL TIME_STAT_ll(XT_START,ZTOT)
!
!
IMI = 1
CALL TIME_HEADER_ll(IMI)
!
CALL TIME_STAT_ll(XT_STORE,ZTOT, ' STORE-FIELDS','=')
CALL TIMING_SEPARATOR('+')
CALL TIMING_SEPARATOR('+')
WRITE(YMI,FMT="(I0)") IMI
CALL TIME_STAT_ll(XT_START,ZTOT, ' MODEL'//YMI,'+')
CALL TIMING_SEPARATOR('+')
CALL TIMING_SEPARATOR('+')
CALL TIMING_SEPARATOR('+')
WRITE(NLUOUT,FMT=*) ' '
WRITE(NLUOUT,FMT=*) '****************************************************'
WRITE(NLUOUT,FMT=*) '* PREP_IDEAL_CASE: PREP_IDEAL_CASE ENDS CORRECTLY. *'
WRITE(NLUOUT,FMT=*) '****************************************************'
!
CALL CLOSE_ll(CLUOUT,IOSTAT=NRESP)
CALL END_PARA_ll(IINFO_ll)
!
!
!callabortstop
!JUAN CALL ABORT
STOP
!
END PROGRAM PREP_IDEAL_CASE