/[MITgcm]/MITgcm/pkg/bulk_force/bulkf_readparms.F
ViewVC logotype

Annotation of /MITgcm/pkg/bulk_force/bulkf_readparms.F

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

Revision 1.1 - (hide annotations) (download)
Thu Nov 21 19:11:42 2002 UTC (21 years, 7 months ago) by cheisey
Branch: MAIN
CVS Tags: checkpoint47e_post, checkpoint47c_post, checkpoint50c_post, checkpoint48e_post, checkpoint50c_pre, checkpoint48i_post, checkpoint50d_pre, checkpoint51, checkpoint50, checkpoint50d_post, checkpoint50b_pre, checkpoint51f_post, checkpoint48b_post, checkpoint51d_post, checkpoint48c_pre, checkpoint47d_pre, checkpoint47a_post, checkpoint48d_pre, checkpoint51j_post, checkpoint47i_post, checkpoint47d_post, checkpoint48d_post, checkpoint48f_post, checkpoint48h_post, checkpoint51b_pre, checkpoint47g_post, checkpoint51h_pre, checkpoint48a_post, checkpoint50f_post, checkpoint50a_post, checkpoint50f_pre, checkpoint47j_post, branch-exfmods-tag, branchpoint-genmake2, checkpoint48c_post, checkpoint51b_post, checkpoint51c_post, checkpoint47b_post, checkpoint50g_post, checkpoint50h_post, checkpoint50e_pre, checkpoint50i_post, checkpoint51i_pre, checkpoint47f_post, checkpoint50e_post, checkpoint51e_post, checkpoint48, checkpoint49, checkpoint51f_pre, checkpoint47h_post, checkpoint51g_post, checkpoint50b_post, checkpoint51a_post, checkpoint48g_post
Branch point for: branch-exfmods-curt, branch-genmake2 
Two packages:  bulk_force (Bulk forcing)
and therm_seaice (thermodynamic_seaice) - adopted from LANL CICE.v2.0.2
Earlier integration from Stephaine Dutkiewicz
and Patrick Heimbach.

Two ifdef statements for compile time,
ALLOW_THERM_SEAICE and ALLOW_BULK_FORCE

Two switches in data.pkg to turn on at run-time:

cat data.pkg
# Packages
 &PACKAGES
 useBulkForce=.TRUE.,
 useThermSeaIce=.TRUE.,
 &

WARNING:  useSEAICE and useThermSEAICE are mutually exclusive.

The bulk package requires an additional parameter file
with two namelists, data.ice and data.blk.

c ADAPTED FROM:
c LANL CICE.v2.0.2
c-----------------------------------------------------------------------
c.. thermodynamics (vertical physics) based on M. Winton 3-layer model
c.. See Bitz, C. M. and W. H. Lipscomb, 1999:  "An energy-conserving
c..       thermodynamic sea ice model for climate study."  J. Geophys.
c..       Res., 104, 15669 - 15677.
c..     Winton, M., 1999:  "A reformulated three-layer sea ice model."
c..       Submitted to J. Atmos. Ocean. Technol.

c.. authors Elizabeth C. Hunke and William Lipscomb
c..         Fluid Dynamics Group, Los Alamos National Laboratory
c-----------------------------------------------------------------------
1 cheisey 1.1
2     #include "CPP_OPTIONS.h"
3    
4     SUBROUTINE BULKF_READPARMS( myThid )
5     C /==========================================================\
6     C | SUBROUTINE BULKF_READPARMS |
7     C | o Routine to initialize BULKF variables and constants. |
8     C |==========================================================|
9     C | Initialize BULKF parameters, read in data.rbc |
10     C \==========================================================/
11     IMPLICIT NONE
12    
13     C === Global variables ===
14     #include "SIZE.h"
15     #include "EEPARAMS.h"
16     #include "PARAMS.h"
17     #include "GRID.h"
18     cswdblk
19     #ifdef ALLOW_BULK_FORCE
20     #include "BULKF.h"
21     #ifdef CONSERV_BULKF
22     #include "BULKF_CONSERV.h"
23     #endif
24     #endif
25     cswdblk
26    
27     C === Routine arguments ===
28     INTEGER myThid
29    
30     #ifdef ALLOW_BULK_FORCE
31    
32     C-- Bulk Formula parameter
33     NAMELIST /BULKF_PARM01/
34     & AirTempFile, AirHumidityFile, RainFile,
35     & SolarFile, LongwaveFile, UWindFile,
36     & VWindFile, RunoffFile, WSpeedFile, QnetFile,
37     & EmPFile, CloudFile, readwindstress, readsurface
38    
39     #ifdef CONSERV_BULKF
40     c- conserving qnet, empmr
41     NAMELIST /BULKF_PARM02/
42     & qnet_off, empmr_off, conservcycle
43     #endif
44    
45    
46     C === Local variables ===
47     C msgBuf - Informational/error meesage buffer
48     C iUnit - Work variable for IO unit number
49     CHARACTER*(MAX_LEN_MBUF) msgBuf
50     INTEGER iUnit
51    
52     _BEGIN_MASTER(myThid)
53    
54     WRITE(msgBuf,'(A)') ' BULKF_READPARMS: opening data.blk'
55     CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
56     & SQUEEZE_RIGHT , 1)
57    
58    
59     CALL OPEN_COPY_DATA_FILE(
60     I 'data.blk', 'BULKF_READPARMS',
61     O iUnit,
62     I myThid )
63    
64     C-- Default values
65     AirTempFile=' '
66     AirHumidityFile=' '
67     RainFile=' '
68     SolarFile=' '
69     LongwaveFile=' '
70     UWindFile=' '
71     VWindFile=' '
72     WspeedFile=' '
73     RunoffFile=' '
74     QnetFile=' '
75     EmPFile=' '
76     CloudFile=' '
77     SnowFile=' '
78     readwindstress=.TRUE.
79     readsurface=.TRUE.
80    
81     C-- Read parameters from open data file
82     READ(UNIT=iUnit,NML=BULKF_PARM01)
83    
84    
85     WRITE(msgBuf,'(A)') ' BULKF_READPARMS: read BULKF_PARM01'
86     CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
87     & SQUEEZE_RIGHT , 1)
88    
89     #ifdef CONSERV_BULKF
90     c -- default
91     qnet_off=0.d0
92     empmr_off=0.d0
93     READ(UNIT=iUnit,NML=BULKF_PARM02)
94     WRITE(msgBuf,'(A)') ' BULKF_READPARMS: read BULKF_PARM02'
95     CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
96     & SQUEEZE_RIGHT , 1)
97    
98     #endif /* CONSERV_BULKF */
99    
100    
101     C-- Close the open data file
102     CLOSE(iUnit)
103     _END_MASTER(myThid)
104    
105     C-- Everyone else must wait for the parameters to be loaded
106     _BARRIER
107    
108    
109     #endif /* ALLOW_BULK_FORCE */
110    
111     return
112     end
  ViewVC Help
Powered by ViewVC 1.1.22