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

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

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

Revision 1.4 - (show 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 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
3 #include "EXF_OPTIONS.h"
4
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 #include "exf_param.h"
31 #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 & climsststartdate1, climsststartdate2, climsstperiod,
48 & climsssstartdate1, climsssstartdate2, climsssperiod,
49 & climtempfile, climsaltfile, climsstfile,
50 & climsssfile, climsstconst, climsssconst,
51 & exf_clim_iprec, exf_clim_yftype
52 #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
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 c Initialise the date arrays.
89 do i = 1,4
90 climtempstartdate(i) = 0
91 climsaltstartdate(i) = 0
92 climsststartdate(i) = 0
93 climsssstartdate(i) = 0
94 enddo
95
96 c Initialise constant values for relax. to constant SST, SSS
97 climsstconst = 0. _d 0
98 climsssconst = 0. _d 0
99
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