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

Diff of /MITgcm/pkg/exf/exf_readparms.F

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

revision 1.2 by heimbach, Thu Feb 7 20:00:09 2002 UTC revision 1.4 by cheisey, Thu Dec 19 13:44:51 2002 UTC
# Line 1  Line 1 
1  c $Header$  c $Header$
2    
3  #include "EXF_CPPOPTIONS.h"  #include "EXF_CPPOPTIONS.h"
4    #ifdef ALLOW_OBCS
5    # include "OBCS_OPTIONS.h"
6    #endif
7    
8          subroutine exf_readparms( mythid )
       subroutine exf_readparms(  
      I                          mythid  
      &                        )  
9    
10  c     ==================================================================  c     ==================================================================
11  c     SUBROUTINE exf_readparms  c     SUBROUTINE exf_readparms
# Line 32  c              Patrick Heimbach, heimbac Line 32  c              Patrick Heimbach, heimbac
32  c              - added obcs parameters  c              - added obcs parameters
33  c     changed: Virginie Thierry, vthierry@ucsd.edu 04-June-2001  c     changed: Virginie Thierry, vthierry@ucsd.edu 04-June-2001
34  c              - added new obcs parameters (for each boundaries)  c              - added new obcs parameters (for each boundaries)
35    c     included runoff D. Stammer, Nov. 25, 2001
36    c     included pressure forcing. heimbach@mit.edu 05-Nov-2002
37  c  c
38    c     added "repeatPeriod" for cycling of forcing datasets 19-Dec-2002
39  c     ==================================================================  c     ==================================================================
40  c     SUBROUTINE exf_readparms  c     SUBROUTINE exf_readparms
41  c     ==================================================================  c     ==================================================================
# Line 47  c     == global variables == Line 50  c     == global variables ==
50  #include "exf.h"  #include "exf.h"
51  #include "exf_param.h"  #include "exf_param.h"
52  #include "exf_constants.h"  #include "exf_constants.h"
 #ifdef ALLOW_OBCS  
 # include "OBCS_OPTIONS.h"  
 #endif  
