Changeset 881 for branches/dev_001_SBC/NEMO/LIM_SRC/limthd.F90
- Timestamp:
- 2008-04-08T11:45:52+02:00 (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/dev_001_SBC/NEMO/LIM_SRC/limthd.F90
r879 r881 1 MODULE limthd 1 MODULE limthd_2 2 2 !!====================================================================== 3 !! *** MODULE limthd ***3 !! *** MODULE limthd_2 *** 4 4 !! LIM thermo ice model : ice thermodynamic 5 5 !!====================================================================== … … 8 8 !! 2.0 ! 03-08 (C. Ethe) add lim_thd_init 9 9 !!--------------------------------------------------------------------- 10 #if defined key_ ice_lim11 !!---------------------------------------------------------------------- 12 !! 'key_ ice_lim' : LIMsea-ice model13 !!---------------------------------------------------------------------- 14 !! lim_thd : thermodynamic of sea ice15 !! lim_thd_init : initialisation of sea-ice thermodynamic10 #if defined key_lim2 11 !!---------------------------------------------------------------------- 12 !! 'key_lim2' : LIM 2.0 sea-ice model 13 !!---------------------------------------------------------------------- 14 !! lim_thd_2 : thermodynamic of sea ice 15 !! lim_thd_init_2 : initialisation of sea-ice thermodynamic 16 16 !!---------------------------------------------------------------------- 17 17 USE phycst ! physical constants … … 19 19 USE lbclnk 20 20 USE in_out_manager ! I/O manager 21 USE ice 21 USE ice_2 ! LIM sea-ice variables 22 22 USE ice_oce ! sea-ice/ocean variables 23 23 USE sbc_oce ! 24 24 USE sbc_ice ! 25 USE thd_ice 26 USE dom_ice 27 USE iceini 28 USE limthd_zdf 29 USE limthd_lac 30 USE limtab 25 USE thd_ice_2 ! LIM thermodynamic sea-ice variables 26 USE dom_ice_2 ! LIM sea-ice domain 27 USE iceini_2 28 USE limthd_zdf_2 29 USE limthd_lac_2 30 USE limtab_2 31 31 USE prtctl ! Print control 32 32 … … 34 34 PRIVATE 35 35 36 PUBLIC lim_thd 36 PUBLIC lim_thd_2 ! called by lim_step 37 37 38 38 REAL(wp) :: epsi20 = 1.e-20 , & ! constant values … … 53 53 CONTAINS 54 54 55 SUBROUTINE lim_thd ( kt )55 SUBROUTINE lim_thd_2( kt ) 56 56 !!------------------------------------------------------------------- 57 !! *** ROUTINE lim_thd ***57 !! *** ROUTINE lim_thd_2 *** 58 58 !! 59 59 !! ** Purpose : This routine manages the ice thermodynamic. … … 93 93 !!------------------------------------------------------------------- 94 94 95 IF( kt == nit000 ) CALL lim_thd_init ! Initialization (first time-step only)95 IF( kt == nit000 ) CALL lim_thd_init_2 ! Initialization (first time-step only) 96 96 97 97 !-------------------------------------------! … … 245 245 246 246 ! put the variable in a 1-D array for thermodynamics process 247 CALL tab_2d_1d ( nbpb, frld_1d (1:nbpb) , frld , jpi, jpj, npb(1:nbpb) )248 CALL tab_2d_1d ( nbpb, h_ice_1d (1:nbpb) , hicif , jpi, jpj, npb(1:nbpb) )249 CALL tab_2d_1d ( nbpb, h_snow_1d (1:nbpb) , hsnif , jpi, jpj, npb(1:nbpb) )250 CALL tab_2d_1d ( nbpb, sist_1d (1:nbpb) , sist , jpi, jpj, npb(1:nbpb) )251 CALL tab_2d_1d ( nbpb, tbif_1d (1:nbpb , 1 ), tbif(:,:,1), jpi, jpj, npb(1:nbpb) )252 CALL tab_2d_1d ( nbpb, tbif_1d (1:nbpb , 2 ), tbif(:,:,2), jpi, jpj, npb(1:nbpb) )253 CALL tab_2d_1d ( nbpb, tbif_1d (1:nbpb , 3 ), tbif(:,:,3), jpi, jpj, npb(1:nbpb) )254 CALL tab_2d_1d ( nbpb, qsr_ice_1d (1:nbpb) , qsr_ice , jpi, jpj, npb(1:nbpb) )255 CALL tab_2d_1d ( nbpb, fr1_i0_1d (1:nbpb) , fr1_i0 , jpi, jpj, npb(1:nbpb) )256 CALL tab_2d_1d ( nbpb, fr2_i0_1d (1:nbpb) , fr2_i0 , jpi, jpj, npb(1:nbpb) )257 CALL tab_2d_1d ( nbpb, qns_ice_1d (1:nbpb) , qns_ice , jpi, jpj, npb(1:nbpb) )247 CALL tab_2d_1d_2( nbpb, frld_1d (1:nbpb) , frld , jpi, jpj, npb(1:nbpb) ) 248 CALL tab_2d_1d_2( nbpb, h_ice_1d (1:nbpb) , hicif , jpi, jpj, npb(1:nbpb) ) 249 CALL tab_2d_1d_2( nbpb, h_snow_1d (1:nbpb) , hsnif , jpi, jpj, npb(1:nbpb) ) 250 CALL tab_2d_1d_2( nbpb, sist_1d (1:nbpb) , sist , jpi, jpj, npb(1:nbpb) ) 251 CALL tab_2d_1d_2( nbpb, tbif_1d (1:nbpb , 1 ), tbif(:,:,1), jpi, jpj, npb(1:nbpb) ) 252 CALL tab_2d_1d_2( nbpb, tbif_1d (1:nbpb , 2 ), tbif(:,:,2), jpi, jpj, npb(1:nbpb) ) 253 CALL tab_2d_1d_2( nbpb, tbif_1d (1:nbpb , 3 ), tbif(:,:,3), jpi, jpj, npb(1:nbpb) ) 254 CALL tab_2d_1d_2( nbpb, qsr_ice_1d (1:nbpb) , qsr_ice , jpi, jpj, npb(1:nbpb) ) 255 CALL tab_2d_1d_2( nbpb, fr1_i0_1d (1:nbpb) , fr1_i0 , jpi, jpj, npb(1:nbpb) ) 256 CALL tab_2d_1d_2( nbpb, fr2_i0_1d (1:nbpb) , fr2_i0 , jpi, jpj, npb(1:nbpb) ) 257 CALL tab_2d_1d_2( nbpb, qns_ice_1d (1:nbpb) , qns_ice , jpi, jpj, npb(1:nbpb) ) 258 258 #if ! defined key_coupled 259 CALL tab_2d_1d ( nbpb, qla_ice_1d (1:nbpb) , qla_ice , jpi, jpj, npb(1:nbpb) )260 CALL tab_2d_1d ( nbpb, dqla_ice_1d(1:nbpb) , dqla_ice , jpi, jpj, npb(1:nbpb) )259 CALL tab_2d_1d_2( nbpb, qla_ice_1d (1:nbpb) , qla_ice , jpi, jpj, npb(1:nbpb) ) 260 CALL tab_2d_1d_2( nbpb, dqla_ice_1d(1:nbpb) , dqla_ice , jpi, jpj, npb(1:nbpb) ) 261 261 #endif 262 CALL tab_2d_1d ( nbpb, dqns_ice_1d(1:nbpb) , dqns_ice , jpi, jpj, npb(1:nbpb) )263 CALL tab_2d_1d ( nbpb, tfu_1d (1:nbpb) , tfu , jpi, jpj, npb(1:nbpb) )264 CALL tab_2d_1d ( nbpb, sprecip_1d (1:nbpb) , sprecip , jpi, jpj, npb(1:nbpb) )265 CALL tab_2d_1d ( nbpb, fbif_1d (1:nbpb) , fbif , jpi, jpj, npb(1:nbpb) )266 CALL tab_2d_1d ( nbpb, thcm_1d (1:nbpb) , thcm , jpi, jpj, npb(1:nbpb) )267 CALL tab_2d_1d ( nbpb, qldif_1d (1:nbpb) , qldif , jpi, jpj, npb(1:nbpb) )268 CALL tab_2d_1d ( nbpb, qstbif_1d (1:nbpb) , qstoif , jpi, jpj, npb(1:nbpb) )269 CALL tab_2d_1d ( nbpb, rdmicif_1d (1:nbpb) , rdmicif , jpi, jpj, npb(1:nbpb) )270 CALL tab_2d_1d ( nbpb, dmgwi_1d (1:nbpb) , dmgwi , jpi, jpj, npb(1:nbpb) )271 CALL tab_2d_1d ( nbpb, qlbbq_1d (1:nbpb) , zqlbsbq , jpi, jpj, npb(1:nbpb) )262 CALL tab_2d_1d_2( nbpb, dqns_ice_1d(1:nbpb) , dqns_ice , jpi, jpj, npb(1:nbpb) ) 263 CALL tab_2d_1d_2( nbpb, tfu_1d (1:nbpb) , tfu , jpi, jpj, npb(1:nbpb) ) 264 CALL tab_2d_1d_2( nbpb, sprecip_1d (1:nbpb) , sprecip , jpi, jpj, npb(1:nbpb) ) 265 CALL tab_2d_1d_2( nbpb, fbif_1d (1:nbpb) , fbif , jpi, jpj, npb(1:nbpb) ) 266 CALL tab_2d_1d_2( nbpb, thcm_1d (1:nbpb) , thcm , jpi, jpj, npb(1:nbpb) ) 267 CALL tab_2d_1d_2( nbpb, qldif_1d (1:nbpb) , qldif , jpi, jpj, npb(1:nbpb) ) 268 CALL tab_2d_1d_2( nbpb, qstbif_1d (1:nbpb) , qstoif , jpi, jpj, npb(1:nbpb) ) 269 CALL tab_2d_1d_2( nbpb, rdmicif_1d (1:nbpb) , rdmicif , jpi, jpj, npb(1:nbpb) ) 270 CALL tab_2d_1d_2( nbpb, dmgwi_1d (1:nbpb) , dmgwi , jpi, jpj, npb(1:nbpb) ) 271 CALL tab_2d_1d_2( nbpb, qlbbq_1d (1:nbpb) , zqlbsbq , jpi, jpj, npb(1:nbpb) ) 272 272 273 CALL lim_thd_zdf ( 1, nbpb ) ! compute ice growth273 CALL lim_thd_zdf_2( 1, nbpb ) ! compute ice growth 274 274 275 275 ! back to the geographic grid. 276 CALL tab_1d_2d ( nbpb, frld , npb, frld_1d (1:nbpb) , jpi, jpj )277 CALL tab_1d_2d ( nbpb, hicif , npb, h_ice_1d (1:nbpb) , jpi, jpj )278 CALL tab_1d_2d ( nbpb, hsnif , npb, h_snow_1d (1:nbpb) , jpi, jpj )279 CALL tab_1d_2d ( nbpb, sist , npb, sist_1d (1:nbpb) , jpi, jpj )280 CALL tab_1d_2d ( nbpb, tbif(:,:,1), npb, tbif_1d (1:nbpb , 1 ), jpi, jpj )281 CALL tab_1d_2d ( nbpb, tbif(:,:,2), npb, tbif_1d (1:nbpb , 2 ), jpi, jpj )282 CALL tab_1d_2d ( nbpb, tbif(:,:,3), npb, tbif_1d (1:nbpb , 3 ), jpi, jpj )283 CALL tab_1d_2d ( nbpb, fscmbq , npb, fscbq_1d (1:nbpb) , jpi, jpj )284 CALL tab_1d_2d ( nbpb, ffltbif , npb, fltbif_1d (1:nbpb) , jpi, jpj )285 CALL tab_1d_2d ( nbpb, fstric , npb, fstbif_1d (1:nbpb) , jpi, jpj )286 CALL tab_1d_2d ( nbpb, qldif , npb, qldif_1d (1:nbpb) , jpi, jpj )287 CALL tab_1d_2d ( nbpb, qfvbq , npb, qfvbq_1d (1:nbpb) , jpi, jpj )288 CALL tab_1d_2d ( nbpb, qstoif , npb, qstbif_1d (1:nbpb) , jpi, jpj )289 CALL tab_1d_2d ( nbpb, rdmicif , npb, rdmicif_1d(1:nbpb) , jpi, jpj )290 CALL tab_1d_2d ( nbpb, dmgwi , npb, dmgwi_1d (1:nbpb) , jpi, jpj )291 CALL tab_1d_2d ( nbpb, rdmsnif , npb, rdmsnif_1d(1:nbpb) , jpi, jpj )292 CALL tab_1d_2d ( nbpb, rdvosif , npb, dvsbq_1d (1:nbpb) , jpi, jpj )293 CALL tab_1d_2d ( nbpb, rdvobif , npb, dvbbq_1d (1:nbpb) , jpi, jpj )294 CALL tab_1d_2d ( nbpb, fdvolif , npb, dvlbq_1d (1:nbpb) , jpi, jpj )295 CALL tab_1d_2d ( nbpb, rdvonif , npb, dvnbq_1d (1:nbpb) , jpi, jpj )276 CALL tab_1d_2d_2( nbpb, frld , npb, frld_1d (1:nbpb) , jpi, jpj ) 277 CALL tab_1d_2d_2( nbpb, hicif , npb, h_ice_1d (1:nbpb) , jpi, jpj ) 278 CALL tab_1d_2d_2( nbpb, hsnif , npb, h_snow_1d (1:nbpb) , jpi, jpj ) 279 CALL tab_1d_2d_2( nbpb, sist , npb, sist_1d (1:nbpb) , jpi, jpj ) 280 CALL tab_1d_2d_2( nbpb, tbif(:,:,1), npb, tbif_1d (1:nbpb , 1 ), jpi, jpj ) 281 CALL tab_1d_2d_2( nbpb, tbif(:,:,2), npb, tbif_1d (1:nbpb , 2 ), jpi, jpj ) 282 CALL tab_1d_2d_2( nbpb, tbif(:,:,3), npb, tbif_1d (1:nbpb , 3 ), jpi, jpj ) 283 CALL tab_1d_2d_2( nbpb, fscmbq , npb, fscbq_1d (1:nbpb) , jpi, jpj ) 284 CALL tab_1d_2d_2( nbpb, ffltbif , npb, fltbif_1d (1:nbpb) , jpi, jpj ) 285 CALL tab_1d_2d_2( nbpb, fstric , npb, fstbif_1d (1:nbpb) , jpi, jpj ) 286 CALL tab_1d_2d_2( nbpb, qldif , npb, qldif_1d (1:nbpb) , jpi, jpj ) 287 CALL tab_1d_2d_2( nbpb, qfvbq , npb, qfvbq_1d (1:nbpb) , jpi, jpj ) 288 CALL tab_1d_2d_2( nbpb, qstoif , npb, qstbif_1d (1:nbpb) , jpi, jpj ) 289 CALL tab_1d_2d_2( nbpb, rdmicif , npb, rdmicif_1d(1:nbpb) , jpi, jpj ) 290 CALL tab_1d_2d_2( nbpb, dmgwi , npb, dmgwi_1d (1:nbpb) , jpi, jpj ) 291 CALL tab_1d_2d_2( nbpb, rdmsnif , npb, rdmsnif_1d(1:nbpb) , jpi, jpj ) 292 CALL tab_1d_2d_2( nbpb, rdvosif , npb, dvsbq_1d (1:nbpb) , jpi, jpj ) 293 CALL tab_1d_2d_2( nbpb, rdvobif , npb, dvbbq_1d (1:nbpb) , jpi, jpj ) 294 CALL tab_1d_2d_2( nbpb, fdvolif , npb, dvlbq_1d (1:nbpb) , jpi, jpj ) 295 CALL tab_1d_2d_2( nbpb, rdvonif , npb, dvnbq_1d (1:nbpb) , jpi, jpj ) 296 296 297 297 … … 348 348 349 349 !...Put the variable in a 1-D array for lateral accretion 350 CALL tab_2d_1d ( nbpac, frld_1d (1:nbpac) , frld , jpi, jpj, npac(1:nbpac) )351 CALL tab_2d_1d ( nbpac, h_snow_1d (1:nbpac) , hsnif , jpi, jpj, npac(1:nbpac) )352 CALL tab_2d_1d ( nbpac, h_ice_1d (1:nbpac) , hicif , jpi, jpj, npac(1:nbpac) )353 CALL tab_2d_1d ( nbpac, tbif_1d (1:nbpac , 1 ), tbif(:,:,1), jpi, jpj, npac(1:nbpac) )354 CALL tab_2d_1d ( nbpac, tbif_1d (1:nbpac , 2 ), tbif(:,:,2), jpi, jpj, npac(1:nbpac) )355 CALL tab_2d_1d ( nbpac, tbif_1d (1:nbpac , 3 ), tbif(:,:,3), jpi, jpj, npac(1:nbpac) )356 CALL tab_2d_1d ( nbpac, qldif_1d (1:nbpac) , qldif , jpi, jpj, npac(1:nbpac) )357 CALL tab_2d_1d ( nbpac, qcmif_1d (1:nbpac) , qcmif , jpi, jpj, npac(1:nbpac) )358 CALL tab_2d_1d ( nbpac, qstbif_1d (1:nbpac) , qstoif , jpi, jpj, npac(1:nbpac) )359 CALL tab_2d_1d ( nbpac, rdmicif_1d(1:nbpac) , rdmicif , jpi, jpj, npac(1:nbpac) )360 CALL tab_2d_1d ( nbpac, dvlbq_1d (1:nbpac) , fdvolif , jpi, jpj, npac(1:nbpac) )361 CALL tab_2d_1d ( nbpac, tfu_1d (1:nbpac) , tfu , jpi, jpj, npac(1:nbpac) )350 CALL tab_2d_1d_2( nbpac, frld_1d (1:nbpac) , frld , jpi, jpj, npac(1:nbpac) ) 351 CALL tab_2d_1d_2( nbpac, h_snow_1d (1:nbpac) , hsnif , jpi, jpj, npac(1:nbpac) ) 352 CALL tab_2d_1d_2( nbpac, h_ice_1d (1:nbpac) , hicif , jpi, jpj, npac(1:nbpac) ) 353 CALL tab_2d_1d_2( nbpac, tbif_1d (1:nbpac , 1 ), tbif(:,:,1), jpi, jpj, npac(1:nbpac) ) 354 CALL tab_2d_1d_2( nbpac, tbif_1d (1:nbpac , 2 ), tbif(:,:,2), jpi, jpj, npac(1:nbpac) ) 355 CALL tab_2d_1d_2( nbpac, tbif_1d (1:nbpac , 3 ), tbif(:,:,3), jpi, jpj, npac(1:nbpac) ) 356 CALL tab_2d_1d_2( nbpac, qldif_1d (1:nbpac) , qldif , jpi, jpj, npac(1:nbpac) ) 357 CALL tab_2d_1d_2( nbpac, qcmif_1d (1:nbpac) , qcmif , jpi, jpj, npac(1:nbpac) ) 358 CALL tab_2d_1d_2( nbpac, qstbif_1d (1:nbpac) , qstoif , jpi, jpj, npac(1:nbpac) ) 359 CALL tab_2d_1d_2( nbpac, rdmicif_1d(1:nbpac) , rdmicif , jpi, jpj, npac(1:nbpac) ) 360 CALL tab_2d_1d_2( nbpac, dvlbq_1d (1:nbpac) , fdvolif , jpi, jpj, npac(1:nbpac) ) 361 CALL tab_2d_1d_2( nbpac, tfu_1d (1:nbpac) , tfu , jpi, jpj, npac(1:nbpac) ) 362 362 363 363 ! call lateral accretion routine. 364 CALL lim_thd_lac ( 1 , nbpac )364 CALL lim_thd_lac_2( 1 , nbpac ) 365 365 366 366 ! back to the geographic grid 367 CALL tab_1d_2d ( nbpac, frld , npac(1:nbpac), frld_1d (1:nbpac) , jpi, jpj )368 CALL tab_1d_2d ( nbpac, hsnif , npac(1:nbpac), h_snow_1d (1:nbpac) , jpi, jpj )369 CALL tab_1d_2d ( nbpac, hicif , npac(1:nbpac), h_ice_1d (1:nbpac) , jpi, jpj )370 CALL tab_1d_2d ( nbpac, tbif(:,:,1), npac(1:nbpac), tbif_1d (1:nbpac , 1 ), jpi, jpj )371 CALL tab_1d_2d ( nbpac, tbif(:,:,2), npac(1:nbpac), tbif_1d (1:nbpac , 2 ), jpi, jpj )372 CALL tab_1d_2d ( nbpac, tbif(:,:,3), npac(1:nbpac), tbif_1d (1:nbpac , 3 ), jpi, jpj )373 CALL tab_1d_2d ( nbpac, qstoif , npac(1:nbpac), qstbif_1d (1:nbpac) , jpi, jpj )374 CALL tab_1d_2d ( nbpac, rdmicif , npac(1:nbpac), rdmicif_1d(1:nbpac) , jpi, jpj )375 CALL tab_1d_2d ( nbpac, fdvolif , npac(1:nbpac), dvlbq_1d (1:nbpac) , jpi, jpj )367 CALL tab_1d_2d_2( nbpac, frld , npac(1:nbpac), frld_1d (1:nbpac) , jpi, jpj ) 368 CALL tab_1d_2d_2( nbpac, hsnif , npac(1:nbpac), h_snow_1d (1:nbpac) , jpi, jpj ) 369 CALL tab_1d_2d_2( nbpac, hicif , npac(1:nbpac), h_ice_1d (1:nbpac) , jpi, jpj ) 370 CALL tab_1d_2d_2( nbpac, tbif(:,:,1), npac(1:nbpac), tbif_1d (1:nbpac , 1 ), jpi, jpj ) 371 CALL tab_1d_2d_2( nbpac, tbif(:,:,2), npac(1:nbpac), tbif_1d (1:nbpac , 2 ), jpi, jpj ) 372 CALL tab_1d_2d_2( nbpac, tbif(:,:,3), npac(1:nbpac), tbif_1d (1:nbpac , 3 ), jpi, jpj ) 373 CALL tab_1d_2d_2( nbpac, qstoif , npac(1:nbpac), qstbif_1d (1:nbpac) , jpi, jpj ) 374 CALL tab_1d_2d_2( nbpac, rdmicif , npac(1:nbpac), rdmicif_1d(1:nbpac) , jpi, jpj ) 375 CALL tab_1d_2d_2( nbpac, fdvolif , npac(1:nbpac), dvlbq_1d (1:nbpac) , jpi, jpj ) 376 376 377 377 ENDIF … … 391 391 IF(ln_ctl) THEN 392 392 CALL prt_ctl_info(' lim_thd end ') 393 CALL prt_ctl(tab2d_1=hicif , clinfo1=' hicif : ', tab2d_2=hsnif , clinfo2=' hsnif : ') 394 CALL prt_ctl(tab2d_1=frld , clinfo1=' frld : ', tab2d_2=hicifp , clinfo2=' hicifp : ') 395 CALL prt_ctl(tab2d_1=phicif , clinfo1=' phicif : ', tab2d_2=pfrld , clinfo2=' pfrld : ') 396 CALL prt_ctl(tab2d_1=sist , clinfo1=' sist : ', tab2d_2=tbif(:,:,1), clinfo2=' tbif 1 : ') 397 CALL prt_ctl(tab2d_1=tbif(:,:,2), clinfo1=' tbif 2 : ', tab2d_2=tbif(:,:,3), clinfo2=' tbif 3 : ') 398 CALL prt_ctl(tab2d_1=fdtcn , clinfo1=' fdtcn : ', tab2d_2=qdtcn , clinfo2=' qdtcn : ') 399 CALL prt_ctl(tab2d_1=qstoif , clinfo1=' qstoif : ', tab2d_2=fsbbq , clinfo2=' fsbbq : ') 393 CALL prt_ctl(tab2d_1=hicif , clinfo1=' lim_thd: hicif : ', tab2d_2=hsnif , clinfo2=' hsnif : ') 394 CALL prt_ctl(tab2d_1=frld , clinfo1=' lim_thd: frld : ', tab2d_2=hicifp, clinfo2=' hicifp : ') 395 CALL prt_ctl(tab2d_1=phicif, clinfo1=' lim_thd: phicif : ', tab2d_2=pfrld , clinfo2=' pfrld : ') 396 CALL prt_ctl(tab2d_1=sist , clinfo1=' lim_thd: sist : ') 397 CALL prt_ctl(tab2d_1=tbif(:,:,1), clinfo1=' lim_thd: tbif 1 : ') 398 CALL prt_ctl(tab2d_1=tbif(:,:,2), clinfo1=' lim_thd: tbif 2 : ') 399 CALL prt_ctl(tab2d_1=tbif(:,:,3), clinfo1=' lim_thd: tbif 3 : ') 400 CALL prt_ctl(tab2d_1=fdtcn , clinfo1=' lim_thd: fdtcn : ', tab2d_2=qdtcn , clinfo2=' qdtcn : ') 401 CALL prt_ctl(tab2d_1=qstoif, clinfo1=' lim_thd: qstoif : ', tab2d_2=fsbbq , clinfo2=' fsbbq : ') 400 402 ENDIF 401 403 ! 402 END SUBROUTINE lim_thd403 404 405 SUBROUTINE lim_thd_init404 END SUBROUTINE lim_thd_2 405 406 407 SUBROUTINE lim_thd_init_2 406 408 !!------------------------------------------------------------------- 407 !! *** ROUTINE lim_thd_init ***409 !! *** ROUTINE lim_thd_init_2 *** 408 410 !! 409 411 !! ** Purpose : Physical constants and parameters linked to the ice … … 427 429 IF(lwp) THEN 428 430 WRITE(numout,*) 429 WRITE(numout,*)'lim_thd_init 430 WRITE(numout,*)'~~~~~~~~~~~~ '431 WRITE(numout,*)'lim_thd_init_2: ice parameters for ice thermodynamic computation ' 432 WRITE(numout,*)'~~~~~~~~~~~~~~' 431 433 WRITE(numout,*)' maximum melting at the bottom hmelt = ', hmelt 432 434 WRITE(numout,*)' ice thick. for lateral accretion in NH (SH) hiccrit(1/2) = ', hiccrit … … 458 460 ENDIF 459 461 460 END SUBROUTINE lim_thd_init 462 END SUBROUTINE lim_thd_init_2 461 463 462 464 #else 463 465 !!---------------------------------------------------------------------- 464 !! Default option Dummy module NO LIMsea-ice model466 !! Default option Dummy module NO LIM 2.0 sea-ice model 465 467 !!---------------------------------------------------------------------- 466 468 CONTAINS 467 SUBROUTINE lim_thd ! Dummy routine468 END SUBROUTINE lim_thd 469 SUBROUTINE lim_thd_2 ! Dummy routine 470 END SUBROUTINE lim_thd_2 469 471 #endif 470 472 471 473 !!====================================================================== 472 END MODULE limthd 474 END MODULE limthd_2
Note: See TracChangeset
for help on using the changeset viewer.