pkg/atm2d/atm2d_init_fixed.F

C $Header: /u/gcmpack/MITgcm/pkg/atm2d/atm2d_init_fixed.F,v 1.4 2011/06/02 21:33:21 jscott Exp $
C $Name:  $

#include "ctrparam.h"
#ifdef OCEAN_3D
#  include "ATM2D_OPTIONS.h"
#endif
C
#ifndef OCEAN_3D
      SUBROUTINE ATM2D_INIT_FIXED(nCouplePer, myThid)
#else
      SUBROUTINE ATM2D_INIT_FIXED(myThid)
#endif
C     |==========================================================|
C     | Initialization steps prior to any pickup info loaded.    |
C     \==========================================================/
      IMPLICIT NONE

#include "ATMSIZE.h"
#include "DRIVER.h"
#ifdef OCEAN_3D
#  include "SIZE.h"
#  include "EEPARAMS.h"
#  include "PARAMS.h"
#endif

C     !INPUT/OUTPUT PARAMETERS:
C     == Routine arguments ==
C     nCouplePer - total number of coupled periods for this run
C     myThid - thread number for this instance of the routine.
#ifndef OCEAN_3D
      INTEGER nCouplePer    !not a subroutine parm for OCEAN_3D
#endif
      INTEGER myThid

      NAMELIST /COUPLE_PARM/
     &        dtcouple, dtatm, dtocn, startYear, endYear, taveDump

C     === Local variables ===
#ifdef OCEAN_3D
      INTEGER nCouplePer    !not a subroutine parm for OCEAN_3D
# ifdef ATM2D_MPI_ON
      INTEGER msgUnit
# endif
#endif

      CALL CHECK_SETTINGS

C     set default values for these parms in couple.nml
      dtatm = 1
      dtocn = 8
      dtcouple = 8
      startYear = 1
      endYear = 100
      taveDump=100

C     Next lines done in stand-alone ML model, so don't use any
C     MITGCM helper routines; hopefully no unit conflict...
      OPEN(514,file='couple.nml',status='old')
      READ(514,COUPLE_PARM)
      CLOSE(514)

#ifdef OCEAN_3D
      CALL ATM2D_READPARMS(myThid)
#endif

      ncall_atm=dtcouple/dtatm
      ncall_ocean=dtcouple/dtocn
      nCouplePer = (endYear-startYear+1)*24*365/dtcouple

      PRINT *,'Model run years:',startYear,' to ', endYear
      PRINT *,'dcouple=',dtcouple
      PRINT *,'dtatm=',dtatm
      PRINT *,'dtocn=',dtocn
      PRINT *,'Total number of coupled periods:',nCouplePer

#ifdef OCEAN_3D
      nTimeSteps = nCouplePer    !overwrite whatever MITGCM gets from data
#endif

#ifdef CPL_OCEANCO2
      ocupt=0.0
      temuptann=0.0
#endif

#ifdef OCEAN_3D
#  ifdef ATM2D_MPI_ON

      CALL INITIALISE

C     Perform registration with other components
      CALL ACCEPT_COMPONENT_REGISTRATIONS

C     Coordinate the transfer configuration information
C     between components
      msgUnit = standardMessageUnit
      CALL EXCH_COMPONENT_CONFIGS(
     I                             msgUnit )
#  endif

      CALL INIT_ATM2D(dtatm, dtocn, dtcouple, myThid)
#endif

      RETURN
      END
1	C $Header: /u/gcmpack/MITgcm/pkg/atm2d/atm2d_init_fixed.F,v 1.4 2011/06/02 21:33:21 jscott Exp $
2	C $Name: $
3
4	#include "ctrparam.h"
5	#ifdef OCEAN_3D
6	# include "ATM2D_OPTIONS.h"
7	#endif
8	C
9	#ifndef OCEAN_3D
10	SUBROUTINE ATM2D_INIT_FIXED(nCouplePer, myThid)
11	#else
12	SUBROUTINE ATM2D_INIT_FIXED(myThid)
13	#endif
14	C \|==========================================================\|
15	C \| Initialization steps prior to any pickup info loaded. \|
16	C \==========================================================/
17	IMPLICIT NONE
18
19	#include "ATMSIZE.h"
20	#include "DRIVER.h"
21	#ifdef OCEAN_3D
22	# include "SIZE.h"
23	# include "EEPARAMS.h"
24	# include "PARAMS.h"
25	#endif
26
27	C !INPUT/OUTPUT PARAMETERS:
28	C == Routine arguments ==
29	C nCouplePer - total number of coupled periods for this run
30	C myThid - thread number for this instance of the routine.
31	#ifndef OCEAN_3D
32	INTEGER nCouplePer !not a subroutine parm for OCEAN_3D
33	#endif
34	INTEGER myThid
35
36	NAMELIST /COUPLE_PARM/
37	& dtcouple, dtatm, dtocn, startYear, endYear, taveDump
38
39	C === Local variables ===
40	#ifdef OCEAN_3D
41	INTEGER nCouplePer !not a subroutine parm for OCEAN_3D
42	# ifdef ATM2D_MPI_ON
43	INTEGER msgUnit
44	# endif
45	#endif
46
47	CALL CHECK_SETTINGS
48
49	C set default values for these parms in couple.nml
50	dtatm = 1
51	dtocn = 8
52	dtcouple = 8
53	startYear = 1
54	endYear = 100
55	taveDump=100
56
57	C Next lines done in stand-alone ML model, so don't use any
58	C MITGCM helper routines; hopefully no unit conflict...
59	OPEN(514,file='couple.nml',status='old')
60	READ(514,COUPLE_PARM)
61	CLOSE(514)
62
63	#ifdef OCEAN_3D
64	CALL ATM2D_READPARMS(myThid)
65	#endif
66
67	ncall_atm=dtcouple/dtatm
68	ncall_ocean=dtcouple/dtocn
69	nCouplePer = (endYear-startYear+1)24365/dtcouple
70
71	PRINT *,'Model run years:',startYear,' to ', endYear
72	PRINT *,'dcouple=',dtcouple
73	PRINT *,'dtatm=',dtatm
74	PRINT *,'dtocn=',dtocn
75	PRINT *,'Total number of coupled periods:',nCouplePer
76
77	#ifdef OCEAN_3D
78	nTimeSteps = nCouplePer !overwrite whatever MITGCM gets from data
79	#endif
80
81	#ifdef CPL_OCEANCO2
82	ocupt=0.0
83	temuptann=0.0
84	#endif
85
86	#ifdef OCEAN_3D
87	# ifdef ATM2D_MPI_ON
88
89	CALL INITIALISE
90
91	C Perform registration with other components
92	CALL ACCEPT_COMPONENT_REGISTRATIONS
93
94	C Coordinate the transfer configuration information
95	C between components
96	msgUnit = standardMessageUnit
97	CALL EXCH_COMPONENT_CONFIGS(
98	I msgUnit )
99	# endif
100
101	CALL INIT_ATM2D(dtatm, dtocn, dtcouple, myThid)
102	#endif
103
104	RETURN
105	END