Changeset 11660


Ignore:
Timestamp:
2019-10-08T12:09:21+02:00 (12 months ago)
Author:
laurent
Message:

LB: xios output of evaporation in sbcblk.F90

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2019/dev_r11085_ASINTER-05_Brodeau_Advanced_Bulk/src/OCE/SBC/sbcblk.F90

    r11638 r11660  
    644644      IF( ABS( rn_zu - rn_zqt) < 0.01_wp ) THEN 
    645645         !! q_air and t_air are given at 10m (wind reference height) 
    646          zevap(:,:) = rn_efac*MAX( 0._wp,             zqla(:,:)*Ce_atm(:,:)*(zsq(:,:) - zqair(:,:)) ) ! Evaporation, using bulk wind speed 
    647          zqsb (:,:) = cp_air(zqair(:,:))*zqla(:,:)*Ch_atm(:,:)*(zst(:,:) - ztpot(:,:)             )  ! Sensible Heat, using bulk wind speed 
     646         zevap(:,:) = rn_efac*MAX( 0._wp,  zqla(:,:)*Ce_atm(:,:)*(zsq(:,:) - zqair(:,:)) ) ! Evaporation, using bulk wind speed 
     647         zqsb (:,:) = cp_air(zqair(:,:)) * zqla(:,:)*Ch_atm(:,:)*(zst(:,:) - ztpot(:,:)  ) ! Sensible Heat, using bulk wind speed 
    648648      ELSE 
    649649         !! q_air and t_air are not given at 10m (wind reference height) 
    650650         ! Values of temp. and hum. adjusted to height of wind during bulk algorithm iteration must be used!!! 
    651          zevap(:,:) = rn_efac*MAX( 0._wp,             zqla(:,:)*Ce_atm(:,:)*(zsq(:,:) - q_zu(:,:) ) ) ! Evaporation, using bulk wind speed 
    652          zqsb (:,:) = cp_air(zqair(:,:))*zqla(:,:)*Ch_atm(:,:)*(zst(:,:) - t_zu(:,:) )   ! Sensible Heat, using bulk wind speed 
     651         zevap(:,:) = rn_efac*MAX( 0._wp,  zqla(:,:)*Ce_atm(:,:)*(zsq(:,:) - q_zu(:,:)) ) ! Evaporation, using bulk wind speed 
     652         zqsb (:,:) = cp_air(zqair(:,:)) * zqla(:,:)*Ch_atm(:,:)*(zst(:,:) - t_zu(:,:) )   ! Sensible Heat, using bulk wind speed 
    653653      ENDIF 
    654654 
     
    665665 
    666666      IF(ln_ctl) THEN 
     667         CALL prt_ctl( tab2d_1=zevap , clinfo1=' blk_oce: zevap  : ') 
    667668         CALL prt_ctl( tab2d_1=zqla  , clinfo1=' blk_oce: zqla   : ', tab2d_2=Ce_atm , clinfo2=' Ce_oce  : ' ) 
    668669         CALL prt_ctl( tab2d_1=zqsb  , clinfo1=' blk_oce: zqsb   : ', tab2d_2=Ch_atm , clinfo2=' Ch_oce  : ' ) 
     
    697698      ! 
    698699      !!#LB: NO WHY???? IF ( nn_ice == 0 ) THEN 
    699       CALL iom_put( "rho_air" ,   rhoa )                 ! output air density (kg/m^3) !#LB 
    700       CALL iom_put( "qlw_oce" ,   zqlw )                 ! output downward longwave heat over the ocean 
    701       CALL iom_put( "qsb_oce" , - zqsb )                 ! output downward sensible heat over the ocean 
    702       CALL iom_put( "qla_oce" , - zqla )                 ! output downward latent   heat over the ocean 
    703       CALL iom_put( "qemp_oce",   qns-zqlw+zqsb+zqla )   ! output downward heat content of E-P over the ocean 
    704       CALL iom_put( "qns_oce" ,   qns  )                 ! output downward non solar heat over the ocean 
    705       CALL iom_put( "qsr_oce" ,   qsr  )                 ! output downward solar heat over the ocean 
    706       CALL iom_put( "qt_oce"  ,   qns+qsr )              ! output total downward heat over the ocean 
     700      CALL iom_put( "rho_air"  ,   rhoa )                 ! output air density (kg/m^3) !#LB 
     701      CALL iom_put( "qlw_oce"  ,   zqlw )                 ! output downward longwave heat over the ocean 
     702      CALL iom_put( "qsb_oce"  , - zqsb )                 ! output downward sensible heat over the ocean 
     703      CALL iom_put( "qla_oce"  , - zqla )                 ! output downward latent   heat over the ocean 
     704      CALL iom_put( "evap_oce" ,  zevap )                 ! evaporation 
     705      CALL iom_put( "qemp_oce" ,   qns-zqlw+zqsb+zqla )   ! output downward heat content of E-P over the ocean 
     706      CALL iom_put( "qns_oce"  ,   qns  )                 ! output downward non solar heat over the ocean 
     707      CALL iom_put( "qsr_oce"  ,   qsr  )                 ! output downward solar heat over the ocean 
     708      CALL iom_put( "qt_oce"   ,   qns+qsr )              ! output total downward heat over the ocean 
    707709      tprecip(:,:) = sf(jp_prec)%fnow(:,:,1) * rn_pfac * tmask(:,:,1) ! output total precipitation [kg/m2/s] 
    708710      sprecip(:,:) = sf(jp_snow)%fnow(:,:,1) * rn_pfac * tmask(:,:,1) ! output solid precipitation [kg/m2/s] 
Note: See TracChangeset for help on using the changeset viewer.