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

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

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

Revision 1.1 - (show annotations) (download)
Thu Nov 21 19:11:42 2002 UTC (21 years, 5 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
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