jscott/pkg_atm2d/atm2d_init_vars.F

#include "ctrparam.h"
C
      SUBROUTINE ATM2D_INIT_VARS(myThid)
C     |==========================================================|
C     | Initialization steps which depend on pickups loading     |
C     \==========================================================/
      IMPLICIT NONE

#include "ATMSIZE.h"
#include "DRIVER.h"

C     !INPUT/OUTPUT PARAMETERS:
C     == Routine arguments ==
C     myThid - thread number for this instance of the routine.
      INTEGER myThid

C     === Local variables ===

      CALL atmosphere(0.0,0)

#ifdef CLM 
      CALL clm4mit2d
#endif

#ifdef ML_2D 
C    2D mixed layer model  initialization
      CALL mlogeom
      CALL read_ocean_clim
      CALL init_ocean
      CALL init_seaice
      CALL zonmeansocean
      CALL zonmeansice
      CALL mlocn2atm
      CALL seaice2atm
#  ifdef CPL_OCEANCO2
      CALL zerogary
#  endif
#endif
       
#ifdef CPL_TEM 
      PRINT *,' Start of tem_init'
      aocuav=0.0
      nepav=0.0
#  ifdef DATA4TEM
      nfile=1
#  endif
#  ifdef CPL_NEM
      PRINT *,' Reading climate data to run NEM during '
      PRINT *, '  TEM iInitialization'
      CALL rd_rstrt_nem
#  endif
c      call tem_init(inyear,lyear)
      CALL temclm_init
      CALL lcluc_init
      CALL updatelcluc(0)
      CALL tem_init
      adupt=0.0
      PRINT *,' TEM initialization finished'
#  ifdef CPL_NEM
      CALL tem2climate(0,11)
      PRINT *,'CH4'
      PRINT *,(temch4(j),j=1,jm0)
      PRINT *,'N2O'
      PRINT *,(temn2o(j),j=1,jm0)
#  endif
#endif

#ifdef OCEAN_3D
      CALL ATM2D_READ_PICKUP(myThid)
#endif

      RETURN
      END
1	jscott	1.1	#include "ctrparam.h"
2			C
3			SUBROUTINE ATM2D_INIT_VARS(myThid)
4			C \|==========================================================\|
5			C \| Initialization steps which depend on pickups loading \|
6			C \==========================================================/
7			IMPLICIT NONE
8
9			#include "ATMSIZE.h"
10			#include "DRIVER.h"
11
12			C !INPUT/OUTPUT PARAMETERS:
13			C == Routine arguments ==
14			C myThid - thread number for this instance of the routine.
15			INTEGER myThid
16
17			C === Local variables ===
18
19			CALL atmosphere(0.0,0)
20
21			#ifdef CLM
22			CALL clm4mit2d
23			#endif
24
25			#ifdef ML_2D
26			C 2D mixed layer model initialization
27			CALL mlogeom
28			CALL read_ocean_clim
29			CALL init_ocean
30			CALL init_seaice
31			CALL zonmeansocean
32			CALL zonmeansice
33			CALL mlocn2atm
34			CALL seaice2atm
35			# ifdef CPL_OCEANCO2
36			CALL zerogary
37			# endif
38			#endif
39
40			#ifdef CPL_TEM
41			PRINT *,' Start of tem_init'
42			aocuav=0.0
43			nepav=0.0
44			# ifdef DATA4TEM
45			nfile=1
46			# endif
47			# ifdef CPL_NEM
48			PRINT *,' Reading climate data to run NEM during '
49			PRINT *, ' TEM iInitialization'
50			CALL rd_rstrt_nem
51			# endif
52			c call tem_init(inyear,lyear)
53			CALL temclm_init
54			CALL lcluc_init
55			CALL updatelcluc(0)
56			CALL tem_init
57			adupt=0.0
58			PRINT *,' TEM initialization finished'
59			# ifdef CPL_NEM
60			CALL tem2climate(0,11)
61			PRINT *,'CH4'
62			PRINT *,(temch4(j),j=1,jm0)
63			PRINT *,'N2O'
64			PRINT *,(temn2o(j),j=1,jm0)
65			# endif
66			#endif
67
68			#ifdef OCEAN_3D
69			CALL ATM2D_READ_PICKUP(myThid)
70			#endif
71
72			RETURN
73			END