Parent Directory | Revision Log | Revision Graph
Merging changes on MAIN between checkpoint43 and checkpoint43a-release1mods Command: cvs -q update -jcheckpoint43 -jcheckpoint43a-release1mods -d -P These changes are most of the changes between c43 and c44 except those that occured after "12:45 11 Jan 2002". As far as I can tell it is checkpoint43 with the following mods: o fix bug in mom_vi_del2uv o select when filters are applied ; add options to zonal_filter (data.zonfilt) o gmredi: fix Pb in the adiabatic form ; add options (.e.g. Bolus advection) o update AIM experiments (NCEP input files) o improve and extend diagnostics (Monitor, TimeAve with NonLin-FrSurf) o added some stuff for AD o Jamar wet-points This update does not contain the following mods that are in checkpoint44 o bug fix in pkg/generic_advdiff/ - thread related bug, bi,bj arguments in vertical advection routines o some changes to pkg/autodiff, pkg/cost, pkg/exf, pkg/ecco, verification/carbon and model/src/ related to adjoint o some new Matlab scripts for diagnosing model density - utils/matlab/dens_poly3.m and ini_poly3.m The list of exclusions is accurate based on a "cvs diff". The list of inclusions is based on the record in doc/tag-index which may not be complete.
1 | adcroft | 1.1.4.1 | C $Header: /u/gcmpack/MITgcm/verification/aim.5l_cs/code/mom_vi_mask_vort3.F,v 1.1 2002/01/09 00:28:56 jmc Exp $ |
2 | C $Name: checkpoint43a-release1mods $ | ||
3 | jmc | 1.1 | |
4 | #include "CPP_OPTIONS.h" | ||
5 | |||
6 | SUBROUTINE MOM_VI_MASK_VORT3( | ||
7 | I bi,bj,k, | ||
8 | O hFacZ,r_hFacZ, | ||
9 | U vort3, | ||
10 | I myThid) | ||
11 | IMPLICIT NONE | ||
12 | C | ||
13 | C Compute hFactor and recip_hFactor at the corner (Z-point) | ||
14 | C and set vort3 to zero where hFac = 0 (apply a mask) | ||
15 | C | ||
16 | |||
17 | C == Global variables == | ||
18 | #include "SIZE.h" | ||
19 | #include "GRID.h" | ||
20 | |||
21 | C == Routine arguments == | ||
22 | INTEGER bi,bj,k | ||
23 | _RS hFacZ(1-OLx:sNx+OLx,1-OLy:sNy+OLy) | ||
24 | _RS r_hFacZ(1-OLx:sNx+OLx,1-OLy:sNy+OLy) | ||
25 | _RL vort3(1-OLx:sNx+OLx,1-OLy:sNy+OLy) | ||
26 | _RS hFacZOpen | ||
27 | INTEGER myThid | ||
28 | |||
29 | C == Local variables == | ||
30 | INTEGER I,J | ||
31 | |||
32 | C-- Calculate open water fraction at vorticity points | ||
33 | |||
34 | c DO i=1-Olx,sNx+Olx | ||
35 | c hFacZ(i,1-Oly)=0. | ||
36 | c r_hFacZ(i,1-Oly)=0. | ||
37 | c ENDDO | ||
38 | |||
39 | DO j=2-Oly,sNy+Oly | ||
40 | c hFacZ(1-Olx,j)=0. | ||
41 | c r_hFacZ(1-Olx,j)=0. | ||
42 | DO i=2-Olx,sNx+Olx | ||
43 | hFacZOpen=min(_hFacW(i,j,k,bi,bj), | ||
44 | & _hFacW(i,j-1,k,bi,bj)) | ||
45 | hFacZOpen=min(_hFacS(i,j,k,bi,bj),hFacZOpen) | ||
46 | hFacZOpen=min(_hFacS(i-1,j,k,bi,bj),hFacZOpen) | ||
47 | hFacZ(i,j)=hFacZOpen | ||
48 | IF (hFacZ(i,j).EQ.0.) THEN | ||
49 | r_hFacZ(i,j)=0. | ||
50 | vort3(i,j)=0. | ||
51 | ELSE | ||
52 | r_hFacZ(i,j)=1./hFacZ(i,j) | ||
53 | ENDIF | ||
54 | ENDDO | ||
55 | ENDDO | ||
56 | |||
57 | RETURN | ||
58 | END |
ViewVC Help | |
Powered by ViewVC 1.1.22 |