jscott/pkg_atm2d/atm2d_read_pickup.F

#include "ctrparam.h"
#include "ATM2D_OPTIONS.h"

C     !INTERFACE:
      SUBROUTINE ATM2D_READ_PICKUP( myThid )
C     *==========================================================*
C     | Read couple flux pickup, which contains fluxes for the   |
C     | initial ocean step and seaice growth/melt step. Called   |
C     | from the "normal" place (atm2d_init_vars).               |
C     *==========================================================*
        IMPLICIT NONE

#include "ATMSIZE.h"
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"

C     === Global SeaIce Variables ===
#include "THSICE_VARS.h"

C     === Atmos/Ocean/Seaice Interface Variables ===
#include "ATM2D_VARS.h"


C     !INPUT/OUTPUT PARAMETERS:
C     === Routine arguments ===
C     myThid - Thread no. that called this routine.
      INTEGER myThid

C     LOCAL VARIABLES:
      CHARACTER*(MAX_LEN_FNAM) fn
      INTEGER prec

      IF (startTime.NE.baseTime .OR. nIter0.NE.0) THEN

        prec=precFloat64

        WRITE(fn,'(A,I10.10)') 'pickup_atm2d.',nIter0

        CALL MDSREADFIELD(fn,prec,'RL',1,pass_slp, 1, myThid)
        CALL MDSREADFIELD(fn,prec,'RL',1,pass_qnet, 2, myThid)
        CALL MDSREADFIELD(fn,prec,'RL',1,pass_solarnet, 3, myThid)
        CALL MDSREADFIELD(fn,prec,'RL',1,pass_fu, 4, myThid)
        CALL MDSREADFIELD(fn,prec,'RL',1,pass_fv, 5, myThid)
        CALL MDSREADFIELD(fn,prec,'RL',1,pass_precip, 6, myThid)
        CALL MDSREADFIELD(fn,prec,'RL',1,pass_evap, 7, myThid)
        CALL MDSREADFIELD(fn,prec,'RL',1,pass_runoff, 8, myThid)
        CALL MDSREADFIELD(fn,prec,'RL',1,pass_wspeed, 9, myThid)
        CALL MDSREADFIELD(fn,prec,'RL',1,pass_pCO2, 10, myThid)
        CALL MDSREADFIELD(fn,prec,'RL',1,pass_sIceLoad, 11, myThid)

        CALL MDSREADFIELD(fn,prec,'RL',1,sHeating, 12, myThid)
        CALL MDSREADFIELD(fn,prec,'RL',1,flxCndBt, 13, myThid)
        CALL MDSREADFIELD(fn,prec,'RL',1,pass_prcAtm, 14, myThid)
        CALL MDSREADFIELD(fn,prec,'RL',1,snowPrc, 15, myThid)
        CALL MDSREADFIELD(fn,prec,'RL',1,icFrwAtm, 16, myThid)
        CALL MDSREADFIELD(fn,prec,'RL',1,icFlxSw   , 17, myThid)
        CALL MDSREADFIELD(fn,prec,'RL',1,siceAlb  , 18, myThid)

      ENDIF

      RETURN
      END

1	#include "ctrparam.h"
2	#include "ATM2D_OPTIONS.h"
3
4	C !INTERFACE:
5	SUBROUTINE ATM2D_READ_PICKUP( myThid )
6	C ==========================================================
7	C \| Read couple flux pickup, which contains fluxes for the \|
8	C \| initial ocean step and seaice growth/melt step. Called \|
9	C \| from the "normal" place (atm2d_init_vars). \|
10	C ==========================================================
11	IMPLICIT NONE
12
13	#include "ATMSIZE.h"
14	#include "SIZE.h"
15	#include "EEPARAMS.h"
16	#include "PARAMS.h"
17
18	C === Global SeaIce Variables ===
19	#include "THSICE_VARS.h"
20
21	C === Atmos/Ocean/Seaice Interface Variables ===
22	#include "ATM2D_VARS.h"
23
24
25	C !INPUT/OUTPUT PARAMETERS:
26	C === Routine arguments ===
27	C myThid - Thread no. that called this routine.
28	INTEGER myThid
29
30	C LOCAL VARIABLES:
31	CHARACTER*(MAX_LEN_FNAM) fn
32	INTEGER prec
33
34	IF (startTime.NE.baseTime .OR. nIter0.NE.0) THEN
35
36	prec=precFloat64
37
38	WRITE(fn,'(A,I10.10)') 'pickup_atm2d.',nIter0
39
40	CALL MDSREADFIELD(fn,prec,'RL',1,pass_slp, 1, myThid)
41	CALL MDSREADFIELD(fn,prec,'RL',1,pass_qnet, 2, myThid)
42	CALL MDSREADFIELD(fn,prec,'RL',1,pass_solarnet, 3, myThid)
43	CALL MDSREADFIELD(fn,prec,'RL',1,pass_fu, 4, myThid)
44	CALL MDSREADFIELD(fn,prec,'RL',1,pass_fv, 5, myThid)
45	CALL MDSREADFIELD(fn,prec,'RL',1,pass_precip, 6, myThid)
46	CALL MDSREADFIELD(fn,prec,'RL',1,pass_evap, 7, myThid)
47	CALL MDSREADFIELD(fn,prec,'RL',1,pass_runoff, 8, myThid)
48	CALL MDSREADFIELD(fn,prec,'RL',1,pass_wspeed, 9, myThid)
49	CALL MDSREADFIELD(fn,prec,'RL',1,pass_pCO2, 10, myThid)
50	CALL MDSREADFIELD(fn,prec,'RL',1,pass_sIceLoad, 11, myThid)
51
52	CALL MDSREADFIELD(fn,prec,'RL',1,sHeating, 12, myThid)
53	CALL MDSREADFIELD(fn,prec,'RL',1,flxCndBt, 13, myThid)
54	CALL MDSREADFIELD(fn,prec,'RL',1,pass_prcAtm, 14, myThid)
55	CALL MDSREADFIELD(fn,prec,'RL',1,snowPrc, 15, myThid)
56	CALL MDSREADFIELD(fn,prec,'RL',1,icFrwAtm, 16, myThid)
57	CALL MDSREADFIELD(fn,prec,'RL',1,icFlxSw , 17, myThid)
58	CALL MDSREADFIELD(fn,prec,'RL',1,siceAlb , 18, myThid)
59
60	ENDIF
61
62	RETURN
63	END
64