51 |
EXTERNAL DIFFERENT_MULTIPLE |
EXTERNAL DIFFERENT_MULTIPLE |
52 |
|
|
53 |
C === Local arrays === |
C === Local arrays === |
|
COMMON /TDFIELDS/ |
|
|
& taux0, tauy0, Qnet0, EmPmR0, SST0, SSS0, Qsw0, |
|
|
& taux1, tauy1, Qnet1, EmPmR1, SST1, SSS1, Qsw1 |
|
|
_RS taux0 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
|
|
_RS tauy0 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
|
|
_RS Qnet0 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
|
|
_RS EmPmR0 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
|
|
_RS SST0 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
|
|
_RS SSS0 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
|
|
_RS Qsw0 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
|
|
_RS taux1 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
|
|
_RS tauy1 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
|
|
_RS Qnet1 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
|
|
_RS EmPmR1 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
|
|
_RS SST1 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
|
|
_RS SSS1 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
|
|
_RS Qsw1 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
|
54 |
|
|
55 |
C === Local variables === |
C === Local variables === |
56 |
INTEGER bi,bj,i,j,intime0,intime1 |
INTEGER bi,bj,i,j,intime0,intime1 |
62 |
IF ( periodicExternalForcing ) THEN |
IF ( periodicExternalForcing ) THEN |
63 |
|
|
64 |
C First call requires that we initialize everything to zero for safety |
C First call requires that we initialize everything to zero for safety |
65 |
IF ( myIter .EQ. nIter0 ) THEN |
cph transfered to ini_forcing |
66 |
CALL LEF_ZERO( taux0 ,myThid ) |
cph IF ( myIter .EQ. nIter0 ) THEN |
67 |
CALL LEF_ZERO( tauy0 ,myThid ) |
cph CALL LEF_ZERO( taux0 ,myThid ) |
68 |
CALL LEF_ZERO( Qnet0 ,myThid ) |
cph CALL LEF_ZERO( tauy0 ,myThid ) |
69 |
CALL LEF_ZERO( EmPmR0 ,myThid ) |
cph CALL LEF_ZERO( Qnet0 ,myThid ) |
70 |
CALL LEF_ZERO( SST0 ,myThid ) |
cph CALL LEF_ZERO( EmPmR0 ,myThid ) |
71 |
CALL LEF_ZERO( SSS0 ,myThid ) |
cph CALL LEF_ZERO( SST0 ,myThid ) |
72 |
CALL LEF_ZERO( Qsw0 ,myThid ) |
cph CALL LEF_ZERO( SSS0 ,myThid ) |
73 |
CALL LEF_ZERO( taux1 ,myThid ) |
cph CALL LEF_ZERO( Qsw0 ,myThid ) |
74 |
CALL LEF_ZERO( tauy1 ,myThid ) |
cph CALL LEF_ZERO( taux1 ,myThid ) |
75 |
CALL LEF_ZERO( Qnet1 ,myThid ) |
cph CALL LEF_ZERO( tauy1 ,myThid ) |
76 |
CALL LEF_ZERO( EmPmR1 ,myThid ) |
cph CALL LEF_ZERO( Qnet1 ,myThid ) |
77 |
CALL LEF_ZERO( SST1 ,myThid ) |
cph CALL LEF_ZERO( EmPmR1 ,myThid ) |
78 |
CALL LEF_ZERO( SSS1 ,myThid ) |
cph CALL LEF_ZERO( SST1 ,myThid ) |
79 |
CALL LEF_ZERO( Qsw1 ,myThid ) |
cph CALL LEF_ZERO( SSS1 ,myThid ) |
80 |
ENDIF |
cph CALL LEF_ZERO( Qsw1 ,myThid ) |
81 |
|
cph ENDIF |
82 |
|
|
83 |
C Now calculate whether it is time to update the forcing arrays |
C Now calculate whether it is time to update the forcing arrays |
84 |
rdt=1. _d 0 / deltaTclock |
rdt=1. _d 0 / deltaTclock |