- Timestamp:
- 2019-12-05T17:54:17+01:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_ASINTER-01-05_merged/src/OCE/SBC/sbcblk.F90
r12063 r12074 392 392 !!---------------------------------------------------------------------- 393 393 INTEGER, INTENT(in) :: kt ! ocean time step 394 !! 394 !!---------------------------------------------------------------------- 395 395 REAL(wp), DIMENSION(jpi,jpj) :: zssq, zcd_du, zsen, zevp 396 !!--------------------------------------------------------------------- 396 REAL(wp) :: ztmp 397 !!---------------------------------------------------------------------- 397 398 ! 398 399 CALL fld_read( kt, nn_fsbc, sf ) ! input fields provided at the current time-step 399 ! 400 401 ! Sanity/consistence test on humidity at first time step to detect potential screw-up: 402 IF( kt == nit000 ) THEN 403 WRITE(numout,*) '' 404 ztmp = SUM(tmask(:,:,1)) ! number of ocean points on local proc domain 405 IF( ztmp > 8._wp ) THEN ! test only on proc domains with at least 8 ocean points! 406 ztmp = SUM(sf(jp_humi)%fnow(:,:,1)*tmask(:,:,1))/ztmp ! mean humidity over ocean on proc 407 SELECT CASE( nhumi ) 408 CASE( np_humi_sph ) ! specific humidity => expect: 0. <= something < 0.065 [kg/kg] (0.061 is saturation at 45degC !!!) 409 IF( (ztmp < 0._wp) .OR. (ztmp > 0.065) ) ztmp = -1._wp 410 CASE( np_humi_dpt ) ! dew-point temperature => expect: 110. <= something < 320. [K] 411 IF( (ztmp < 110._wp).OR.(ztmp > 320._wp) ) ztmp = -1._wp 412 CASE( np_humi_rlh ) ! relative humidity => expect: 0. <= something < 100. [%] 413 IF( (ztmp < 0._wp) .OR.(ztmp > 100._wp) ) ztmp = -1._wp 414 END SELECT 415 IF(ztmp < 0._wp) THEN 416 WRITE(numout,'(" Mean humidity value found on proc #",i5.5," is: ",f)'), narea, ztmp 417 CALL ctl_stop( 'STOP', 'Something is wrong with air humidity!!!', & 418 & ' ==> check the unit in your input files' , & 419 & ' ==> check consistence of namelist choice: specific? relative? dew-point?', & 420 & ' ==> ln_humi_sph -> [kg/kg] | ln_humi_rlh -> [%] | ln_humi_dpt -> [K] !!!' ) 421 END IF 422 END IF 423 WRITE(numout,*) ' === Sanity/consistence test on air humidity sucessfuly passed! ===' 424 WRITE(numout,*) '' 425 END IF !IF( kt == nit000 ) 426 400 427 ! ! compute the surface ocean fluxes using bulk formulea 401 428 IF( MOD( kt - 1, nn_fsbc ) == 0 ) THEN
Note: See TracChangeset
for help on using the changeset viewer.