53    
54  c     == routine arguments ==  c     == routine arguments ==
55    
# Line 68  c     == end of interface == Line 68  c     == end of interface ==
68    
69  c     Surface flux data.  c     Surface flux data.
70        namelist /exf_nml/        namelist /exf_nml/
71         &    repeatPeriod,
72       &    hfluxstartdate1,    hfluxstartdate2,   hfluxperiod,       &    hfluxstartdate1,    hfluxstartdate2,   hfluxperiod,
73       &    atempstartdate1,    atempstartdate2,   atempperiod,       &    atempstartdate1,    atempstartdate2,   atempperiod,
74       &      aqhstartdate1,      aqhstartdate2,     aqhperiod,       &      aqhstartdate1,      aqhstartdate2,     aqhperiod,
75       &    sfluxstartdate1,    sfluxstartdate2,   sfluxperiod,       &    sfluxstartdate1,    sfluxstartdate2,   sfluxperiod,
76       &   precipstartdate1,   precipstartdate2,  precipperiod,       &   precipstartdate1,   precipstartdate2,  precipperiod,
77         &   runoffstartdate1,   runoffstartdate2,  runoffperiod,
78       &  ustressstartdate1,  ustressstartdate2, ustressperiod,       &  ustressstartdate1,  ustressstartdate2, ustressperiod,
79       &  vstressstartdate1,  vstressstartdate2, vstressperiod,       &  vstressstartdate1,  vstressstartdate2, vstressperiod,
80       &    uwindstartdate1,    uwindstartdate2,   uwindperiod,       &    uwindstartdate1,    uwindstartdate2,   uwindperiod,
# Line 83  c     Surface flux data. Line 85  c     Surface flux data.
85       &    obcsSstartdate1,    obcsSstartdate2,   obcsSperiod,       &    obcsSstartdate1,    obcsSstartdate2,   obcsSperiod,
86       &    obcsEstartdate1,    obcsEstartdate2,   obcsEperiod,       &    obcsEstartdate1,    obcsEstartdate2,   obcsEperiod,
87       &    obcsWstartdate1,    obcsWstartdate2,   obcsWperiod,       &    obcsWstartdate1,    obcsWstartdate2,   obcsWperiod,
88         &apressurestartdate1,apressurestartdate2,apressureperiod,
89       &          hfluxfile,          atempfile,       aqhfile,       &          hfluxfile,          atempfile,       aqhfile,
90       &         precipfile,          sfluxfile,   ustressfile,       &          sfluxfile,         precipfile,    runofffile,
91       &        vstressfile,          uwindfile,     vwindfile,       &        ustressfile,        vstressfile,
92       &         swfluxfile,         lwfluxfile,       &          uwindfile,          vwindfile,
93       &          exf_iprec,         exf_yftype       &         swfluxfile,         lwfluxfile, apressurefile,
94         &          exf_iprec,         exf_yftype,
95         &           scal_hfl,           scal_ust,       scal_vst,
96         &           scal_swf,           scal_sst,       scal_sss,
97         &     scal_apressure,           scal_prc,       scal_sfl
98    
99        _BEGIN_MASTER(mythid)        _BEGIN_MASTER(mythid)
100    
# Line 114  c     Calendar data. Line 121  c     Calendar data.
121        precipstartdate2   = 0        precipstartdate2   = 0
122        precipperiod       = 0.0 _d 0        precipperiod       = 0.0 _d 0
123    
124          runoffstartdate1   = 0
125          runoffstartdate2   = 0
126          runoffperiod       = 0.0 _d 0
127    
128        ustressstartdate1  = 0        ustressstartdate1  = 0
129        ustressstartdate2  = 0        ustressstartdate2  = 0
130        ustressperiod      = 0.0 _d 0        ustressperiod      = 0.0 _d 0
# Line 154  c     Calendar data. Line 165  c     Calendar data.
165        obcsWstartdate2    = 0        obcsWstartdate2    = 0
166        obcsWperiod        = 0.0 _d 0        obcsWperiod        = 0.0 _d 0
167    
168          apressurestartdate1    = 0
169          apressurestartdate2    = 0
170          apressureperiod        = 0.0 _d 0
171    
172          repeatPeriod           = 0.0 _d 0
173    
174    
175  c     Data files.  c     Data files.
176        hfluxfile          = ' '        hfluxfile          = ' '
177        atempfile          = ' '        atempfile          = ' '
178        aqhfile            = ' '        aqhfile            = ' '
179        precipfile         = ' '        precipfile         = ' '
180        sfluxfile          = ' '        sfluxfile          = ' '
181          runofffile         = ' '
182        ustressfile        = ' '        ustressfile        = ' '
183        vstressfile        = ' '        vstressfile        = ' '
184        uwindfile          = ' '        uwindfile          = ' '
185        vwindfile          = ' '        vwindfile          = ' '
186        swfluxfile         = ' '        swfluxfile         = ' '
187        lwfluxfile         = ' '        lwfluxfile         = ' '
188          apressurefile      = ' '
189    
190  c     Initialise the date arrays.  c     Initialise the date arrays.
191        do i = 1,4        do i = 1,4
# Line 174  c     Initialise the date arrays. Line 194  c     Initialise the date arrays.
194           aqhstartdate(i)      = 0           aqhstartdate(i)      = 0
195           precipstartdate(i)   = 0           precipstartdate(i)   = 0
196           sfluxstartdate(i)    = 0           sfluxstartdate(i)    = 0
197             runoffstartdate(i)   = 0
198           ustressstartdate(i)  = 0           ustressstartdate(i)  = 0
199           vstressstartdate(i)  = 0           vstressstartdate(i)  = 0
200           uwindstartdate(i)    = 0           uwindstartdate(i)    = 0
# Line 184  c     Initialise the date arrays. Line 205  c     Initialise the date arrays.
205           obcsSstartdate(i)    = 0           obcsSstartdate(i)    = 0
206           obcsEstartdate(i)    = 0           obcsEstartdate(i)    = 0
207           obcsWstartdate(i)    = 0           obcsWstartdate(i)    = 0
208             apressurestartdate(i)= 0
209        enddo        enddo
210    
211  c     Initialise file type and field precision  c     Initialise file type and field precision
212        exf_iprec       = 32        exf_iprec       = 32
213        exf_yftype      = 'RL'        exf_yftype      = 'RL'
214    
215    c     scaling between exf units and MITgcm units
216          scal_hfl       =  1. _d 0
217          scal_ust       =  1. _d 0
218          scal_vst       =  1. _d 0
219          scal_swf       =  1. _d 0
220          scal_sst       =  1. _d 0
221          scal_sss       =  1. _d 0
222          scal_apressure =  1. _d 0
223    #if (defined (ALLOW_BULKFORMULAE) && defined (ALLOW_ATM_TEMP))
224          scal_prc       =  1. _d 0
225    #else
226          scal_sfl       =  1. _d 0    
227    #endif
228    
229  c     Check for the availability of the right calendar version.  c     Check for the availability of the right calendar version.
230        if ( calendarversion .ne. usescalendarversion ) then        if ( calendarversion .ne. usescalendarversion ) then
231           print*,' exf_readparms: You are not using the appropriate'           print*,' exf_readparms: You are not using the appropriate'
# Line 224  c     check for consistency Line 260  c     check for consistency
260           stop 'stop in exf_readparms: value of exf_yftype not allowed'           stop 'stop in exf_readparms: value of exf_yftype not allowed'
261        end if        end if
262    
263    #ifdef ALLOW_RUNOFF
264          call cal_FullDate(  runoffstartdate1,  runoffstartdate2,
265         &                    runoffstartdate,            mythid )
266    #endif
267    
268  #ifdef ALLOW_BULKFORMULAE  #ifdef ALLOW_BULKFORMULAE
269    
270  #ifdef ALLOW_ATM_TEMP  #ifdef ALLOW_ATM_TEMP
# Line 296  c     check for consistency Line 337  c     check for consistency
337  #endif  #endif
338  #endif  #endif
339    
340    #ifdef ATMOSPHERIC_LOADING
341          call cal_FullDate(   apressurestartdate1, apressurestartdate2,
342         &                     apressurestartdate,          mythid )
343    #endif
344    
345        _END_MASTER( mythid )        _END_MASTER( mythid )
346    
347        _BARRIER        _BARRIER
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22