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 13080 for NEMO/branches/UKMO/NEMO_4.0.1_fix_cpl_retain_melt/src/ICE/iceistate.F90 – NEMO

Ignore:
Timestamp:
2020-06-09T18:45:11+02:00 (4 years ago)
Author:
dancopsey
Message:

Merge in existing fix_cpl branch.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/UKMO/NEMO_4.0.1_fix_cpl_retain_melt/src/ICE/iceistate.F90

    r11715 r13080  
    4141   !                             !! ** namelist (namini) ** 
    4242   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 
    4447   REAL(wp) ::   rn_thres_sst 
    4548   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 
     
    4851   REAL(wp) ::   rn_apd_ini_s, rn_hpd_ini_s 
    4952   ! 
    50    !                              ! if ln_iceini_file = T 
     53   !                              ! if nn_iceini_file = 1 
    5154   INTEGER , PARAMETER ::   jpfldi = 9           ! maximum number of files to read 
    5255   INTEGER , PARAMETER ::   jp_hti = 1           ! index of ice thickness    (m) 
     
    155158      a_ip     (:,:,:) = 0._wp 
    156159      v_ip     (:,:,:) = 0._wp 
     160      lh_ip    (:,:,:) = 0._wp 
    157161      a_ip_frac(:,:,:) = 0._wp 
     162      a_ip_eff (:,:,:) = 0._wp 
    158163      h_ip     (:,:,:) = 0._wp 
    159164      ! 
     
    167172      IF( ln_iceini ) THEN 
    168173         !                             !---------------! 
    169          IF( ln_iceini_file )THEN      ! Read a file   ! 
     174         IF( nn_iceini_file == 1 )THEN      ! Read a file   ! 
    170175            !                          !---------------! 
    171176            WHERE( ff_t(:,:) >= 0._wp )   ;   zswitch(:,:) = 1._wp 
     
    362367            a_ip_frac(:,:,:) = 0._wp 
    363368         END WHERE 
     369         a_ip_eff(:,:,:) = a_ip_frac(:,:,:) 
    364370         v_ip(:,:,:) = h_ip(:,:,:) * a_ip(:,:,:) 
    365371           
     
    466472      TYPE(FLD_N), DIMENSION(jpfldi) ::   slf_i                 ! array of namelist informations on the fields to read 
    467473      ! 
    468       NAMELIST/namini/ ln_iceini, ln_iceini_file, rn_thres_sst, & 
     474      NAMELIST/namini/ ln_iceini, nn_iceini_file, rn_thres_sst, & 
    469475         &             rn_hti_ini_n, rn_hti_ini_s, rn_hts_ini_n, rn_hts_ini_s, & 
    470476         &             rn_ati_ini_n, rn_ati_ini_s, rn_smi_ini_n, rn_smi_ini_s, & 
     
    493499         WRITE(numout,*) '   Namelist namini:' 
    494500         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_file 
     501         WRITE(numout,*) '      ice initialization from a netcdf file            nn_iceini_file = ', nn_iceini_file 
    496502         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 ) THEN 
     503         IF( ln_iceini .AND. nn_iceini_file == 0 ) THEN 
    498504            WRITE(numout,*) '      initial snw thickness in the north-south         rn_hts_ini     = ', rn_hts_ini_n,rn_hts_ini_s  
    499505            WRITE(numout,*) '      initial ice thickness in the north-south         rn_hti_ini     = ', rn_hti_ini_n,rn_hti_ini_s 
     
    508514      ENDIF 
    509515      ! 
    510       IF( ln_iceini_file ) THEN                      ! Ice initialization using input file 
     516      IF( nn_iceini_file == 1 ) THEN                      ! Ice initialization using input file 
    511517         ! 
    512518         ! set si structure 
Note: See TracChangeset for help on using the changeset viewer.