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 13004 – NEMO

Changeset 13004


Ignore:
Timestamp:
2020-06-02T14:22:16+02:00 (3 years ago)
Author:
cetlod
Message:

r4.0-HEAD_r12713 : bugfix on cloud fraction in coupled mode, see ticket #2473

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2020/r4.0-HEAD_r12713_clem_dan_fixcpl/src/OCE/SBC/sbccpl.F90

    r12955 r13004  
    11291129      REAL(wp) ::   zcdrag = 1.5e-3        ! drag coefficient 
    11301130      REAL(wp) ::   zzx, zzy               ! temporary variables 
    1131       REAL(wp), DIMENSION(jpi,jpj) ::   ztx, zty, zmsk, zemp, zqns, zqsr 
     1131      REAL(wp), DIMENSION(jpi,jpj) ::   ztx, zty, zmsk, zemp, zqns, zqsr, zcloud_fra 
    11321132      !!---------------------------------------------------------------------- 
    11331133      ! 
     
    12471247         ENDIF 
    12481248      ENDIF 
    1249  
     1249!!$      !                                                      ! ========================= ! 
     1250!!$      SELECT CASE( TRIM( sn_rcv_clouds%cldes ) )             !       cloud fraction      ! 
     1251!!$      !                                                      ! ========================= ! 
     1252!!$      cloud_fra(:,:) = frcv(jpr_clfra)*z3(:,:,1) 
     1253!!$      END SELECT 
     1254!!$ 
     1255      zcloud_fra(:,:) = pp_cldf   ! should be real cloud fraction instead (as in the bulk) but needs to be read from atm. 
     1256      IF( ln_mixcpl ) THEN 
     1257         cloud_fra(:,:) = cloud_fra(:,:) * xcplmask(:,:,0) + zcloud_fra(:,:)* zmsk(:,:) 
     1258      ELSE 
     1259         cloud_fra(:,:) = zcloud_fra(:,:) 
     1260      ENDIF 
     1261      !                                                      ! ========================= ! 
    12501262      ! u(v)tau and taum will be modified by ice model 
    12511263      ! -> need to be reset before each call of the ice/fsbc       
     
    16491661      REAL(wp), DIMENSION(jpi,jpj)     ::   zevap_ice_total 
    16501662      REAL(wp), DIMENSION(jpi,jpj,jpl) ::   zqns_ice, zqsr_ice, zdqns_ice, zqevap_ice, zevap_ice, zqtr_ice_top, ztsu 
    1651       REAL(wp), DIMENSION(jpi,jpj)     ::   ztri, zcloud_fra 
     1663      REAL(wp), DIMENSION(jpi,jpj)     ::   ztri 
    16521664      !!---------------------------------------------------------------------- 
    16531665      ! 
     
    20852097         ENDIF 
    20862098      END SELECT 
    2087 !!$      !                                                      ! ========================= ! 
    2088 !!$      SELECT CASE( TRIM( sn_rcv_clouds%cldes ) )             !       cloud fraction      ! 
    2089 !!$      !                                                      ! ========================= ! 
    2090 !!$         cloud_fra(:,:) = frcv(jpr_clfra)*z3(:,:,1) 
    2091 !!$      END SELECT 
    2092       zcloud_fra(:,:) = pp_cldf   ! should be real cloud fraction instead (as in the bulk) but needs to be read from atm. 
    2093       IF( ln_mixcpl ) THEN 
    2094          cloud_fra(:,:) = cloud_fra(:,:) * xcplmask(:,:,0) + zcloud_fra(:,:)* zmsk(:,:) 
    2095       ELSE 
    2096          cloud_fra(:,:) = zcloud_fra(:,:) 
    2097       ENDIF 
    20982099      !                                                      ! ========================= ! 
    20992100      !                                                      !      Transmitted Qsr      !   [W/m2] 
     
    21082109            !    2) is 0 when there is any snow 
    21092110            !    3) tends to 1 for thin ice 
    2110             ztri(:,:) = 0.18 * ( 1.0 - zcloud_fra(:,:) ) + 0.35 * zcloud_fra(:,:)  ! surface transmission when hi>10cm 
     2111            ztri(:,:) = 0.18 * ( 1.0 - cloud_fra(:,:) ) + 0.35 * cloud_fra(:,:)  ! surface transmission when hi>10cm 
    21112112            DO jl = 1, jpl 
    21122113               WHERE    ( phs(:,:,jl) <= 0._wp .AND. phi(:,:,jl) <  0.1_wp )       ! linear decrease from hi=0 to 10cm   
Note: See TracChangeset for help on using the changeset viewer.