Changeset 921 for trunk/NEMO/LIM_SRC_3/iceini.F90
- Timestamp:
- 2008-05-13T10:28:52+02:00 (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/LIM_SRC_3/iceini.F90
r888 r921 32 32 33 33 !! * Share Module variables 34 LOGICAL , PUBLIC :: & !!! ** init namelist (namicerun) **35 ln_limdyn = .TRUE., & !: flag for ice dynamics (T) or not (F)36 ln_nicep = .TRUE. !: flag for sea-ice points output (T) or not (F)37 34 INTEGER , PUBLIC :: & !: 38 35 nstart , & !: iteration number of the begining of the run … … 41 38 numit !: iteration number 42 39 REAL(wp), PUBLIC :: & !: 43 hsndif = 0.e0 , & !: computation of temp. in snow (0) or not (9999) 44 hicdif = 0.e0 , & !: computation of temp. in ice (0) or not (9999) 45 tpstot , & !: time of the run in seconds 46 cai = 1.40e-3 , & !: atmospheric drag over sea ice 47 cao = 1.00e-3 !: atmospheric drag over ocean 48 REAL(wp), PUBLIC, DIMENSION(2) :: & !: 49 acrit = (/ 1.e-06 , 1.e-06 /) !: minimum fraction for leads in 50 ! ! north and south hemisphere 40 tpstot !: time of the run in seconds 51 41 !!---------------------------------------------------------------------- 52 42 !! LIM 3.0, UCL-ASTR-LOCEAN-IPSL (2008) … … 72 62 73 63 CALL ice_run ! read in namelist some run parameters 74 64 75 65 ! Louvain la Neuve Ice model 76 66 IF( nacc == 1 ) THEN 77 78 67 dtsd2 = nn_fsbc * rdtmin * 0.5 68 rdt_ice = nn_fsbc * rdtmin 79 69 ELSE 80 81 70 dtsd2 = nn_fsbc * rdt * 0.5 71 rdt_ice = nn_fsbc * rdt 82 72 ENDIF 83 73 84 74 CALL lim_msh ! ice mesh initialization 85 75 86 76 CALL lim_itd_ini ! initialize the ice thickness 87 77 ! distribution 88 78 ! Initial sea-ice state 89 79 IF( .NOT.ln_rstart ) THEN … … 92 82 CALL lim_istate ! start from rest: sea-ice deduced from sst 93 83 CALL lim_var_agg(1) ! aggregate category variables in 94 84 ! bulk variables 95 85 CALL lim_var_glo2eqv ! convert global variables in equivalent 96 86 ! variables 97 87 ELSE 98 88 CALL lim_rst_read ! start from a restart file … … 108 98 alb_ice(:,:,:) = albege(:,:) ! sea-ice albedo 109 99 # endif 110 100 111 101 nstart = numit + nn_fsbc 112 102 nitrun = nitend - nit000 + 1 … … 138 128 REWIND ( numnam_ice ) 139 129 READ ( numnam_ice , namicerun ) 130 ln_nicep = ln_nicep .AND. lwp 140 131 IF(lwp) THEN 141 132 WRITE(numout,*) … … 150 141 WRITE(numout,*) ' Several ice points in the ice or not in ocean.output = ', ln_nicep 151 142 ENDIF 152 143 153 144 END SUBROUTINE ice_run 154 145 155 146 SUBROUTINE lim_itd_ini 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 !!-- End of declarations195 !!------------------------------------------------------------------------------196 197 !------------------------------------------------------------------------------!198 ! 1) Ice thickness distribution parameters initialization199 !------------------------------------------------------------------------------!147 !!------------------------------------------------------------------ 148 !! *** ROUTINE lim_itd_ini *** 149 !! ** Purpose : 150 !! Initializes the ice thickness distribution 151 !! ** Method : 152 !! Very simple. Currently there are no ice types in the 153 !! model... 154 !! 155 !! ** Arguments : 156 !! kideb , kiut : Starting and ending points on which the 157 !! the computation is applied 158 !! 159 !! ** Inputs / Ouputs : (global commons) 160 !! 161 !! ** External : 162 !! 163 !! ** References : 164 !! 165 !! ** History : 166 !! (12-2005) Martin Vancoppenolle 167 !! 168 !!------------------------------------------------------------------ 169 !! * Arguments 170 171 !! * Local variables 172 INTEGER :: jl, & ! ice category dummy loop index 173 jm ! ice types dummy loop index 174 175 REAL(wp) :: & ! constant values 176 zeps = 1.0e-10, & ! 177 zc1 , & ! 178 zc2 , & ! 179 zc3 , & ! 180 zx1 181 182 WRITE(numout,*) 'lim_itd_ini : Initialization of ice thickness distribution ' 183 WRITE(numout,*) '~~~~~~~~~~~~' 184 185 !!-- End of declarations 186 !!------------------------------------------------------------------------------ 187 188 !------------------------------------------------------------------------------! 189 ! 1) Ice thickness distribution parameters initialization 190 !------------------------------------------------------------------------------! 200 191 201 192 !- Types boundaries (integer) … … 266 257 tn_ice(:,:,:) = t_su(:,:,:) 267 258 268 259 END SUBROUTINE lim_itd_ini 269 260 270 261 #else
Note: See TracChangeset
for help on using the changeset viewer.