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