- Timestamp:
- 2020-04-25T14:32:42+02:00 (4 years ago)
- Location:
- NEMO/branches/2020/r4.0-HEAD_r12713_clem_dan_fixcpl/src/ICE
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/r4.0-HEAD_r12713_clem_dan_fixcpl/src/ICE/icedyn_adv_pra.F90
r12744 r12816 870 870 ! 871 871 IF( ln_pnd_H12 ) THEN ! melt pond fraction 872 CALL iom_get( numrir, jpdom_autoglo, 'sxap' , sxap ) 873 CALL iom_get( numrir, jpdom_autoglo, 'syap' , syap ) 874 CALL iom_get( numrir, jpdom_autoglo, 'sxxap', sxxap ) 875 CALL iom_get( numrir, jpdom_autoglo, 'syyap', syyap ) 876 CALL iom_get( numrir, jpdom_autoglo, 'sxyap', sxyap ) 877 ! ! melt pond volume 878 CALL iom_get( numrir, jpdom_autoglo, 'sxvp' , sxvp ) 879 CALL iom_get( numrir, jpdom_autoglo, 'syvp' , syvp ) 880 CALL iom_get( numrir, jpdom_autoglo, 'sxxvp', sxxvp ) 881 CALL iom_get( numrir, jpdom_autoglo, 'syyvp', syyvp ) 882 CALL iom_get( numrir, jpdom_autoglo, 'sxyvp', sxyvp ) 883 ! 872 IF( iom_varid( numror, 'sxap', ldstop = .FALSE. ) > 0 ) THEN 873 CALL iom_get( numrir, jpdom_autoglo, 'sxap' , sxap ) 874 CALL iom_get( numrir, jpdom_autoglo, 'syap' , syap ) 875 CALL iom_get( numrir, jpdom_autoglo, 'sxxap', sxxap ) 876 CALL iom_get( numrir, jpdom_autoglo, 'syyap', syyap ) 877 CALL iom_get( numrir, jpdom_autoglo, 'sxyap', sxyap ) 878 ! ! melt pond volume 879 CALL iom_get( numrir, jpdom_autoglo, 'sxvp' , sxvp ) 880 CALL iom_get( numrir, jpdom_autoglo, 'syvp' , syvp ) 881 CALL iom_get( numrir, jpdom_autoglo, 'sxxvp', sxxvp ) 882 CALL iom_get( numrir, jpdom_autoglo, 'syyvp', syyvp ) 883 CALL iom_get( numrir, jpdom_autoglo, 'sxyvp', sxyvp ) 884 ELSE 885 sxap = 0._wp ; syap = 0._wp ; sxxap = 0._wp ; syyap = 0._wp ; sxyap = 0._wp ! melt pond fraction 886 sxvp = 0._wp ; syvp = 0._wp ; sxxvp = 0._wp ; syyvp = 0._wp ; sxyvp = 0._wp ! melt pond volume 887 ENDIF 888 ! 884 889 IF ( ln_pnd_lids ) THEN ! melt pond lid volume 885 CALL iom_get( numrir, jpdom_autoglo, 'sxvl' , sxvl ) 886 CALL iom_get( numrir, jpdom_autoglo, 'syvl' , syvl ) 887 CALL iom_get( numrir, jpdom_autoglo, 'sxxvl', sxxvl ) 888 CALL iom_get( numrir, jpdom_autoglo, 'syyvl', syyvl ) 889 CALL iom_get( numrir, jpdom_autoglo, 'sxyvl', sxyvl ) 890 IF( iom_varid( numror, 'sxvl', ldstop = .FALSE. ) > 0 ) THEN 891 CALL iom_get( numrir, jpdom_autoglo, 'sxvl' , sxvl ) 892 CALL iom_get( numrir, jpdom_autoglo, 'syvl' , syvl ) 893 CALL iom_get( numrir, jpdom_autoglo, 'sxxvl', sxxvl ) 894 CALL iom_get( numrir, jpdom_autoglo, 'syyvl', syyvl ) 895 CALL iom_get( numrir, jpdom_autoglo, 'sxyvl', sxyvl ) 896 ELSE 897 sxvl = 0._wp; syvl = 0._wp ; sxxvl = 0._wp ; syyvl = 0._wp ; sxyvl = 0._wp ! melt pond lid volume 898 ENDIF 890 899 ENDIF 891 900 ENDIF -
NEMO/branches/2020/r4.0-HEAD_r12713_clem_dan_fixcpl/src/ICE/iceistate.F90
r12744 r12816 182 182 ! 183 183 ! -- mandatory fields -- ! 184 zht_i_ini(:,:) = si(jp_hti)%fnow(:,:,1) 185 zht_s_ini(:,:) = si(jp_hts)%fnow(:,:,1) 186 zat_i_ini(:,:) = si(jp_ati)%fnow(:,:,1) 184 zht_i_ini(:,:) = si(jp_hti)%fnow(:,:,1) * tmask(:,:,1) 185 zht_s_ini(:,:) = si(jp_hts)%fnow(:,:,1) * tmask(:,:,1) 186 zat_i_ini(:,:) = si(jp_ati)%fnow(:,:,1) * tmask(:,:,1) 187 187 188 188 ! -- optional fields -- ! … … 226 226 & si(jp_hld)%fnow(:,:,1) = ( rn_hld_ini_n * zswitch + rn_hld_ini_s * (1._wp - zswitch) ) * tmask(:,:,1) 227 227 ! 228 zsm_i_ini(:,:) = si(jp_smi)%fnow(:,:,1) 229 ztm_i_ini(:,:) = si(jp_tmi)%fnow(:,:,1) 230 zt_su_ini(:,:) = si(jp_tsu)%fnow(:,:,1) 231 ztm_s_ini(:,:) = si(jp_tms)%fnow(:,:,1) 232 zapnd_ini(:,:) = si(jp_apd)%fnow(:,:,1) 233 zhpnd_ini(:,:) = si(jp_hpd)%fnow(:,:,1) 234 zhlid_ini(:,:) = si(jp_hld)%fnow(:,:,1) 228 zsm_i_ini(:,:) = si(jp_smi)%fnow(:,:,1) * tmask(:,:,1) 229 ztm_i_ini(:,:) = si(jp_tmi)%fnow(:,:,1) * tmask(:,:,1) 230 zt_su_ini(:,:) = si(jp_tsu)%fnow(:,:,1) * tmask(:,:,1) 231 ztm_s_ini(:,:) = si(jp_tms)%fnow(:,:,1) * tmask(:,:,1) 232 zapnd_ini(:,:) = si(jp_apd)%fnow(:,:,1) * tmask(:,:,1) 233 zhpnd_ini(:,:) = si(jp_hpd)%fnow(:,:,1) * tmask(:,:,1) 234 zhlid_ini(:,:) = si(jp_hld)%fnow(:,:,1) * tmask(:,:,1) 235 235 ! 236 236 ! change the switch for the following
Note: See TracChangeset
for help on using the changeset viewer.