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 |