- Timestamp:
- 2017-09-06T19:02:43+02:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/icedia.F90
r8486 r8505 31 31 PRIVATE 32 32 33 PUBLIC ice_dia ! routinecalled by icestp.F9034 PUBLIC ice_dia_init ! routinecalled in icestp.F9033 PUBLIC ice_dia ! called by icestp.F90 34 PUBLIC ice_dia_init ! called in icestp.F90 35 35 36 36 REAL(wp), DIMENSION(:,:), ALLOCATABLE :: vol_loc_ini, sal_loc_ini, tem_loc_ini ! initial volume, salt and heat contents … … 162 162 !! - Compute coefficients for conversion 163 163 !!--------------------------------------------------------------------------- 164 INTEGER :: i error ! local integer165 !! 166 !!NAMELIST/namicehsb/ blabla164 INTEGER :: ios, ierror ! local integer 165 !! 166 NAMELIST/namicediag/ ln_limdiachk, ln_limdiahsb, ln_limctl, iiceprt, jiceprt 167 167 !!---------------------------------------------------------------------- 168 168 ! 169 !!REWIND ( numnam_ice ) ! Read Namelist namicehsb 170 !!READ ( numnam_ice, namicehsb ) 171 ! 172 IF(lwp) THEN ! Control print 169 REWIND( numnam_ice_ref ) ! Namelist namicediag in reference namelist : Parameters for ice 170 READ ( numnam_ice_ref, namicediag, IOSTAT = ios, ERR = 901) 171 901 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namicediag in reference namelist', lwp ) 172 173 REWIND( numnam_ice_cfg ) ! Namelist namicediag in configuration namelist : Parameters for ice 174 READ ( numnam_ice_cfg, namicediag, IOSTAT = ios, ERR = 902 ) 175 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namicediag in configuration namelist', lwp ) 176 IF(lwm) WRITE ( numoni, namicediag ) 177 ! 178 IF(lwp) THEN ! control print 173 179 WRITE(numout,*) 174 WRITE(numout,*) 'ice_dia_init : check the heat and salt budgets' 175 WRITE(numout,*) '~~~~~~~~~~~~' 180 WRITE(numout,*) 'ice_dia_init: ice diagnostics' 181 WRITE(numout,*) ' ~~~~~~~~~~~' 182 WRITE(numout,*) ' Namelist namicediag : ' 183 WRITE(numout,*) ' Diagnose online heat/mass/salt budget ln_limdiachk = ', ln_limdiachk 184 WRITE(numout,*) ' Output heat/mass/salt budget ln_limdiahsb = ', ln_limdiahsb 185 WRITE(numout,*) ' control prints for a given grid point ln_limctl = ', ln_limctl 186 WRITE(numout,*) ' chosen grid point position (iiceprt,jiceprt) = (', iiceprt,',', jiceprt,')' 176 187 ENDIF 177 188 ! 178 ALLOCATE( vol_loc_ini(jpi,jpj), sal_loc_ini(jpi,jpj), tem_loc_ini(jpi,jpj), STAT=ierror ) 179 IF( ierror > 0 ) THEN 180 CALL ctl_stop( 'ice_dia: unable to allocate vol_loc_ini' ) 181 RETURN 182 ENDIF 183 ! 184 CALL ice_dia_rst( 'READ' ) !* read or initialize all required files 189 IF( ln_limdiahsb ) THEN 190 ALLOCATE( vol_loc_ini(jpi,jpj), sal_loc_ini(jpi,jpj), tem_loc_ini(jpi,jpj), STAT=ierror ) 191 IF( ierror > 0 ) THEN 192 CALL ctl_stop( 'ice_dia: unable to allocate vol_loc_ini' ) 193 RETURN 194 ENDIF 195 ! 196 CALL ice_dia_rst( 'READ' ) !* read or initialize all required files 197 ENDIF 185 198 ! 186 199 END SUBROUTINE ice_dia_init
Note: See TracChangeset
for help on using the changeset viewer.