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 |