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

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

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

Revision 1.23 - (show annotations) (download)
Fri Nov 15 19:58:21 2002 UTC (21 years, 6 months ago) by cheisey
Branch: MAIN
Changes since 1.22: +34 -1 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/do_the_model_io.F,v 1.22 2001/09/26 18:09:14 cnh Exp $
2 C $Name: $
3
4 #include "CPP_OPTIONS.h"
5
6 CBOP
7 C !ROUTINE: DO_THE_MODEL_IO
8 C !INTERFACE:
9 SUBROUTINE DO_THE_MODEL_IO(myCurrentTime, myIter, myThid)
10 C !DESCRIPTION: \bv
11 C *==========================================================*
12 C | SUBROUTINE DO_THE_MODEL_IO
13 C | o Controlling routine for IO in model main time-stepping
14 C | loop.
15 C *==========================================================*
16 C | Many systems do not have thread safe IO so it is easier
17 C | to lump everything together and do dumping of fields
18 C | and updating of forcing terms in a single place.
19 C | The approach to IO used here is that writes are only
20 C | performed by thread 1 and that a process only writes out
21 C | its data ( it does not know about anyone elses data!)
22 C | Reading on the other hand is assumed to be from a file
23 C | containing all the data for all the processes. Only the
24 C | portion of data of interest to this process is actually
25 C | loaded. To work well this assumes the existence of some
26 C | reliable tool to join datasets together at the end of a
27 C | run - see joinds.p
28 C *==========================================================*
29 C \ev
30
31 C !USES:
32 IMPLICIT NONE
33 C == Global variables ===
34 #include "SIZE.h"
35 #include "EEPARAMS.h"
36 #include "PARAMS.h"
37 #include "DYNVARS.h"
38 cswdblk --- add ---
39 #ifdef ALLOW_BULKFORMULA
40 #include "BULKF.h"
41 #include "FFIELDS.h"
42 #endif
43 cswdblk -- end add ---
44 cswdice -- add ---
45 #ifdef ALLOW_TSEAICE
46 #include "ICE.h"
47 #include "ICE_DIAGS.h"
48 #endif
49 cswdice -- end add ---
50
51
52
53 LOGICAL DIFFERENT_MULTIPLE
54 EXTERNAL DIFFERENT_MULTIPLE
55
56 C !INPUT/OUTPUT PARAMETERS:
57 C == Routine arguments ==
58 C myThid - Thread number for this instance of the routine.
59 C myIter - Iteration number
60 C myCurrentTime - Current time of simulation ( s )
61 INTEGER myThid
62 INTEGER myIter
63 _RL myCurrentTime
64
65 CEOP
66
67 C-- Generaly only thread 1 does IO here. It can not start until
68 C-- all threads fields are ready.
69 IF (
70 & DIFFERENT_MULTIPLE(dumpFreq,myCurrentTime,
71 & myCurrentTime-deltaTClock)
72 & ) THEN
73
74 _BARRIER
75
76 C-- Write "text-plots" of certain fields
77 CALL PLOT_FIELD_XYZRL( uVel , 'Current uVel ' ,
78 & Nr, myIter, myThid )
79 CALL PLOT_FIELD_XYZRL( vVel , 'Current vVel ' ,
80 & Nr, myIter, myThid )
81 CALL PLOT_FIELD_XYZRL( theta, 'Current theta ' ,
82 & Nr, myIter, myThid )
83 CALL PLOT_FIELD_XYRL( etaN , 'Current etaN ' ,
84 & myIter, myThid )
85
86 ENDIF
87
88 C-- Write model state to binary file
89 CALL WRITE_STATE( myCurrentTime, myIter, myThid )
90
91 #ifdef ALLOW_TIMEAVE
92 C-- Do time averages
93 IF (taveFreq.GT.0.) THEN
94 CALL TIMEAVE_STATV_WRITE( myCurrentTime, myIter, myThid )
95 #ifdef ALLOW_AIM
96 IF (useAIM)
97 & CALL AIM_WRITE_DIAGS( myCurrentTime, myIter, myThid )
98 #endif
99 ENDIF
100 #endif
101
102 #ifdef ALLOW_GMREDI
103 C-- Do KPP diagnostics.
104 IF (useGMRedi)
105 & CALL GMREDI_DIAGS( myCurrentTime, myIter, myThid )
106 #endif
107
108 #ifdef ALLOW_KPP
109 C-- Do KPP diagnostics.
110 IF (useKPP)
111 & CALL KPP_DO_DIAGS( myCurrentTime, myIter, myThid )
112 #endif
113
114 cswdblk --- add ---
115 #ifdef ALLOW_BULKFORMULA
116 C-- Do bulkf diagnostics.
117 IF (useBulkf)
118 & CALL BULKF_DIAGS( myCurrentTime, myIter, myThid )
119 #endif
120 cswdblk -- end add ---
121
122 cswdice --- add ---
123 #ifdef ALLOW_TSEAICE
124 C-- Do seaice diagnostics.
125 IF (useSeaice)
126 & CALL ICE_DIAGS( myCurrentTime, myIter, myThid )
127 #endif
128 cswdice -- end add ---
129
130
131
132 RETURN
133 END
134
  ViewVC Help
Powered by ViewVC 1.1.22