Ignore:
Timestamp:
2020-04-24T17:20:27+02:00 (12 months ago)
Message:

debug a restartability issue. All sette tests passed now

Location:
Files:
8 edited

### Legend:

Unmodified
Removed

 r12726 h_ip_1d(ji)      = 0._wp h_il_1d(ji)      = 0._wp ! ! clem: problem with conservation or not ? !                                                         !--------------------------------! ELSE                                                         ! Case ice thickness >= rn_himin ! zhp = ( h_i_1d(ji) * ( rau0 - rhoi ) + h_ip_1d(ji) * ( rau0 - rhow * a_ip_1d(ji) / a_i_1d(ji) ) ) * r1_rau0 ! Calculate the permeability of the ice (Assur 1958) ! Calculate the permeability of the ice (Assur 1958, see Flocco 2010) DO jk = 1, nlay_i zsbr = - 1.2_wp                                  & &   - 21.8_wp    * ( t_i_1d(ji,jk) - rt0 )    & &   - 0.919_wp   * ( t_i_1d(ji,jk) - rt0 )**2 & &   - 0.0178_wp  * ( t_i_1d(ji,jk) - rt0 )**3 ! clem: error here the factor was 0.01878 instead of 0.0178 (cf Flocco 2010) &   - 0.0178_wp  * ( t_i_1d(ji,jk) - rt0 )**3 ztmp(jk) = sz_i_1d(ji,jk) / zsbr END DO

 r12785 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:)   ::   wndm_ice       !: wind speed module at T-point                 [m/s] REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:)   ::   sstfrz         !: sea surface freezing temperature            [degC] REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:)   ::   cloud_fra      !: cloud cover                                    [-] #endif #endif REAL(wp), PUBLIC, SAVE ::   cldf_ice = 0.81    !: cloud fraction over sea ice, summer CLIO value   [-] REAL(wp), PUBLIC, SAVE ::   pp_cldf = 0.81    !: cloud fraction over sea ice, summer CLIO value   [-] !! arrays relating to embedding ice in the ocean &      qemp_ice(jpi,jpj)     , qevap_ice(jpi,jpj,jpl) , qemp_oce   (jpi,jpj)     ,   & &      qns_oce (jpi,jpj)     , qsr_oce  (jpi,jpj)     , emp_oce    (jpi,jpj)     ,   & &      emp_ice (jpi,jpj)     , sstfrz   (jpi,jpj)     , cloud_fra  (jpi,jpj)     , STAT= ierr(2) ) &      emp_ice (jpi,jpj)     , sstfrz   (jpi,jpj)     , STAT= ierr(2) ) #endif LOGICAL         , PUBLIC, PARAMETER ::   lk_si3     = .FALSE.  !: no SI3 ice model LOGICAL         , PUBLIC, PARAMETER ::   lk_cice    = .FALSE.  !: no CICE ice model REAL(wp)        , PUBLIC, PARAMETER ::   cldf_ice = 0.81       !: cloud fraction over sea ice, summer CLIO value   [-] REAL(wp)        , PUBLIC, PARAMETER ::   pp_cldf    = 0.81     !: cloud fraction over sea ice, summer CLIO value   [-] INTEGER         , PUBLIC, PARAMETER ::   jpl = 1 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:)   ::   u_ice, v_ice                        ! jpi, jpj

 r12132 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   atm_co2           !: atmospheric pCO2                             [ppm] REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: xcplmask          !: coupling mask for ln_mixcpl (warning: allocated in sbccpl) REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   cloud_fra         !: cloud cover (fraction of cloud in a gridcell) [-] !!---------------------------------------------------------------------- &      fwficb  (jpi,jpj), fwficb_b(jpi,jpj), STAT=ierr(3) ) ! ALLOCATE( tprecip(jpi,jpj) , sprecip(jpi,jpj) , fr_i(jpi,jpj) ,     & &      atm_co2(jpi,jpj) ,                                        & &      ssu_m  (jpi,jpj) , sst_m(jpi,jpj) , frq_m(jpi,jpj) ,      & &      ssv_m  (jpi,jpj) , sss_m(jpi,jpj) , ssh_m(jpi,jpj) , STAT=ierr(4) ) ALLOCATE( tprecip(jpi,jpj) , sprecip  (jpi,jpj) , fr_i(jpi,jpj) ,   & &      atm_co2(jpi,jpj) , cloud_fra(jpi,jpj) ,                   & &      ssu_m  (jpi,jpj) , sst_m    (jpi,jpj) , frq_m(jpi,jpj) ,  & &      ssv_m  (jpi,jpj) , sss_m    (jpi,jpj) , ssh_m(jpi,jpj) , STAT=ierr(4) ) ! ALLOCATE( e3t_m(jpi,jpj) , STAT=ierr(5) )

 r12785 ENDDO ! fill cloud cover array with constant value if "not used" IF( TRIM(sf(jp_cc)%clrootname) == 'NOT USED' )   sf(jp_cc)%fnow(:,:,1) = cldf_ice IF( TRIM(sf(jp_cc)%clrootname) == 'NOT USED' )   sf(jp_cc)%fnow(:,:,1) = pp_cldf IF ( ln_wave ) THEN zst(:,:) = pst(:,:) + rt0      ! convert SST from Celcius to Kelvin (and set minimum value far above 0 K) ! --- cloud cover --- ! cloud_fra(:,:) = sf(jp_cc)%fnow(:,:,1) ! ----------------------------------------------------------------------------- ! !      0   Wind components and module at T-point relative to the moving ocean   ! END DO ! --- cloud cover --- ! cloud_fra(:,:) = sf(jp_cc)%fnow(:,:,1) ! --- shortwave radiation transmitted below the surface (W/m2, see Grenfell Maykut 77) --- ! ztri(:,:) = 0.18 * ( 1.0 - cloud_fra(:,:) ) + 0.35 * cloud_fra(:,:)  ! surface transmission when hi>10cm