- Timestamp:
- 2020-06-09T18:45:11+02:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/UKMO/NEMO_4.0.1_fix_cpl_retain_melt/src/ICE/iceistate.F90
r11715 r13080 41 41 ! !! ** namelist (namini) ** 42 42 LOGICAL, PUBLIC :: ln_iceini !: Ice initialization or not 43 LOGICAL, PUBLIC :: ln_iceini_file !: Ice initialization from 2D netcdf file 43 INTEGER, PUBLIC :: nn_iceini_file ! Ice initialization: 44 ! 0 = Initialise sea ice based on SSTs 45 ! 1 = Initialise sea ice from single category netcdf file 46 ! 2 = Initialise sea ice from multi category restart file 44 47 REAL(wp) :: rn_thres_sst 45 48 REAL(wp) :: rn_hti_ini_n, rn_hts_ini_n, rn_ati_ini_n, rn_smi_ini_n, rn_tmi_ini_n, rn_tsu_ini_n, rn_tms_ini_n … … 48 51 REAL(wp) :: rn_apd_ini_s, rn_hpd_ini_s 49 52 ! 50 ! ! if ln_iceini_file = T53 ! ! if nn_iceini_file = 1 51 54 INTEGER , PARAMETER :: jpfldi = 9 ! maximum number of files to read 52 55 INTEGER , PARAMETER :: jp_hti = 1 ! index of ice thickness (m) … … 155 158 a_ip (:,:,:) = 0._wp 156 159 v_ip (:,:,:) = 0._wp 160 lh_ip (:,:,:) = 0._wp 157 161 a_ip_frac(:,:,:) = 0._wp 162 a_ip_eff (:,:,:) = 0._wp 158 163 h_ip (:,:,:) = 0._wp 159 164 ! … … 167 172 IF( ln_iceini ) THEN 168 173 ! !---------------! 169 IF( ln_iceini_file)THEN ! Read a file !174 IF( nn_iceini_file == 1 )THEN ! Read a file ! 170 175 ! !---------------! 171 176 WHERE( ff_t(:,:) >= 0._wp ) ; zswitch(:,:) = 1._wp … … 362 367 a_ip_frac(:,:,:) = 0._wp 363 368 END WHERE 369 a_ip_eff(:,:,:) = a_ip_frac(:,:,:) 364 370 v_ip(:,:,:) = h_ip(:,:,:) * a_ip(:,:,:) 365 371 … … 466 472 TYPE(FLD_N), DIMENSION(jpfldi) :: slf_i ! array of namelist informations on the fields to read 467 473 ! 468 NAMELIST/namini/ ln_iceini, ln_iceini_file, rn_thres_sst, &474 NAMELIST/namini/ ln_iceini, nn_iceini_file, rn_thres_sst, & 469 475 & rn_hti_ini_n, rn_hti_ini_s, rn_hts_ini_n, rn_hts_ini_s, & 470 476 & rn_ati_ini_n, rn_ati_ini_s, rn_smi_ini_n, rn_smi_ini_s, & … … 493 499 WRITE(numout,*) ' Namelist namini:' 494 500 WRITE(numout,*) ' ice initialization (T) or not (F) ln_iceini = ', ln_iceini 495 WRITE(numout,*) ' ice initialization from a netcdf file ln_iceini_file = ', ln_iceini_file501 WRITE(numout,*) ' ice initialization from a netcdf file nn_iceini_file = ', nn_iceini_file 496 502 WRITE(numout,*) ' max ocean temp. above Tfreeze with initial ice rn_thres_sst = ', rn_thres_sst 497 IF( ln_iceini .AND. .NOT.ln_iceini_file) THEN503 IF( ln_iceini .AND. nn_iceini_file == 0 ) THEN 498 504 WRITE(numout,*) ' initial snw thickness in the north-south rn_hts_ini = ', rn_hts_ini_n,rn_hts_ini_s 499 505 WRITE(numout,*) ' initial ice thickness in the north-south rn_hti_ini = ', rn_hti_ini_n,rn_hti_ini_s … … 508 514 ENDIF 509 515 ! 510 IF( ln_iceini_file) THEN ! Ice initialization using input file516 IF( nn_iceini_file == 1 ) THEN ! Ice initialization using input file 511 517 ! 512 518 ! set si structure
Note: See TracChangeset
for help on using the changeset viewer.