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 | adcroft | 1.2 | C $Header: /u/gcmpack/models/MITgcmUV/model/src/ini_vertical_grid.F,v 1.1 1998/07/02 14:17:11 adcroft Exp $ |
2 | adcroft | 1.1 | |
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 | adcroft | 1.2 | zFace(K+1) = zFace(K)-delZ(K) |
34 | adcroft | 1.1 | 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 | adcroft | 1.2 | zC(K) = zC(K) - dzC(K) |
40 | adcroft | 1.1 | 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 |