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 ===
      INTEGER j

      CALL atmosphere(0.0,0)   ! initialize the atmospheric model

#ifdef CLM 
      CALL clm4mit2d           ! initialize CLM
#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	#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	INTEGER j
19
20	CALL atmosphere(0.0,0) ! initialize the atmospheric model
21
22	#ifdef CLM
23	CALL clm4mit2d ! initialize CLM
24	#endif
25
26	#ifdef ML_2D
27	C 2D mixed layer model initialization
28	CALL mlogeom
29	CALL read_ocean_clim
30	CALL init_ocean
31	CALL init_seaice
32	CALL zonmeansocean
33	CALL zonmeansice
34	CALL mlocn2atm
35	CALL seaice2atm
36	# ifdef CPL_OCEANCO2
37	CALL zerogary
38	# endif
39	#endif
40
41	#ifdef CPL_TEM
42	PRINT *,' Start of tem_init'
43	aocuav=0.0
44	nepav=0.0
45	# ifdef DATA4TEM
46	nfile=1
47	# endif
48	# ifdef CPL_NEM
49	PRINT *,' Reading climate data to run NEM during '
50	PRINT *, ' TEM iInitialization'
51	CALL rd_rstrt_nem
52	# endif
53	c call tem_init(inyear,lyear)
54	CALL temclm_init
55	CALL lcluc_init
56	CALL updatelcluc(0)
57	CALL tem_init
58	adupt=0.0
59	PRINT *,' TEM initialization finished'
60	# ifdef CPL_NEM
61	CALL tem2climate(0,11)
62	PRINT *,'CH4'
63	PRINT *,(temch4(j),j=1,jm0)
64	PRINT *,'N2O'
65	PRINT *,(temn2o(j),j=1,jm0)
66	# endif
67	#endif
68
69	#ifdef OCEAN_3D
70	CALL ATM2D_READ_PICKUP(myThid)
71	#endif
72
73	RETURN
74	END