/[MITgcm]/MITgcm/pkg/exf/exf_clim_readparms.F
ViewVC logotype

Annotation of /MITgcm/pkg/exf/exf_clim_readparms.F

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

Revision 1.4 - (hide annotations) (download)
Thu Oct 9 04:19:19 2003 UTC (20 years, 7 months ago) by edhill
Branch: MAIN
CVS Tags: checkpoint51k_post, checkpoint51i_post, checkpoint51l_pre
Branch point for: tg2-branch
Changes since 1.3: +2 -2 lines 
 o first check-in for the "branch-genmake2" merge
 o verification suite as run on shelley (gcc 3.2.2):

Wed Oct  8 23:42:29 EDT 2003
                T           S           U           V
G D M    c        m  s        m  s        m  s        m  s
E p a R  g  m  m  e  .  m  m  e  .  m  m  e  .  m  m  e  .
N n k u  2  i  a  a  d  i  a  a  d  i  a  a  d  i  a  a  d
2 d e n  d  n  x  n  .  n  x  n  .  n  x  n  .  n  x  n  .

OPTFILE=NONE

Y Y Y Y 13 16 16 16  0 16 16 16 16 16 16 16 16 13 12  0  0 pass  adjustment.128x64x1
Y Y Y Y 16 16 16 16  0 16 16 16 16 16 16  0  0 16 16  0  0 pass  adjustment.cs-32x32x1
Y Y Y Y 16 16 16 16  0 16 16 16 16 16 16 22  0 16 16 22  0 pass  adjust_nlfs.cs-32x32x1
Y Y Y Y -- 13 13 16 16 13 13 13 13 16 16 16 16 16 16 16 16 N/O   advect_cs
Y Y Y Y -- 22 16 16 16 16 16 16 13 16 16 16 16 16 16 16 16 N/O   advect_xy
Y Y Y Y -- 13 16 13 16 16 16 16 16 16 16 22 16 16 16 16 16 N/O   advect_xz
Y Y Y Y 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 pass  aim.5l_cs
Y Y Y Y 14 16 16 16 16 16 16 16 16 13 16 16 16 16 16 13 16 pass  aim.5l_Equatorial_Channel
Y Y Y Y 16 16 16 16 16 16 16 16 16 16 16 13 16 16 13 13 16 pass  aim.5l_LatLon
Y Y Y Y 13 16 16 16 16 16 16 16 16 16 13 12 13 13 16 13 16 pass  exp0
Y Y Y Y 14 16 16 16 16 16 16 16 22 16 16 16 13 16 16 22 16 pass  exp1
Y Y Y Y 13 13 16 13 16 16 16 16 16 13 13 16 16 13 13 13 13 pass  exp2
Y Y Y Y 16 16 16 16 16 16 16 16 22 16 16 16 16 16 16 16 16 pass  exp4
Y Y Y Y 16 16 16 16 16 16 16 16 16 16 16 22 16 16 16 22 16 pass  exp5
Y Y Y Y 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 pass  front_relax
Y Y Y Y 14 16 16 13 13 16 16 13 13 16 13 13 16 12 13 13 16 pass  global_ocean.90x40x15
Y Y Y Y 10 16 16 13 13 16 13 16 16 13 13 13 13 16 16 13 16 FAIL  global_ocean.cs32x15
Y Y Y Y  6 11 12 13 13 12 13 16 13  9  9  9  9 10  9  9 11 FAIL  global_ocean_pressure
Y Y Y Y 14 16 16 13 16 16 16 13 13 13 13 13 16 12 16 13 16 pass  global_with_exf
Y Y Y Y 14 16 16 16 16 16 16 16 16 11 13 22 13 16 16  9 16 pass  hs94.128x64x5
Y Y Y Y 13 16 16 16 16 16 16 16 16 11 16 16 16 13 16 22 13 pass  hs94.1x64x5
Y Y Y Y 14 16 16 16 16 16 16 16 16 13 16 13 13 16 16 22 13 pass  hs94.cs-32x32x5
Y Y Y Y 10 10 16 13 13 16 16 16 22 16 13 13 13 13 13 22 13 FAIL  ideal_2D_oce
Y Y Y Y  8 16 16 16 16 16 16 16 16 13 13  8 16 16 16 16 16 FAIL  internal_wave
Y Y Y Y 14 16 16 16 16 16 16 16 16 13 13 22 13 13 13 22 16 pass  inverted_barometer
Y Y Y Y 12 16 16 16 16 16 16 16 16 16 13 12 13 13 13 13 13 FAIL  lab_sea
Y Y Y Y 11 16 16 16 16 16 16 16 13 13 13 12 13 16 13 12 13 FAIL  natl_box
Y Y Y Y 16 16 16 16 16 16 16 16 22 16 16 16 16 16 16 16 16 pass  plume_on_slope
Y Y Y Y 13 16 16 16 16 13 16 16 16 16 16 16 16 13 16 16 16 pass  solid-body.cs-32x32x1
1 edhill 1.4 c $Header: /u/u3/gcmpack/MITgcm/pkg/exf/exf_clim_readparms.F,v 1.3.2.1 2003/10/02 18:30:07 adcroft Exp $
2 heimbach 1.1
3 edhill 1.4 #include "EXF_OPTIONS.h"
4 heimbach 1.1
5    
6     subroutine exf_clim_readparms(
7     I mythid
8     & )
9    
10     c ==================================================================
11     c SUBROUTINE exf_clim_readparms
12     c ==================================================================
13     c
14     c o This routine initialises the climatologic forcing
15     c
16     c started: Ralf.Giering@FastOpt.de 25-Mai-20000
17     c
18     c ==================================================================
19     c SUBROUTINE exf_clim_readparms
20     c ==================================================================
21    
22     implicit none
23    
24     c == global variables ==
25    
26     #include "EEPARAMS.h"
27     #include "SIZE.h"
28     #include "cal.h"
29     #include "exf.h"
30 dimitri 1.3 #include "exf_param.h"
31 heimbach 1.1 #include "exf_clim_param.h"
32    
33     c == routine arguments ==
34    
35     integer mythid
36    
37     c == local variables ==
38    
39     integer i
40    
41     c == end of interface ==
42    
43     c Surface flux data.
44     namelist /exf_clim_nml/
45     & climtempstartdate1, climtempstartdate2, climtempperiod,
46     & climsaltstartdate1, climsaltstartdate2, climsaltperiod,
47 heimbach 1.2 & climsststartdate1, climsststartdate2, climsstperiod,
48     & climsssstartdate1, climsssstartdate2, climsssperiod,
49     & climtempfile, climsaltfile, climsstfile,
50     & climsssfile, climsstconst, climsssconst,
51 heimbach 1.1 & exf_clim_iprec, exf_clim_yftype
52 dimitri 1.3 #ifdef USE_EXF_INTERPOLATION
53     & ,climsst_lon0, climsst_lon_inc,
54     & climsst_lat0, climsst_lat_inc,
55     & climsst_nlon, climsst_nlat,
56     & climsss_lon0, climsss_lon_inc,
57     & climsss_lat0, climsss_lat_inc,
58     & climsss_nlon, climsss_nlat
59     #endif
60 heimbach 1.1
61     _BEGIN_MASTER(mythid)
62    
63     c Set default values.
64    
65     c Calendar data.
66     climtempstartdate1 = 0
67     climtempstartdate2 = 0
68     climtempperiod = 0
69    
70     climsaltstartdate1 = 0
71     climsaltstartdate2 = 0
72     climsaltperiod = 0
73    
74     climsststartdate1 = 0
75     climsststartdate2 = 0
76     climsstperiod = 0
77    
78     climsssstartdate1 = 0
79     climsssstartdate2 = 0
80     climsssperiod = 0
81    
82     c Data files.
83     climtempfile = ' '
84     climsaltfile = ' '
85     climsstfile = ' '
86     climsssfile = ' '
87    
88 heimbach 1.2 c Initialise the date arrays.
89 heimbach 1.1 do i = 1,4
90     climtempstartdate(i) = 0
91     climsaltstartdate(i) = 0
92     climsststartdate(i) = 0
93     climsssstartdate(i) = 0
94     enddo
95 heimbach 1.2
96     c Initialise constant values for relax. to constant SST, SSS
97     climsstconst = 0. _d 0
98     climsssconst = 0. _d 0
99 heimbach 1.1
100     c Initialise file type and field precision
101     exf_clim_iprec = 32
102     exf_clim_yftype = 'RL'
103    
104     c Check for the availability of the right calendar version.
105     if ( calendarversion .ne. usescalendarversion ) then
106     print*,' exf_Init: You are not using the appropriate'
107     print*,' version of the calendar package.'
108     print*
109     print*,' Please use Calendar version: ',
110     & usescalendarversion
111     stop ' stopped in exf_Init.'
112     endif
113    
114     c Next, read the forcing data file.
115     call nml_filter( 'data.exf_clim', scrunit1, myThid )
116     if (scrunit1 .eq. 0) then
117     stop 'exf_clim_readparms: reading namelist failed'
118     end if
119     read( scrunit1, nml = exf_clim_nml )
120     close( scrunit1 )
121    
122     c Complete the start date specifications for the forcing
123     c fields to get a complete calendar date array.
124    
125     c check for consistency
126    
127     if (.NOT. (exf_clim_iprec .EQ. 32
128     & .OR. exf_clim_iprec .EQ. 64)) then
129     stop 'stopped in exf_clim_readparms: value of iprec not allowed'
130     else if
131     & (.NOT. (exf_clim_yftype .EQ. 'RS'
132     & .OR. exf_clim_yftype .EQ. 'RL'))
133     & then
134     stop 'stopped in exf_clim_readparms: value of yftype not allowed'
135     end if
136    
137     #ifdef ALLOW_CLIMTEMP_RELAXATION
138     call cal_FullDate( climtempstartdate1, climtempstartdate2,
139     & climtempstartdate, mythid )
140     #endif
141    
142     #ifdef ALLOW_CLIMSALT_RELAXATION
143     call cal_FullDate( climsaltstartdate1, climsaltstartdate2,
144     & climsaltstartdate, mythid )
145     #endif
146    
147     #ifdef ALLOW_CLIMSST_RELAXATION
148     call cal_FullDate( climsststartdate1, climsststartdate2,
149     & climsststartdate, mythid )
150     #endif
151    
152     #ifdef ALLOW_CLIMSSS_RELAXATION
153     call cal_FullDate( climsssstartdate1, climsssstartdate2,
154     & climsssstartdate, mythid )
155     #endif
156    
157     _END_MASTER( mythid )
158    
159     _BARRIER
160    
161     end
  ViewVC Help
Powered by ViewVC 1.1.22