Parent Directory
| 
Revision Log
| 
Revision Graph
Changed Z co-ordinate and lopping code: o Z now increases upwards (ie. for the ocean, z<0) This needed changes only in ini_vertical_grid() where zC and zFace are calculated and in ini_masks_etc() where hFacC is calculated using IF's (eg. IF H(i,j,bi,bj).LT.zFace(K) ...). o The lopping code needed a correction to the hydrostatic pressure calculation. o Two new parameters control the lopping: hFacMin - is the lower limit for water fraction in a cell hFacMinDz - is the lower limit for water thickness in lopped cells Both default to 0 so that pure lopping is the default. To force full cells, using continuous bathymetry, set hFacMin=1
| 1 | C $Header: /u/gcmpack/models/MITgcmUV/model/src/ini_vertical_grid.F,v 1.1 1998/07/02 14:17:11 adcroft Exp $ |
| 2 | |
| 3 | #include "CPP_EEOPTIONS.h" |
| 4 | |
| 5 | CStartOfInterface |
| 6 | SUBROUTINE INI_VERTICAL_GRID( myThid ) |
| 7 | C /==========================================================\ |
| 8 | C | SUBROUTINE INI_VERTICAL_GRID | |
| 9 | C | o Initialise vertical gridding arrays | |
| 10 | C |==========================================================| |
| 11 | C | | |
| 12 | C \==========================================================/ |
| 13 | |
| 14 | C === Global variables === |
| 15 | #include "SIZE.h" |
| 16 | #include "EEPARAMS.h" |
| 17 | #include "PARAMS.h" |
| 18 | #include "GRID.h" |
| 19 | |
| 20 | C == Routine arguments == |
| 21 | C myThid - Number of this instance of INI_DEPTHS |
| 22 | INTEGER myThid |
| 23 | CEndOfInterface |
| 24 | |
| 25 | C == Local variables == |
| 26 | C K |
| 27 | INTEGER K |
| 28 | |
| 29 | C Calculate depths of centers and interfaces |
| 30 | zFace(1) = 0. _d 0 |
| 31 | DO K=1,Nz |
| 32 | dzF(K) = delZ(K) |
| 33 | zFace(K+1) = zFace(K)-delZ(K) |
| 34 | ENDDO |
| 35 | dzC(1) = delZ(1) * 0.5 _d 0 |
| 36 | zC(1) = delZ(1) * 0.5 _d 0 |
| 37 | DO K=2,Nz |
| 38 | dzC(K) = 0.5 _d 0 *(delZ(K-1)+delZ(K)) |
| 39 | zC(K) = zC(K) - dzC(K) |
| 40 | ENDDO |
| 41 | DO K=1,Nz |
| 42 | saFac(K) = 1. D0 |
| 43 | rDzC(K) = 1. _d 0/dzC(K) |
| 44 | rDzF(K) = 1. _d 0/dzF(K) |
| 45 | ENDDO |
| 46 | C |
| 47 | RETURN |
| 48 | END |
| ViewVC Help | |
| Powered by ViewVC 1.1.22 |