@@ -1272,32 +1272,12 @@ SUBROUTINE open_met_file(dels,koffset,kend,spinup, TFRZ)
1272
1272
' in ' // TRIM (filename% met)// ' (SUBROUTINE open_met_data)' )
1273
1273
END IF
1274
1274
ELSE
1275
- ! block below inserted by rk4417 - phase2
1276
- ok = NF90_INQ_VARID(ncid_met,' dlwrf' ,id% LWdown) ! MMY ! For Princeton
1277
- IF (ok == NF90_NOERR) THEN ! If inquiry is okay ! MMY
1278
- exists% LWdown = .TRUE. ! LWdown is present in met file ! MMY
1279
- ! Get LWdown units and check okay:
1280
- ok = NF90_GET_ATT(ncid_met,id% LWdown,' units' ,metunits% LWdown) ! MMY
1281
- IF (ok /= NF90_NOERR) CALL nc_abort & ! MMY
1282
- (ok,' Error finding LWdown units in met data file ' & ! MMY
1283
- // TRIM (filename% met)// ' (SUBROUTINE open_met_file)' ) ! MMY
1284
- IF (.NOT. (metunits% LWdown(1 :4 )/= ' W/m2' .OR. metunits% LWdown(1 :5 ) & ! MMY
1285
- /= ' W/m^2' .OR. metunits% LWdown(1 :5 )/= ' Wm^-2' & ! MMY
1286
- .OR. metunits% LWdown(1 :4 )/= ' Wm-2' .or. metunits% SWdown(1 :5 ) /= ' W m-2' )) THEN ! MMY
1287
-
1288
- WRITE (* ,* ) metunits% LWdown ! MMY
1289
- CALL abort(' Unknown units for LWdown' // & ! MMY
1290
- ' in ' // TRIM (filename% met)// ' (SUBROUTINE open_met_data)' ) ! MMY
1291
- END IF ! MMY
1292
- ELSE ! MMY
1293
- ! end of block - rk4417 - phase2
1294
1275
exists% LWdown = .FALSE. ! LWdown is not present in met file
1295
1276
all_met= .FALSE. ! not all met variables are present in file
1296
1277
! Note this in log file:
1297
1278
WRITE (logn,* ) ' LWdown not present in met file; ' , &
1298
1279
' values will be synthesised based on air temperature.'
1299
1280
END IF
1300
- END IF ! inserted by rk4417 - phase2
1301
1281
1302
1282
! Look for PSurf (can be synthesised):- - - - - - - - - - - - - - - -
1303
1283
IF (ncciy > 0 ) ncid_met = ncid_ps
0 commit comments