/[MITgcm]/MITgcm/model/src/ini_grid.F
ViewVC logotype

Annotation of /MITgcm/model/src/ini_grid.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph

Revision 1.10.2.1 - (hide annotations) (download)
Tue Feb 26 16:04:47 2002 UTC (22 years, 2 months ago) by adcroft
Branch: release1
Changes since 1.10: +8 -2 lines 
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.10.2.1 C $Header: /u/gcmpack/MITgcm/model/src/ini_grid.F,v 1.11 2001/11/26 15:43:28 adcroft Exp $
2 adcroft 1.8 C $Name: $
3 cnh 1.1
4 cnh 1.5 #include "CPP_OPTIONS.h"
5 cnh 1.1
6 cnh 1.10 CBOP
7     C !ROUTINE: INI_GRID
8     C !INTERFACE:
9 cnh 1.1 SUBROUTINE INI_GRID( myThid )
10 cnh 1.10 C !DESCRIPTION: \bv
11     C *==========================================================*
12     C | SUBROUTINE INI_GRID
13     C | o Initialise model grid
14     C *==========================================================*
15     C | These arrays are used throughout the code in evaluating
16     C | gradients, integrals and spatial avarages. This routine
17     C | is called separately by each thread and initialise only
18     C | the region of the domain it is "responsible" for.
19     C | Notes:
20     C | Two examples are shown in this code. One illustrates the
21     C | initialisation of a cartesian grid. The other shows the
22     C | inialisation of a spherical polar grid. Other orthonormal
23     C | grids can be fitted into this design. In this case
24     C | custom metric terms also need adding to account for the
25     C | projections of velocity vectors onto these grids.
26     C | The structure used here also makes it possible to
27     C | implement less regular grid mappings. In particular
28     C | o Schemes which leave out blocks of the domain that are
29     C | all land could be supported.
30     C | o Multi-level schemes such as icosohedral or cubic
31     C | grid projectedions onto a sphere can also be fitted
32     C | within the strategy we use.
33     C | Both of the above also require modifying the support
34     C | routines that map computational blocks to simulation
35     C | domain blocks.
36     C *==========================================================*
37     C \ev
38    
39     C !USES:
40 adcroft 1.6 IMPLICIT NONE
41 cnh 1.1 C === Global variables ===
42     #include "SIZE.h"
43     #include "EEPARAMS.h"
44     #include "PARAMS.h"
45     #include "GRID.h"
46    
47 cnh 1.10 C !INPUT/OUTPUT PARAMETERS:
48 cnh 1.1 C == Routine arguments ==
49     C myThid - Number of this instance of INI_GRID
50     INTEGER myThid
51    
52 cnh 1.10 C !LOCAL VARIABLES:
53 cnh 1.1 C == Local variables ==
54     C msgBuf - Used for informational I/O.
55     CHARACTER*(MAX_LEN_MBUF) msgBuf
56 cnh 1.10 CEOP
57 cnh 1.1
58 adcroft 1.4 C-- Set up vertical grid and coordinate system
59     CALL INI_VERTICAL_GRID( myThid )
60    
61     C-- Set up horizontal grid and coordinate system
62 cnh 1.1 IF ( usingCartesianGrid ) THEN
63     CALL INI_CARTESIAN_GRID( myThid )
64     ELSEIF ( usingSphericalPolarGrid ) THEN
65     CALL INI_SPHERICAL_POLAR_GRID( myThid )
66 adcroft 1.8 ELSEIF ( usingCurvilinearGrid ) THEN
67     CALL INI_CURVILINEAR_GRID( myThid )
68 cnh 1.1 ELSE
69     _BEGIN_MASTER(myThid)
70     WRITE(msgBuf,'(A)')
71     & 'S/R INI_GRID: No grid coordinate system has been selected'
72     CALL PRINT_ERROR( msgBuf , myThid)
73     STOP 'ABNORMAL END: S/R INI_GRID'
74     _END_MASTER(myThid)
75     ENDIF
76 adcroft 1.8
77     C-- Write certain grid data to files (useful for creating netCDF
78     C and general post-analysis)
79     CALL WRITE_FLD_XY_RS( 'XC',' ',XC,0,myThid)
80     CALL WRITE_FLD_XY_RS( 'YC',' ',YC,0,myThid)
81     CALL WRITE_FLD_XY_RS( 'XG',' ',XG,0,myThid)
82     CALL WRITE_FLD_XY_RS( 'YG',' ',YG,0,myThid)
83 adcroft 1.10.2.1 CALL WRITE_FLD_XY_RS( 'RAC',' ',rA,0,myThid)
84     CALL WRITE_FLD_XY_RS( 'RAW',' ',rAw,0,myThid)
85     CALL WRITE_FLD_XY_RS( 'RAS',' ',rAs,0,myThid)
86     CALL WRITE_FLD_XY_RS( 'DXG',' ',DXG,0,myThid)
87     CALL WRITE_FLD_XY_RS( 'DYG',' ',DYG,0,myThid)
88     CALL WRITE_FLD_XY_RS( 'DXC',' ',DXC,0,myThid)
89     CALL WRITE_FLD_XY_RS( 'DYC',' ',DYC,0,myThid)
90 adcroft 1.9
91     C-- Print out statistics of each horizontal grid array (helps when debugging)
92     CALL MON_PRINTSTATS_RS(1,XC,'XC',myThid)
93     CALL MON_PRINTSTATS_RS(1,XG,'XG',myThid)
94     CALL MON_PRINTSTATS_RS(1,DXC,'DXC',myThid)
95     CALL MON_PRINTSTATS_RS(1,DXF,'DXF',myThid)
96     CALL MON_PRINTSTATS_RS(1,DXG,'DXG',myThid)
97     CALL MON_PRINTSTATS_RS(1,DXV,'DXV',myThid)
98     CALL MON_PRINTSTATS_RS(1,YC,'YC',myThid)
99     CALL MON_PRINTSTATS_RS(1,YG,'YG',myThid)
100     CALL MON_PRINTSTATS_RS(1,DYC,'DYC',myThid)
101     CALL MON_PRINTSTATS_RS(1,DYF,'DYF',myThid)
102     CALL MON_PRINTSTATS_RS(1,DYG,'DYG',myThid)
103     CALL MON_PRINTSTATS_RS(1,DYU,'DYU',myThid)
104     CALL MON_PRINTSTATS_RS(1,RA,'RA',myThid)
105     CALL MON_PRINTSTATS_RS(1,RAW,'RAW',myThid)
106     CALL MON_PRINTSTATS_RS(1,RAS,'RAS',myThid)
107     CALL MON_PRINTSTATS_RS(1,RAZ,'RAZ',myThid)
108 adcroft 1.8
109 cnh 1.1 RETURN
110     END
  ViewVC Help
Powered by ViewVC 1.1.22