Changeset 4990 for trunk/NEMOGCM/NEMO/LIM_SRC_3/limistate.F90
- Timestamp:
- 2014-12-15T17:42:49+01:00 (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMOGCM/NEMO/LIM_SRC_3/limistate.F90
r4908 r4990 6 6 !! History : 2.0 ! 2004-01 (C. Ethe, G. Madec) Original code 7 7 !! 4.0 ! 2011-02 (G. Madec) dynamical allocation 8 !! - ! 201 2 (C. Rousset) add par_oce (for jp_sal)...bug?8 !! - ! 2014 (C. Rousset) add N/S initializations 9 9 !!---------------------------------------------------------------------- 10 10 #if defined key_lim3 … … 29 29 USE lib_fortran ! Fortran utilities (allows no signed zero when 'key_nosignedzero' defined) 30 30 USE wrk_nemo ! work arrays 31 USE cpl_oasis3, ONLY : lk_cpl32 31 33 32 IMPLICIT NONE … … 36 35 PUBLIC lim_istate ! routine called by lim_init.F90 37 36 38 !! * Module variables39 37 ! !!** init namelist (namiceini) ** 40 38 REAL(wp) :: thres_sst ! threshold water temperature for initial sea ice … … 56 54 !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt) 57 55 !!---------------------------------------------------------------------- 58 59 56 CONTAINS 60 57 … … 80 77 !! 81 78 !! ** Notes : o_i, t_su, t_s, t_i, s_i must be filled everywhere, even 82 !! where there is no ice (clem: I do not know why but it is mandatory)79 !! where there is no ice (clem: I do not know why, is it mandatory?) 83 80 !! 84 81 !! History : … … 116 113 CALL lim_istate_init ! reading the initials parameters of the ice 117 114 118 # if defined key_coupled119 albege(:,:) = 0.8 * tms(:,:)120 # endif121 122 115 ! surface temperature 123 116 DO jl = 1, jpl ! loop over categories … … 125 118 tn_ice(:,:,jl) = rtt * tms(:,:) 126 119 END DO 127 ! Basal temperature is set to the freezing point of seawater in Kelvin 128 t_bo(:,:) = ( tfreez( tsn(:,:,1,jp_sal) ) + rt0 ) * tms(:,:) 120 121 ! basal temperature (considered at freezing point) 122 t_bo(:,:) = ( eos_fzp( tsn(:,:,1,jp_sal) ) + rt0 ) * tms(:,:) 129 123 130 124 IF( ln_limini ) THEN … … 133 127 ! 2) Basal temperature, ice mask and hemispheric index 134 128 !-------------------------------------------------------------------- 135 ! ice if sst <= t-freez + thres_sst 136 DO jj = 1, jpj 129 130 DO jj = 1, jpj ! ice if sst <= t-freez + ttest 137 131 DO ji = 1, jpi 138 132 IF( ( tsn(ji,jj,1,jp_tem) - ( t_bo(ji,jj) - rt0 ) ) * tms(ji,jj) >= thres_sst ) THEN … … 146 140 147 141 ! Hemispheric index 148 ! MV 2011 new initialization149 142 DO jj = 1, jpj 150 143 DO ji = 1, jpi … … 156 149 END DO 157 150 END DO 158 ! END MV 2011 new initialization159 151 160 152 !--------------------------------------------------------------------
Note: See TracChangeset
for help on using the changeset viewer.