- Timestamp:
- 2013-11-20T19:24:53+01:00 (10 years ago)
- Location:
- branches/2013/dev_MERGE_2013/NEMOGCM/NEMO/LIM_SRC_3
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2013/dev_MERGE_2013/NEMOGCM/NEMO/LIM_SRC_3/iceini.F90
r4293 r4298 71 71 & 'use more ocean levels or less ice/snow layers/categories.' ) 72 72 73 ! ! Open the namelist file 74 CALL ctl_opn( numnam_ice, 'namelist_ice', 'OLD', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp ) 73 ! Open the reference and configuration namelist files and namelist output file 74 CALL ctl_opn( numnam_ice_ref, 'namelist_ice_ref', 'OLD', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp ) 75 CALL ctl_opn( numnam_ice_cfg, 'namelist_ice_cfg', 'OLD', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp ) 76 CALL ctl_opn( numoni, 'output.namelist.ice', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp ) 75 77 ! 76 78 CALL ice_run ! set some ice run parameters … … 129 131 !!------------------------------------------------------------------- 130 132 NAMELIST/namicerun/ cn_icerst_in, cn_icerst_out, ln_limdyn, amax, cai, cao, ln_nicep, ln_limdiahsb, ln_limdiaout 133 INTEGER :: ios ! Local integer output status for namelist read 131 134 !!------------------------------------------------------------------- 132 135 ! 133 REWIND( numnam_ice ) ! Read Namelist namicerun 134 READ ( numnam_ice , namicerun ) 135 ! 136 !IF( lk_mpp .AND. ln_nicep ) THEN 137 ! ln_nicep = .FALSE. 138 ! CALL ctl_warn( 'ice_run : specific control print for LIM3 desactivated with MPI' ) 139 !ENDIF 136 REWIND( numnam_ice_ref ) ! Namelist namicerun in reference namelist : Parameters for ice 137 READ ( numnam_ice_ref, namicerun, IOSTAT = ios, ERR = 901) 138 901 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namicerun in reference namelist', lwp ) 139 140 REWIND( numnam_ice_cfg ) ! Namelist namicerun in configuration namelist : Parameters for ice 141 READ ( numnam_ice_cfg, namicerun, IOSTAT = ios, ERR = 902 ) 142 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namicerun in configuration namelist', lwp ) 143 WRITE ( numoni, namicerun ) 144 ! 145 IF( lk_mpp .AND. ln_nicep ) THEN 146 ln_nicep = .FALSE. 147 CALL ctl_warn( 'ice_run : specific control print for LIM3 desactivated with MPI' ) 148 ENDIF 140 149 ! 141 150 IF(lwp) THEN ! control print -
branches/2013/dev_MERGE_2013/NEMOGCM/NEMO/LIM_SRC_3/limcat_1D.F90
r4269 r4298 53 53 !-------------------------------------------------------------------- 54 54 ijpij = SIZE(zhti,1) 55 zht_i(1:ijpij,1:jpl) = 0. d056 zht_s(1:ijpij,1:jpl) = 0. d057 za_i (1:ijpij,1:jpl) = 0. d055 zht_i(1:ijpij,1:jpl) = 0._wp 56 zht_s(1:ijpij,1:jpl) = 0._wp 57 za_i (1:ijpij,1:jpl) = 0._wp 58 58 59 59 !------------------------------------------------------------------------------------ … … 101 101 102 102 ! initialisation of ice variables for each try 103 zht_i(ji,1:jpl) = 0. d0104 za_i (ji,1:jpl) = 0. d0103 zht_i(ji,1:jpl) = 0._wp 104 za_i (ji,1:jpl) = 0._wp 105 105 106 106 ! *** case very thin ice: fill only category 1 … … 114 114 ! Fill ice thicknesses except the last one (i_fill) by (hmax-hmin)/2 115 115 DO jl = 1, i_fill - 1 116 zht_i(ji,jl) = ( zhi_max(jl) + zhi_max(jl-1) ) / 2.116 zht_i(ji,jl) = ( zhi_max(jl) + zhi_max(jl-1) ) * 0.5_wp 117 117 END DO 118 118 … … 130 130 DO jl = 1, i_fill - 1 131 131 IF ( jl == jl0 ) CYCLE 132 zarg = ( zht_i(ji,jl) - zhti(ji) ) / ( zhti(ji) / 2.)132 zarg = ( zht_i(ji,jl) - zhti(ji) ) / ( zhti(ji) * 0.5_wp ) 133 133 za_i(ji,jl) = za_i (ji,jl0) * EXP(-zarg**2) 134 134 END DO … … 160 160 ztest_4 = 1 161 161 DO jl = 1, i_fill 162 IF ( za_i(ji,jl) < 0. 0d0) ztest_4 = 0162 IF ( za_i(ji,jl) < 0._wp ) ztest_4 = 0 163 163 END DO 164 164 -
branches/2013/dev_MERGE_2013/NEMOGCM/NEMO/LIM_SRC_3/limitd_me.F90
r4293 r4298 1477 1477 !! ** input : Namelist namiceitdme 1478 1478 !!------------------------------------------------------------------- 1479 INTEGER :: ios ! Local integer output status for namelist read 1479 1480 NAMELIST/namiceitdme/ ridge_scheme_swi, Cs, Cf, fsnowrdg, fsnowrft,& 1480 1481 Gstar, astar, & … … 1484 1485 !!------------------------------------------------------------------- 1485 1486 ! 1486 REWIND( numnam_ice ) ! read namiceitdme namelist 1487 READ ( numnam_ice , namiceitdme) 1487 REWIND( numnam_ice_ref ) ! Namelist namicetdme in reference namelist : Ice mechanical ice redistribution 1488 READ ( numnam_ice_ref, namiceitdme, IOSTAT = ios, ERR = 901) 1489 901 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namiceitdme in reference namelist', lwp ) 1490 1491 REWIND( numnam_ice_cfg ) ! Namelist namiceitdme in configuration namelist : Ice mechanical ice redistribution 1492 READ ( numnam_ice_cfg, namiceitdme, IOSTAT = ios, ERR = 902 ) 1493 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namiceitdme in configuration namelist', lwp ) 1494 WRITE ( numoni, namiceitdme ) 1488 1495 ! 1489 1496 IF (lwp) THEN ! control print
Note: See TracChangeset
for help on using the changeset viewer.