Ignore:
Timestamp:
2018-05-10T16:26:43+02:00 (3 years ago)
Author:
mathiot
Message:

apply tmask(:,:,1) to erp, qrp, precip, snow, rnf, qns to fix #2083

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2015/nemo_v3_6_STABLE/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk_core.F90

    r8529 r9567  
    394394         &     * ( sf(jp_tair)%fnow(:,:,1) - rt0 ) * rcp                          & 
    395395         &     + sf(jp_snow)%fnow(:,:,1) * rn_pfac                                &   ! add solid  precip heat content at min(Tair,Tsnow) 
    396          &     * ( MIN( sf(jp_tair)%fnow(:,:,1), rt0_snow ) - rt0 ) * cpic * tmask(:,:,1) 
     396         &     * ( MIN( sf(jp_tair)%fnow(:,:,1), rt0_snow ) - rt0 ) * cpic 
     397      qns(:,:) = qns(:,:) * tmask(:,:,1) 
    397398      ! 
    398399#if defined key_lim3 
     
    409410         CALL iom_put( "qsr_oce" ,   qsr  )                 ! output downward solar heat over the ocean 
    410411         CALL iom_put( "qt_oce"  ,   qns+qsr )              ! output total downward heat over the ocean 
    411          tprecip(:,:) = sf(jp_prec)%fnow(:,:,1) * rn_pfac   ! output total precipitation [kg/m2/s] 
    412          sprecip(:,:) = sf(jp_snow)%fnow(:,:,1) * rn_pfac   ! output solid precipitation [kg/m2/s] 
     412         tprecip(:,:) = sf(jp_prec)%fnow(:,:,1) * rn_pfac * tmask(:,:,1) ! output total precipitation [kg/m2/s] 
     413         sprecip(:,:) = sf(jp_snow)%fnow(:,:,1) * rn_pfac * tmask(:,:,1) ! output solid precipitation [kg/m2/s] 
    413414         CALL iom_put( 'snowpre', sprecip * 86400. )        ! Snow 
    414415         CALL iom_put( 'precip' , tprecip * 86400. )        ! Total precipitation 
     
    608609      END DO 
    609610      ! 
    610       tprecip(:,:) = sf(jp_prec)%fnow(:,:,1) * rn_pfac      ! total precipitation [kg/m2/s] 
    611       sprecip(:,:) = sf(jp_snow)%fnow(:,:,1) * rn_pfac      ! solid precipitation [kg/m2/s] 
     611      tprecip(:,:) = sf(jp_prec)%fnow(:,:,1) * rn_pfac * tmask(:,:,1)  ! total precipitation [kg/m2/s] 
     612      sprecip(:,:) = sf(jp_snow)%fnow(:,:,1) * rn_pfac * tmask(:,:,1)  ! solid precipitation [kg/m2/s] 
    612613      CALL iom_put( 'snowpre', sprecip * 86400. )                  ! Snow precipitation 
    613614      CALL iom_put( 'precip' , tprecip * 86400. )                  ! Total precipitation 
Note: See TracChangeset for help on using the changeset viewer.