New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 12024 for NEMO/branches/2019/dev_r11085_ASINTER-05_Brodeau_Advanced_Bulk/src – NEMO

Ignore:
Timestamp:
2019-12-02T12:33:48+01:00 (4 years ago)
Author:
laurent
Message:

Small fixes...

File:
1 edited

Legend:

Unmodified
Added
Removed

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

    r12020 r12024  
    641641#endif 
    642642      ! 
    643       !!#LB: NO WHY???? IF( nn_ice == 0 ) THEN 
    644643      CALL iom_put( "rho_air"  ,   rhoa )                 ! output air density (kg/m^3) !#LB 
    645       CALL iom_put( "qlw_oce"  ,   zqlw )                 ! output downward longwave heat over the ocean 
    646       CALL iom_put( "qsb_oce"  ,   zqsb )                 ! output downward sensible heat over the ocean 
    647       CALL iom_put( "qla_oce"  ,   zqla )                 ! output downward latent   heat over the ocean 
    648644      CALL iom_put( "evap_oce" ,  zevap )                 ! evaporation 
    649       CALL iom_put( "qemp_oce" ,   qns-zqlw-zqsb-zqla )   ! output downward heat content of E-P over the ocean 
    650       CALL iom_put( "qns_oce"  ,   qns  )                 ! output downward non solar heat over the ocean 
    651       CALL iom_put( "qsr_oce"  ,   qsr  )                 ! output downward solar heat over the ocean 
    652       CALL iom_put( "qt_oce"   ,   qns+qsr )              ! output total downward heat over the ocean 
    653       tprecip(:,:) = sf(jp_prec)%fnow(:,:,1) * rn_pfac * tmask(:,:,1) ! output total precipitation [kg/m2/s] 
    654       sprecip(:,:) = sf(jp_snow)%fnow(:,:,1) * rn_pfac * tmask(:,:,1) ! output solid precipitation [kg/m2/s] 
    655       CALL iom_put( 'snowpre', sprecip )                 ! Snow 
    656       CALL iom_put( 'precip' , tprecip )                 ! Total precipitation 
     645      ! 
     646      IF ( nn_ice == 0 ) THEN 
     647         CALL iom_put( "qlw_oce"  ,   zqlw )                 ! output downward longwave heat over the ocean 
     648         CALL iom_put( "qsb_oce"  ,   zqsb )                 ! output downward sensible heat over the ocean 
     649         CALL iom_put( "qla_oce"  ,   zqla )                 ! output downward latent   heat over the ocean 
     650         CALL iom_put( "qemp_oce" ,   qns-zqlw-zqsb-zqla )   ! output downward heat content of E-P over the ocean 
     651         CALL iom_put( "qns_oce"  ,   qns  )                 ! output downward non solar heat over the ocean 
     652         CALL iom_put( "qsr_oce"  ,   qsr  )                 ! output downward solar heat over the ocean 
     653         CALL iom_put( "qt_oce"   ,   qns+qsr )              ! output total downward heat over the ocean 
     654         tprecip(:,:) = sf(jp_prec)%fnow(:,:,1) * rn_pfac * tmask(:,:,1) ! output total precipitation [kg/m2/s] 
     655         sprecip(:,:) = sf(jp_snow)%fnow(:,:,1) * rn_pfac * tmask(:,:,1) ! output solid precipitation [kg/m2/s] 
     656         CALL iom_put( 'snowpre', sprecip )                 ! Snow 
     657         CALL iom_put( 'precip' , tprecip )                 ! Total precipitation 
     658      ENDIF 
    657659      IF( ln_skin_cs .OR. ln_skin_wl ) THEN 
    658660         CALL iom_put( "t_skin" ,  (zst - rt0) * tmask(:,:,1) )           ! T_skin in Celsius 
Note: See TracChangeset for help on using the changeset viewer.