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

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

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

Revision 1.14 - (show annotations) (download)
Wed Sep 18 16:38:02 2002 UTC (21 years, 8 months ago) by mlosch
Branch: MAIN
CVS Tags: checkpoint46n_post, checkpoint47e_post, checkpoint46l_post, checkpoint47c_post, checkpoint50c_post, checkpoint48e_post, checkpoint50c_pre, checkpoint48i_post, checkpoint46l_pre, checkpoint50, checkpoint50d_post, checkpoint50b_pre, checkpoint48b_post, checkpoint48c_pre, checkpoint47d_pre, checkpoint47a_post, checkpoint48d_pre, checkpoint47i_post, checkpoint47d_post, checkpoint48d_post, checkpoint48f_post, checkpoint46j_pre, checkpoint48h_post, checkpoint47g_post, checkpoint46j_post, checkpoint46k_post, checkpoint48a_post, checkpoint50f_post, checkpoint50a_post, checkpoint50f_pre, checkpoint47j_post, branch-exfmods-tag, checkpoint48c_post, checkpoint47b_post, checkpoint46h_pre, checkpoint46m_post, checkpoint50g_post, checkpoint46g_post, checkpoint50e_pre, checkpoint47f_post, checkpoint50e_post, checkpoint46i_post, checkpoint50d_pre, checkpoint47, checkpoint48, checkpoint49, checkpoint46h_post, checkpoint48g_post, checkpoint47h_post, checkpoint50b_post
Branch point for: branch-exfmods-curt
Changes since 1.13: +2 -1 lines 
o Include a new diagnostic variable phiHydLow for the ocean model
  - in z-coordinates, it is the bottom pressure anomaly
  - in p-coordinates, it is the sea surface elevation
  - in both cases, these variable have global drift, reflecting the mass
    drift in z-coordinates and the volume drift in p-coordinates
  - included time averaging for phiHydLow, be aware of the drift!
o depth-dependent computation of Bo_surf for pressure coordinates
  in the ocean (buoyancyRelation='OCEANICP')
  - requires a new routine (FIND_RHO_SCALAR) to compute density with only
    Theta, Salinity, and Pressure in the parameter list. This routine is
    presently contained in find_rho.F. This routine does not give the
    correct density for 'POLY3', which would be a z-dependent reference
    density.
o cleaned up find_rho
  - removed obsolete 'eqn' from the parameter list.
o added two new verification experiments: gop and goz
  (4x4 degree global ocean, 15 layers in pressure and height coordinates)
1 C $Header: /u/gcmpack/MITgcm/model/src/ini_fields.F,v 1.13 2002/07/13 04:59:42 heimbach Exp $
2 C $Name: $
3
4 #include "CPP_OPTIONS.h"
5
6 CBOP
7 C !ROUTINE: INI_FIELDS
8 C !INTERFACE:
9 SUBROUTINE INI_FIELDS( myThid )
10 C !DESCRIPTION: \bv
11 C *==========================================================*
12 C | SUBROUTINE INI_FIELDS
13 C | o Initialise model fields
14 C *==========================================================*
15 C | Initial conditions of the model state variables are
16 C | set here. For continuation runs a restart file will be
17 C | used. For original runs custom rules can be included
18 C | here.
19 C *==========================================================*
20 C \ev
21
22 C !USES:
23 IMPLICIT NONE
24 C === Global variables ===
25 #include "SIZE.h"
26 #include "EEPARAMS.h"
27 #include "PARAMS.h"
28
29 C !INPUT/OUTPUT PARAMETERS:
30 C == Routine arguments ==
31 C myThid - Number of this instance of INI_FIELDS
32 INTEGER myThid
33 CEOP
34
35 C-- Set model variables to initial/restart conditions
36 IF ( startTime .EQ. 0. .AND. nIter0 .EQ. 0 ) THEN
37 CALL INI_VEL( myThid )
38 CALL INI_THETA( myThid )
39 CALL INI_SALT( myThid )
40 CALL INI_PSURF( myThid )
41 CALL INI_PRESSURE( myThid )
42 #ifdef INCLUDE_EP_FORCING_CODE
43 CALL INI_EP( myThid )
44 #endif
45 IF (nonHydrostatic) CALL INI_PNH( myThid )
46 ELSE
47 IF ( nIter0 .EQ. 0 ) THEN
48 nIter0 = INT(startTime/deltaTClock )
49 ENDIF
50 CALL READ_CHECKPOINT( nIter0, myThid )
51 ENDIF
52
53 RETURN
54 END
  ViewVC Help
Powered by ViewVC 1.1.22