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

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

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

Revision 1.14 - (hide 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 mlosch 1.14 C $Header: /u/gcmpack/MITgcm/model/src/ini_fields.F,v 1.13 2002/07/13 04:59:42 heimbach Exp $
2 adcroft 1.10 C $Name: $
3 cnh 1.1
4 cnh 1.4 #include "CPP_OPTIONS.h"
5 cnh 1.1
6 cnh 1.12 CBOP
7     C !ROUTINE: INI_FIELDS
8     C !INTERFACE:
9 cnh 1.1 SUBROUTINE INI_FIELDS( myThid )
10 cnh 1.12 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 adcroft 1.5 IMPLICIT NONE
24 cnh 1.1 C === Global variables ===
25     #include "SIZE.h"
26     #include "EEPARAMS.h"
27     #include "PARAMS.h"
28    
29 cnh 1.12 C !INPUT/OUTPUT PARAMETERS:
30 cnh 1.1 C == Routine arguments ==
31     C myThid - Number of this instance of INI_FIELDS
32     INTEGER myThid
33 cnh 1.12 CEOP
34 cnh 1.1
35     C-- Set model variables to initial/restart conditions
36 cnh 1.4 IF ( startTime .EQ. 0. .AND. nIter0 .EQ. 0 ) THEN
37 adcroft 1.10 CALL INI_VEL( myThid )
38 cnh 1.1 CALL INI_THETA( myThid )
39     CALL INI_SALT( myThid )
40 adcroft 1.6 CALL INI_PSURF( myThid )
41 mlosch 1.14 CALL INI_PRESSURE( myThid )
42 heimbach 1.13 #ifdef INCLUDE_EP_FORCING_CODE
43     CALL INI_EP( myThid )
44     #endif
45 adcroft 1.11 IF (nonHydrostatic) CALL INI_PNH( myThid )
46 cnh 1.1 ELSE
47 cnh 1.3 IF ( nIter0 .EQ. 0 ) THEN
48     nIter0 = INT(startTime/deltaTClock )
49     ENDIF
50     CALL READ_CHECKPOINT( nIter0, myThid )
51 cnh 1.1 ENDIF
52    
53     RETURN
54     END
  ViewVC Help
Powered by ViewVC 1.1.22