Changeset 5055
- Timestamp:
- 2015-02-04T11:10:53+01:00 (10 years ago)
- Location:
- branches/2015/dev_r5044_CNRS_LIM3CLEAN/NEMOGCM/NEMO
- Files:
-
- 1 deleted
- 30 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/dev_r5044_CNRS_LIM3CLEAN/NEMOGCM/NEMO/LIM_SRC_3/dom_ice.F90
r4161 r5055 7 7 !! 4.0 ! 2011-02 (G. Madec) dynamical allocation 8 8 !!---------------------------------------------------------------------- 9 USE par_ice ! LIM-3 parameter10 9 USE in_out_manager ! I/O manager 11 10 USE lib_mpp ! MPP library -
branches/2015/dev_r5044_CNRS_LIM3CLEAN/NEMOGCM/NEMO/LIM_SRC_3/ice.F90
r5053 r5055 11 11 !! 'key_lim3' LIM-3 sea-ice model 12 12 !!---------------------------------------------------------------------- 13 USE par_ice ! LIM sea-ice parameters14 13 USE in_out_manager ! I/O manager 15 14 USE lib_mpp ! MPP library … … 366 365 !!-------------------------------------------------------------------------- 367 366 ! !!: ** Namelist namicerun read in sbc_lim_init ** 368 !clemINTEGER , PUBLIC :: jpl !: number of ice categories369 !clemINTEGER , PUBLIC :: nlay_i !: number of ice layers370 !clemINTEGER , PUBLIC :: nlay_s !: number of snow layers367 INTEGER , PUBLIC :: jpl !: number of ice categories 368 INTEGER , PUBLIC :: nlay_i !: number of ice layers 369 INTEGER , PUBLIC :: nlay_s !: number of snow layers 371 370 CHARACTER(len=32) , PUBLIC :: cn_icerst_in !: suffix of ice restart name (input) 372 371 CHARACTER(len=32) , PUBLIC :: cn_icerst_out !: suffix of ice restart name (output) -
branches/2015/dev_r5044_CNRS_LIM3CLEAN/NEMOGCM/NEMO/LIM_SRC_3/limcons.F90
r4873 r5055 16 16 !!---------------------------------------------------------------------- 17 17 USE phycst ! physical constants 18 USE par_ice ! LIM-3 parameter19 18 USE ice ! LIM-3 variables 20 19 USE dom_ice ! LIM-3 domain -
branches/2015/dev_r5044_CNRS_LIM3CLEAN/NEMOGCM/NEMO/LIM_SRC_3/limctl.F90
r5051 r5055 15 15 USE oce ! ocean dynamics and tracers 16 16 USE dom_oce ! ocean space and time domain 17 USE par_ice ! sea-ice parameters18 17 USE ice ! LIM-3: ice variables 19 18 USE thd_ice ! LIM-3: thermodynamical variables -
branches/2015/dev_r5044_CNRS_LIM3CLEAN/NEMOGCM/NEMO/LIM_SRC_3/limdiahsb.F90
r5051 r5055 14 14 !!---------------------------------------------------------------------- 15 15 USE ice ! LIM-3: sea-ice variable 16 USE par_ice ! LIM-3: ice parameters17 16 USE dom_ice ! LIM-3: sea-ice domain 18 17 USE dom_oce ! ocean domain -
branches/2015/dev_r5044_CNRS_LIM3CLEAN/NEMOGCM/NEMO/LIM_SRC_3/limdyn.F90
r5053 r5055 20 20 USE sbc_ice ! Surface boundary condition: ice fields 21 21 USE ice ! LIM-3 variables 22 USE par_ice ! LIM-3 parameters23 22 USE dom_ice ! LIM-3 domain 24 23 USE limrhg ! LIM-3 rheology -
branches/2015/dev_r5044_CNRS_LIM3CLEAN/NEMOGCM/NEMO/LIM_SRC_3/limistate.F90
r5051 r5055 22 22 USE eosbn2 ! equation of state 23 23 USE ice ! sea-ice variables 24 USE par_ice ! ice parameters25 24 USE par_oce ! ocean parameters 26 25 USE dom_ice ! sea-ice domain -
branches/2015/dev_r5044_CNRS_LIM3CLEAN/NEMOGCM/NEMO/LIM_SRC_3/limitd_me.F90
r5053 r5055 18 18 USE thd_ice ! LIM thermodynamics 19 19 USE ice ! LIM variables 20 USE par_ice ! LIM parameters21 20 USE dom_ice ! LIM domain 22 21 USE limthd_lac ! LIM -
branches/2015/dev_r5044_CNRS_LIM3CLEAN/NEMOGCM/NEMO/LIM_SRC_3/limitd_th.F90
r5053 r5055 24 24 USE thd_ice ! LIM-3 thermodynamic variables 25 25 USE ice ! LIM-3 variables 26 USE par_ice ! LIM-3 parameters27 26 USE limvar ! LIM-3 variables 28 27 USE limcons ! LIM-3 conservation -
branches/2015/dev_r5044_CNRS_LIM3CLEAN/NEMOGCM/NEMO/LIM_SRC_3/limrst.F90
r5053 r5055 18 18 USE ice ! sea-ice variables 19 19 USE oce , ONLY : snwice_mass, snwice_mass_b 20 USE par_ice ! sea-ice parameters21 20 USE dom_oce ! ocean domain 22 21 USE sbc_oce ! Surface boundary condition: ocean fields -
branches/2015/dev_r5044_CNRS_LIM3CLEAN/NEMOGCM/NEMO/LIM_SRC_3/limsbc.F90
r5053 r5055 25 25 USE par_oce ! ocean parameters 26 26 USE phycst ! physical constants 27 USE par_ice ! ice parameters28 27 USE dom_oce ! ocean domain 29 28 USE dom_ice, ONLY : tms, area -
branches/2015/dev_r5044_CNRS_LIM3CLEAN/NEMOGCM/NEMO/LIM_SRC_3/limthd.F90
r5053 r5055 24 24 USE oce , ONLY : fraqsr_1lev 25 25 USE ice ! LIM: sea-ice variables 26 USE par_ice ! LIM: sea-ice parameters27 26 USE sbc_oce ! Surface boundary condition: ocean fields 28 27 USE sbc_ice ! Surface boundary condition: ice fields -
branches/2015/dev_r5044_CNRS_LIM3CLEAN/NEMOGCM/NEMO/LIM_SRC_3/limthd_dh.F90
r5053 r5055 20 20 USE sbc_oce ! Surface boundary condition: ocean fields 21 21 USE ice ! LIM variables 22 USE par_ice ! LIM parameters23 22 USE thd_ice ! LIM thermodynamics 24 23 USE in_out_manager ! I/O manager -
branches/2015/dev_r5044_CNRS_LIM3CLEAN/NEMOGCM/NEMO/LIM_SRC_3/limthd_dif.F90
r5051 r5055 19 19 USE phycst ! physical constants (ocean directory) 20 20 USE ice ! LIM-3 variables 21 USE par_ice ! LIM-3 parameters22 21 USE thd_ice ! LIM-3: thermodynamics 23 22 USE in_out_manager ! I/O manager -
branches/2015/dev_r5044_CNRS_LIM3CLEAN/NEMOGCM/NEMO/LIM_SRC_3/limthd_ent.F90
r5051 r5055 25 25 USE sbc_oce ! Surface boundary condition: ocean fields 26 26 USE ice ! LIM variables 27 USE par_ice ! LIM parameters28 27 USE thd_ice ! LIM thermodynamics 29 28 USE limvar ! LIM variables -
branches/2015/dev_r5044_CNRS_LIM3CLEAN/NEMOGCM/NEMO/LIM_SRC_3/limthd_lac.F90
r5047 r5055 22 22 USE thd_ice ! LIM thermodynamics 23 23 USE dom_ice ! LIM domain 24 USE par_ice ! LIM parameters25 24 USE ice ! LIM variables 26 25 USE limtab ! LIM 2D <==> 1D -
branches/2015/dev_r5044_CNRS_LIM3CLEAN/NEMOGCM/NEMO/LIM_SRC_3/limthd_sal.F90
r5051 r5055 18 18 USE sbc_oce ! Surface boundary condition: ocean fields 19 19 USE ice ! LIM variables 20 USE par_ice ! LIM parameters21 20 USE thd_ice ! LIM thermodynamics 22 21 USE limvar ! LIM variables -
branches/2015/dev_r5044_CNRS_LIM3CLEAN/NEMOGCM/NEMO/LIM_SRC_3/limtrp.F90
r5053 r5055 17 17 USE dom_oce ! ocean domain 18 18 USE sbc_oce ! ocean surface boundary condition 19 USE par_ice ! ice parameter20 19 USE dom_ice ! ice domain 21 20 USE ice ! ice variables -
branches/2015/dev_r5044_CNRS_LIM3CLEAN/NEMOGCM/NEMO/LIM_SRC_3/limupdate1.F90
r5053 r5055 31 31 USE limrst 32 32 USE thd_ice ! LIM thermodynamic sea-ice variables 33 USE par_ice34 33 USE limitd_th 35 34 USE limitd_me -
branches/2015/dev_r5044_CNRS_LIM3CLEAN/NEMOGCM/NEMO/LIM_SRC_3/limupdate2.F90
r5053 r5055 31 31 USE limrst 32 32 USE thd_ice ! LIM thermodynamic sea-ice variables 33 USE par_ice34 33 USE limitd_th 35 34 USE limitd_me -
branches/2015/dev_r5044_CNRS_LIM3CLEAN/NEMOGCM/NEMO/LIM_SRC_3/limvar.F90
r5053 r5055 40 40 USE sbc_oce ! Surface boundary condition: ocean fields 41 41 USE ice ! ice variables 42 USE par_ice ! ice parameters43 42 USE thd_ice ! ice variables (thermodynamics) 44 43 USE dom_ice ! ice domain -
branches/2015/dev_r5044_CNRS_LIM3CLEAN/NEMOGCM/NEMO/LIM_SRC_3/limwri.F90
r5051 r5055 24 24 USE lib_mpp ! MPP library 25 25 USE wrk_nemo ! work arrays 26 USE par_ice27 26 USE iom 28 27 USE timing ! Timing -
branches/2015/dev_r5044_CNRS_LIM3CLEAN/NEMOGCM/NEMO/LIM_SRC_3/thd_ice.F90
r5053 r5055 6 6 !! History : 3.0 ! 2002-11 (C. Ethe) F90: Free form and module 7 7 !!---------------------------------------------------------------------- 8 USE par_ice ! LIM-3 parameters9 8 USE in_out_manager ! I/O manager 10 9 USE lib_mpp ! MPP library 10 USE ice, ONLY : nlay_i, nlay_s 11 11 12 12 IMPLICIT NONE -
branches/2015/dev_r5044_CNRS_LIM3CLEAN/NEMOGCM/NEMO/OPA_SRC/BDY/bdydta.F90
r5053 r5055 33 33 USE ice_2 34 34 #elif defined key_lim3 35 USE par_ice36 35 USE ice 37 36 USE limvar ! redistribute ice input into categories -
branches/2015/dev_r5044_CNRS_LIM3CLEAN/NEMOGCM/NEMO/OPA_SRC/BDY/bdyice_lim.F90
r5053 r5055 26 26 USE dom_ice_2 ! sea-ice domain 27 27 #elif defined key_lim3 28 USE par_ice29 28 USE ice ! LIM_3 ice variables 30 29 USE dom_ice ! sea-ice domain -
branches/2015/dev_r5044_CNRS_LIM3CLEAN/NEMOGCM/NEMO/OPA_SRC/DIA/diadct.F90
r4624 r5055 42 42 #endif 43 43 #if defined key_lim3 44 USE par_ice45 44 USE ice 46 45 #endif -
branches/2015/dev_r5044_CNRS_LIM3CLEAN/NEMOGCM/NEMO/OPA_SRC/IOM/iom.F90
r5043 r5055 33 33 USE icb_oce, ONLY : nclasses, class_num ! !: iceberg classes 34 34 #if defined key_lim3 35 USE par_ice35 USE ice , ONLY : jpl 36 36 #elif defined key_lim2 37 37 USE par_ice_2 -
branches/2015/dev_r5044_CNRS_LIM3CLEAN/NEMOGCM/NEMO/OPA_SRC/SBC/sbc_ice.F90
r5009 r5055 16 16 USE sbc_oce ! surface boundary condition: ocean 17 17 # if defined key_lim3 18 USE par_ice! LIM-3 parameters18 USE ice ! LIM-3 parameters 19 19 # endif 20 20 # if defined key_lim2 -
branches/2015/dev_r5044_CNRS_LIM3CLEAN/NEMOGCM/NEMO/OPA_SRC/SBC/sbccpl.F90
r5009 r5055 24 24 USE phycst ! physical constants 25 25 #if defined key_lim3 26 USE par_ice ! ice parameters27 26 USE ice ! ice variables 28 27 #endif -
branches/2015/dev_r5044_CNRS_LIM3CLEAN/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_lim.F90
r5053 r5055 22 22 USE oce ! ocean dynamics and tracers 23 23 USE dom_oce ! ocean space and time domain 24 USE par_ice ! sea-ice parameters25 24 USE ice ! LIM-3: ice variables 26 25 USE thd_ice ! LIM-3: thermodynamical variables … … 304 303 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~ via Louvain la Neuve Ice Model (LIM-3) time stepping' 305 304 ! 305 ! Open the reference and configuration namelist files and namelist output file 306 CALL ctl_opn( numnam_ice_ref, 'namelist_ice_ref', 'OLD', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp ) 307 CALL ctl_opn( numnam_ice_cfg, 'namelist_ice_cfg', 'OLD', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp ) 308 IF(lwm) CALL ctl_opn( numoni, 'output.namelist.ice', 'UNKNOWN', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp, 1 ) 309 310 CALL ice_run ! set some ice run parameters including jpl, nlay_i and nlay_s 306 311 ! ! Allocate the ice arrays 307 312 ierr = ice_alloc () ! ice variables … … 316 321 ! ! adequation jpk versus ice/snow layers/categories 317 322 IF( jpl > jpk .OR. (nlay_i+1) > jpk .OR. nlay_s > jpk ) & 318 & CALL ctl_stop( 'STOP', &323 & CALL ctl_stop( 'STOP', & 319 324 & 'sbc_lim_init: the 3rd dimension of workspace arrays is too small.', & 320 325 & 'use more ocean levels or less ice/snow layers/categories.' ) 321 322 ! Open the reference and configuration namelist files and namelist output file 323 CALL ctl_opn( numnam_ice_ref, 'namelist_ice_ref', 'OLD', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp ) 324 CALL ctl_opn( numnam_ice_cfg, 'namelist_ice_cfg', 'OLD', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp ) 325 IF(lwm) CALL ctl_opn( numoni, 'output.namelist.ice', 'UNKNOWN', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp, 1 ) 326 ! 327 CALL ice_run ! set some ice run parameters 326 ! 328 327 ! 329 328 CALL lim_thd_init ! set ice thermodynics parameters … … 335 334 CALL lim_itd_init ! ice thickness distribution initialization 336 335 ! 337 CALL lim_itd_me_init ! ice thickness distribution initialization 336 CALL lim_itd_me_init ! ice thickness distribution initialization for mecanical deformation 338 337 ! ! Initial sea-ice state 339 338 IF( .NOT. ln_rstart ) THEN ! start from rest: sea-ice deduced from sst … … 373 372 !! ** input : Namelist namicerun 374 373 !!------------------------------------------------------------------- 375 !clem NAMELIST/namicerun/ jpl, nlay_i, nlay_s, cn_icerst_in, cn_icerst_out, & 376 !clem & ln_limdyn, amax, ln_nicep, ln_limdiahsb, ln_limdiaout 377 NAMELIST/namicerun/ cn_icerst_in, cn_icerst_out, & 374 INTEGER :: ios ! Local integer output status for namelist read 375 NAMELIST/namicerun/ jpl, nlay_i, nlay_s, cn_icerst_in, cn_icerst_out, & 378 376 & ln_limdyn, amax, ln_nicep, ln_limdiahsb, ln_limdiaout 379 INTEGER :: ios ! Local integer output status for namelist read380 377 !!------------------------------------------------------------------- 381 378 ! … … 394 391 WRITE(numout,*) 'ice_run : ice share parameters for dynamics/advection/thermo of sea-ice' 395 392 WRITE(numout,*) ' ~~~~~~' 396 !clemWRITE(numout,*) ' number of ice categories = ', jpl397 !clemWRITE(numout,*) ' number of ice layers = ', nlay_i398 !clemWRITE(numout,*) ' number of snow layers = ', nlay_s393 WRITE(numout,*) ' number of ice categories = ', jpl 394 WRITE(numout,*) ' number of ice layers = ', nlay_i 395 WRITE(numout,*) ' number of snow layers = ', nlay_s 399 396 WRITE(numout,*) ' switch for ice dynamics (1) or not (0) ln_limdyn = ', ln_limdyn 400 397 WRITE(numout,*) ' maximum ice concentration = ', amax
Note: See TracChangeset
for help on using the changeset viewer.