- Timestamp:
- 2015-12-08T12:39:53+01:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/icebergs_restart_single_file/NEMOGCM/NEMO/OPA_SRC/CRS/crsini.F90
r6019 r6020 7 7 !!---------------------------------------------------------------------- 8 8 9 USE timing ! Timing 9 !!---------------------------------------------------------------------- 10 !! crs_init : 11 !!---------------------------------------------------------------------- 12 USE par_kind, ONLY: wp 10 13 USE par_oce ! For parameter jpi,jpj,jphgr_msh 11 14 USE dom_oce ! For parameters in par_oce (jperio, lk_vvl) 12 USE crs ! Coarse grid domain15 USE crs ! Coarse grid domain 13 16 USE phycst, ONLY: omega, rad ! physical constants 14 USE wrk_nemo15 USE in_out_manager16 USE par_kind, ONLY: wp17 USE iom18 17 USE crsdom 19 18 USE crsdomwri 20 19 USE crslbclnk 20 ! 21 USE iom 22 USE in_out_manager 21 23 USE lib_mpp 24 USE wrk_nemo 25 USE timing ! Timing 22 26 23 27 IMPLICIT NONE 24 28 PRIVATE 25 29 26 PUBLIC crs_init30 PUBLIC crs_init ! called by nemogcm.F90 module 27 31 28 32 !! * Substitutions 29 33 # include "domzgr_substitute.h90" 30 34 !!---------------------------------------------------------------------- 31 35 !! $Id$ 36 !!---------------------------------------------------------------------- 32 37 CONTAINS 33 38 … … 65 70 !! - Read in pertinent data ? 66 71 !!------------------------------------------------------------------- 67 !! Local variables68 72 INTEGER :: ji,jj,jk ! dummy indices 69 73 INTEGER :: ierr ! allocation error status … … 183 187 184 188 ! 185 CALL wrk_alloc( jpi, jpj, jpk,zfse3t, zfse3u, zfse3v, zfse3w )189 CALL wrk_alloc( jpi,jpj,jpk, zfse3t, zfse3u, zfse3v, zfse3w ) 186 190 ! 187 191 zfse3t(:,:,:) = fse3t(:,:,:) … … 200 204 ! 3.d.3 Vertical scale factors 201 205 ! 202 203 204 206 CALL crs_dom_e3( e1t, e2t, zfse3t, e1e2w_crs, 'T', tmask, e3t_crs, e3t_max_crs) 205 207 CALL crs_dom_e3( e1u, e2u, zfse3u, e2e3u_crs, 'U', umask, e3u_crs, e3u_max_crs) … … 207 209 CALL crs_dom_e3( e1t, e2t, zfse3w, e1e2w_crs, 'W', tmask, e3w_crs, e3w_max_crs) 208 210 209 ! Re set 0 toe3t_0 or e3w_0211 ! Replace 0 by e3t_0 or e3w_0 210 212 DO jk = 1, jpk 211 213 DO ji = 1, jpi_crs … … 247 249 ENDIF 248 250 249 !--------------------------------------------------------- 250 ! 7. Finish and clean-up 251 !--------------------------------------------------------- 252 CALL wrk_dealloc(jpi, jpj, jpk, zfse3t, zfse3u, zfse3v, zfse3w ) 253 254 251 !--------------------------------------------------------- 252 ! 7. Finish and clean-up 253 !--------------------------------------------------------- 254 CALL wrk_dealloc( jpi,jpj,jpk, zfse3t, zfse3u, zfse3v, zfse3w ) 255 ! 255 256 END SUBROUTINE crs_init 256 257 257 258 !!====================================================================== 258 259 259 END MODULE crsini
Note: See TracChangeset
for help on using the changeset viewer.