- Timestamp:
- 2020-04-09T17:06:57+02:00 (4 years ago)
- Location:
- NEMO/branches/2020/r4.0-HEAD_r12713_clem_dan_fixcpl/src/ICE
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/r4.0-HEAD_r12713_clem_dan_fixcpl/src/ICE/iceistate.F90
r12726 r12728 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, rn_hld_ini_s 49 52 ! 50 ! ! if ln_iceini_file = T53 ! ! if nn_iceini_file = 1 51 54 INTEGER , PARAMETER :: jpfldi = 10 ! maximum number of files to read 52 55 INTEGER , PARAMETER :: jp_hti = 1 ! index of ice thickness (m) … … 170 173 IF( ln_iceini ) THEN 171 174 ! !---------------! 172 IF( ln_iceini_file )THEN! Read a file !175 IF( nn_iceini_file == 1 )THEN ! Read a file ! 173 176 ! !---------------! 174 177 WHERE( ff_t(:,:) >= 0._wp ) ; zswitch(:,:) = 1._wp … … 483 486 TYPE(FLD_N), DIMENSION(jpfldi) :: slf_i ! array of namelist informations on the fields to read 484 487 ! 485 NAMELIST/namini/ ln_iceini, ln_iceini_file, rn_thres_sst, &488 NAMELIST/namini/ ln_iceini, nn_iceini_file, rn_thres_sst, & 486 489 & rn_hti_ini_n, rn_hti_ini_s, rn_hts_ini_n, rn_hts_ini_s, & 487 490 & rn_ati_ini_n, rn_ati_ini_s, rn_smi_ini_n, rn_smi_ini_s, & … … 510 513 WRITE(numout,*) ' Namelist namini:' 511 514 WRITE(numout,*) ' ice initialization (T) or not (F) ln_iceini = ', ln_iceini 512 WRITE(numout,*) ' ice initialization from a netcdf file ln_iceini_file = ', ln_iceini_file515 WRITE(numout,*) ' ice initialization from a netcdf file nn_iceini_file = ', nn_iceini_file 513 516 WRITE(numout,*) ' max ocean temp. above Tfreeze with initial ice rn_thres_sst = ', rn_thres_sst 514 IF( ln_iceini .AND. .NOT.ln_iceini_file) THEN517 IF( ln_iceini .AND. nn_iceini_file == 0 ) THEN 515 518 WRITE(numout,*) ' initial snw thickness in the north-south rn_hts_ini = ', rn_hts_ini_n,rn_hts_ini_s 516 519 WRITE(numout,*) ' initial ice thickness in the north-south rn_hti_ini = ', rn_hti_ini_n,rn_hti_ini_s … … 526 529 ENDIF 527 530 ! 528 IF( ln_iceini_file) THEN ! Ice initialization using input file531 IF( nn_iceini_file == 1 ) THEN ! Ice initialization using input file 529 532 ! 530 533 ! set si structure -
NEMO/branches/2020/r4.0-HEAD_r12713_clem_dan_fixcpl/src/ICE/icerst.F90
r12725 r12728 283 283 CALL ice_istate( nit000 ) 284 284 ! 285 IF( .NOT.ln_iceini .OR. .NOT.ln_iceini_file) &286 & CALL ctl_stop('STOP', 'ice_rst_read: you need ln_ice_ini=T and ln_iceini_file=T')285 IF( .NOT.ln_iceini .OR. nn_iceini_file == 0 ) & 286 & CALL ctl_stop('STOP', 'ice_rst_read: you need ln_ice_ini=T and nn_iceini_file=0') 287 287 ! 288 288 ENDIF -
NEMO/branches/2020/r4.0-HEAD_r12713_clem_dan_fixcpl/src/ICE/icestp.F90
r12726 r12728 252 252 ! 253 253 ! ! Initial sea-ice state 254 IF( .NOT. ln_rstart ) THEN ! start from rest: sea-ice deduced from sst 254 255 IF ( ln_rstart .OR. nn_iceini_file == 2 ) THEN 256 CALL ice_rst_read ! start from a restart file 257 ELSE 255 258 CALL ice_istate_init 256 CALL ice_istate( nit000 ) 257 ELSE ! start from a restart file 258 CALL ice_rst_read 259 CALL ice_istate( nit000 ) ! start from rest: sea-ice deduced from sst 259 260 ENDIF 260 261 CALL ice_var_glo2eqv
Note: See TracChangeset
for help on using the changeset viewer.