Ignore:
Timestamp:
2017-09-08T18:19:17+02:00 (3 years ago)
Author:
clem
Message:

changes in style - part5 - very nearly finished

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/iceistate.F90

    r8514 r8515  
    3636   PRIVATE 
    3737 
    38    PUBLIC   ice_istate      ! called by icestp.F90 
     38   PUBLIC   ice_istate        ! called by icestp.F90 
     39   PUBLIC   ice_istate_init   ! called by icestp.F90 
    3940 
    4041   INTEGER , PARAMETER ::   jpfldi = 6           ! maximum number of files to read 
     
    102103 
    103104      IF(lwp) WRITE(numout,*) 
    104       IF(lwp) WRITE(numout,*) 'ice_istate : sea-ice initialization ' 
    105       IF(lwp) WRITE(numout,*) '~~~~~~~~~~ ' 
     105      IF(lwp) WRITE(numout,*) 'ice_istate: sea-ice initialization ' 
     106      IF(lwp) WRITE(numout,*) '~~~~~~~~~~' 
    106107 
    107108      !-------------------------------------------------------------------- 
     
    109110      !-------------------------------------------------------------------- 
    110111      ! 
    111       CALL ice_istate_init 
    112  
    113112      ! init surface temperature 
    114113      DO jl = 1, jpl 
     
    154153            !----------------------------- 
    155154            ! assign initial thickness, concentration, snow depth and salinity to an hemisphere-dependent array 
    156             DO jj = 1, jpj 
    157                DO ji = 1, jpi 
    158                   IF( ff_t(ji,jj) >= 0._wp ) THEN 
    159                      zht_i_ini(ji,jj) = rn_hti_ini_n * zswitch(ji,jj) 
    160                      zht_s_ini(ji,jj) = rn_hts_ini_n * zswitch(ji,jj) 
    161                      zat_i_ini(ji,jj) = rn_ati_ini_n * zswitch(ji,jj) 
    162                      zts_u_ini(ji,jj) = rn_tmi_ini_n * zswitch(ji,jj) 
    163                      zsm_i_ini(ji,jj) = rn_smi_ini_n * zswitch(ji,jj) 
    164                      ztm_i_ini(ji,jj) = rn_tmi_ini_n * zswitch(ji,jj) 
    165                   ELSE 
    166                      zht_i_ini(ji,jj) = rn_hti_ini_s * zswitch(ji,jj) 
    167                      zht_s_ini(ji,jj) = rn_hts_ini_s * zswitch(ji,jj) 
    168                      zat_i_ini(ji,jj) = rn_ati_ini_s * zswitch(ji,jj) 
    169                      zts_u_ini(ji,jj) = rn_tmi_ini_s * zswitch(ji,jj) 
    170                      zsm_i_ini(ji,jj) = rn_smi_ini_s * zswitch(ji,jj) 
    171                      ztm_i_ini(ji,jj) = rn_tmi_ini_s * zswitch(ji,jj) 
    172                   ENDIF 
    173                END DO 
    174             END DO 
     155            WHERE( ff_t(:,:) >= 0._wp ) 
     156               zht_i_ini(:,:) = rn_hti_ini_n * zswitch(:,:) 
     157               zht_s_ini(:,:) = rn_hts_ini_n * zswitch(:,:) 
     158               zat_i_ini(:,:) = rn_ati_ini_n * zswitch(:,:) 
     159               zts_u_ini(:,:) = rn_tmi_ini_n * zswitch(:,:) 
     160               zsm_i_ini(:,:) = rn_smi_ini_n * zswitch(:,:) 
     161               ztm_i_ini(:,:) = rn_tmi_ini_n * zswitch(:,:) 
     162            ELSEWHERE 
     163               zht_i_ini(:,:) = rn_hti_ini_s * zswitch(:,:) 
     164               zht_s_ini(:,:) = rn_hts_ini_s * zswitch(:,:) 
     165               zat_i_ini(:,:) = rn_ati_ini_s * zswitch(:,:) 
     166               zts_u_ini(:,:) = rn_tmi_ini_s * zswitch(:,:) 
     167               zsm_i_ini(:,:) = rn_smi_ini_s * zswitch(:,:) 
     168               ztm_i_ini(:,:) = rn_tmi_ini_s * zswitch(:,:) 
     169            END WHERE 
    175170            ! 
    176171         ENDIF ! ln_iceini_file 
     
    554549      slf_i(jp_ati) = sn_ati  ;  slf_i(jp_tsu) = sn_tsu 
    555550      slf_i(jp_tmi) = sn_tmi  ;  slf_i(jp_smi) = sn_smi 
    556  
    557       ! Define the initial parameters 
    558       ! ------------------------- 
    559  
    560       IF(lwp) THEN 
     551      ! 
     552      ! 
     553      IF(lwp) THEN                          ! control print 
    561554         WRITE(numout,*) 
    562555         WRITE(numout,*) 'ice_istate_init: ice parameters inititialisation ' 
    563556         WRITE(numout,*) '~~~~~~~~~~~~~~~' 
    564557         WRITE(numout,*) '   Namelist namice_ini' 
    565          WRITE(numout,*) '      initialization with ice (T) or not (F)       ln_iceini     = ', ln_iceini 
    566          WRITE(numout,*) '      ice initialization from a netcdf file      ln_iceini_file  = ', ln_iceini_file 
    567          WRITE(numout,*) '      threshold water temp. for initial sea-ice    rn_thres_sst  = ', rn_thres_sst 
    568          WRITE(numout,*) '      initial snow thickness in the north          rn_hts_ini_n  = ', rn_hts_ini_n 
    569          WRITE(numout,*) '      initial snow thickness in the south          rn_hts_ini_s  = ', rn_hts_ini_s  
    570          WRITE(numout,*) '      initial ice thickness  in the north          rn_hti_ini_n  = ', rn_hti_ini_n 
    571          WRITE(numout,*) '      initial ice thickness  in the south          rn_hti_ini_s  = ', rn_hti_ini_s 
    572          WRITE(numout,*) '      initial ice concentr.  in the north          rn_ati_ini_n  = ', rn_ati_ini_n 
    573          WRITE(numout,*) '      initial ice concentr.  in the north          rn_ati_ini_s  = ', rn_ati_ini_s 
    574          WRITE(numout,*) '      initial  ice salinity  in the north          rn_smi_ini_n  = ', rn_smi_ini_n 
    575          WRITE(numout,*) '      initial  ice salinity  in the south          rn_smi_ini_s  = ', rn_smi_ini_s 
    576          WRITE(numout,*) '      initial  ice/snw temp  in the north          rn_tmi_ini_n  = ', rn_tmi_ini_n 
    577          WRITE(numout,*) '      initial  ice/snw temp  in the south          rn_tmi_ini_s  = ', rn_tmi_ini_s 
     558         WRITE(numout,*) '      initialization with ice (T) or not (F)                 ln_iceini     = ', ln_iceini 
     559         WRITE(numout,*) '      ice initialization from a netcdf file                ln_iceini_file  = ', ln_iceini_file 
     560         WRITE(numout,*) '      max delta ocean temp. above Tfreeze with initial ice   rn_thres_sst  = ', rn_thres_sst 
     561         WRITE(numout,*) '      initial snow thickness in the north                    rn_hts_ini_n  = ', rn_hts_ini_n 
     562         WRITE(numout,*) '      initial snow thickness in the south                    rn_hts_ini_s  = ', rn_hts_ini_s  
     563         WRITE(numout,*) '      initial ice thickness  in the north                    rn_hti_ini_n  = ', rn_hti_ini_n 
     564         WRITE(numout,*) '      initial ice thickness  in the south                    rn_hti_ini_s  = ', rn_hti_ini_s 
     565         WRITE(numout,*) '      initial ice concentr.  in the north                    rn_ati_ini_n  = ', rn_ati_ini_n 
     566         WRITE(numout,*) '      initial ice concentr.  in the north                    rn_ati_ini_s  = ', rn_ati_ini_s 
     567         WRITE(numout,*) '      initial  ice salinity  in the north                    rn_smi_ini_n  = ', rn_smi_ini_n 
     568         WRITE(numout,*) '      initial  ice salinity  in the south                    rn_smi_ini_s  = ', rn_smi_ini_s 
     569         WRITE(numout,*) '      initial  ice/snw temp  in the north                    rn_tmi_ini_n  = ', rn_tmi_ini_n 
     570         WRITE(numout,*) '      initial  ice/snw temp  in the south                    rn_tmi_ini_s  = ', rn_tmi_ini_s 
    578571      ENDIF 
    579572 
Note: See TracChangeset for help on using the changeset viewer.