/[MITgcm]/MITgcm/model/src/packages_boot.F
ViewVC logotype

Contents of /MITgcm/model/src/packages_boot.F

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


Revision 1.11 - (show annotations) (download)
Fri Nov 15 19:58:21 2002 UTC (21 years, 6 months ago) by cheisey
Branch: MAIN
Changes since 1.10: +18 -3 lines
Two packages:  bulk_forcing (Bulk forcing)
and thermodynamic_seaice (adapted from LANL CICE.v2.0.2)
Earlier integration from Stephaine Dutkiewicz
and Patrick Heimbach.

Two ifdef statements for compile time,
ALLOW_TSEAICE and ALLOW_BULKFORMULA

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

cat data.pkg
# Packages
 &PACKAGES
 useBulkf=.TRUE.,
 useTSeaIce=.TRUE.,
 &

The bulk package requires an additional parameter file
with two namelists.

 cat data.blk
 &BULKF_PARM01
 RainFile=       'ncep_precip_m_cubed.bin',
 SolarFile=      'ncep_downsolr_cubed.bin',
 AirTempFile=    'ncep_tair_cubed.bin',
 AirhumidityFile='ncep_qair_g_cubed.bin',
 LongwaveFile=   'ncep_netlw_cubed.bin',
 UWindFile=      'ncep_uwind_cubed.bin',
 VWindFile=      'ncep_vwind_cubed.bin',
 WspeedFile=    ' ',
 RunoffFile=    ' ',
 QnetFile=       ' ',
 EmPFile=        'ncep_emp_calc_cubed.bin',
 CloudFile=      'ncep_totalcloud_cubed.bin',
 &

 &BULKF_PARM02
 qnet_off=0.0,
 empmr_off=0.0,
 conservcycle=311040000.,
 &



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-----------------------------------------------------------------------

WARNING:  useSEAICE and useTSEAICE are mutually exclusive.

todo: thermodynamic.F should be reviewed and cleaned up a bit.

1 C $Header: /u/u0/gcmpack/MITgcm/model/src/packages_boot.F,v 1.10 2002/11/12 20:45:41 heimbach Exp $
2 C $Name: $
3
4 #include "CPP_OPTIONS.h"
5
6 CBOP
7 C !ROUTINE: PACKAGES_BOOT
8 C !INTERFACE:
9 SUBROUTINE PACKAGES_BOOT( myThid )
10
11 C !DESCRIPTION: \bv
12 C *==========================================================*
13 C | SUBROUTINE PACKAGES_BOOT
14 C | o Routine to parse runtime package selection file
15 C *==========================================================*
16 C | Routine reads a flat file which contains a single flag
17 C | for each known "package". Flag can enables runtime
18 C | activation or deactivation of a package. The data is read
19 C | from a file called data.pkg.
20 C | Note - If a runtime activation flag is set on but the
21 C | package code was not included in the build then
22 C | the default behavior is to stop with an error.
23 C *==========================================================*
24 C \ev
25
26 C !USES:
27 IMPLICIT NONE
28 C === Global variables ===
29 #include "SIZE.h"
30 #include "EEPARAMS.h"
31 #include "PARAMS.h"
32
33 C !INPUT/OUTPUT PARAMETERS:
34 C === Routine arguments ===
35 C myThid - Number of this instances
36 INTEGER myThid
37
38 C !LOCAL VARIABLES:
39 C === Local variables ===
40 C msgBuf - Informational/error meesage buffer
41 C iUnit - Work variable for IO unit number
42 CHARACTER*(MAX_LEN_MBUF) msgBuf
43 INTEGER iUnit
44
45 C data.pkg namelists
46 C data.pkg namelists
47 cswdblk -- add useBulkf to NAMELIST
48 cswdice -- addd useTSeaIce to NAMELIST
49
50 NAMELIST /PACKAGES/
51 & useKPP,
52 & useGMRedi,
53 & useOBCS,
54 & useAIM,
55 & useGrdchk,
56 & useECCO,
57 & useSHAP_FILT,
58 & useZONAL_FILT,
59 & useFLT,
60 & usePTRACERS,
61 & useSEAICE,
62 & useBulkf,
63 & useTSEAICE
64
65
66 CEOP
67
68 _BEGIN_MASTER(myThid)
69
70 WRITE(msgBuf,'(A)') ' PACKAGES_BOOT: opening data.pkg'
71 CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
72 & SQUEEZE_RIGHT , 1)
73
74
75 CALL OPEN_COPY_DATA_FILE(
76 I 'data.pkg', 'PACKAGES_BOOT',
77 O iUnit,
78 I myThid )
79
80 C-- Default package configuration
81 useKPP=.FALSE.
82 useGMRedi=.FALSE.
83 useOBCS=.FALSE.
84 useAIM=.FALSE.
85 useGrdchk=.FALSE.
86 useECCO=.FALSE.
87 useSHAP_FILT=.FALSE.
88 useZONAL_FILT=.FALSE.
89 useFLT=.FALSE.
90 usePTRACERS=.FALSE.
91 useSeaIce=.FALSE.
92 cswdblk -- add
93 useBulkf=.FALSE.
94 cswdblk -- end add ---
95 cswdice --- add ---
96 useTSeaIce=.FALSE.
97 cswdice --- end add ---
98
99
100 C-- Read parameters from open data file
101 READ(UNIT=iUnit,NML=PACKAGES)
102
103
104 WRITE(msgBuf,'(A)') ' PACKAGES_BOOT: finished reading data.pkg'
105 CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
106 & SQUEEZE_RIGHT , 1)
107
108 C-- Close the open data file
109 CLOSE(iUnit)
110 _END_MASTER(myThid)
111
112 C-- Everyone else must wait for the parameters to be loaded
113 _BARRIER
114
115 return
116 end

  ViewVC Help
Powered by ViewVC 1.1.22