35 |
#include "GRID.h" |
#include "GRID.h" |
36 |
#include "OCNCPL.h" |
#include "OCNCPL.h" |
37 |
C == Routine arguments == |
C == Routine arguments == |
38 |
_RL atm_HeatFlux( 1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RL atm_HeatFlux( 1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
39 |
_RL atm_TauX( 1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RL atm_TauX( 1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
40 |
_RL atm_TauY( 1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RL atm_TauY( 1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
41 |
_RL atm_Qlatent( 1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RL atm_Qlatent( 1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
42 |
_RL atm_Qsensible( 1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RL atm_Qsensible( 1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
43 |
_RL atm_Qlongwave( 1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RL atm_Qlongwave( 1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
44 |
_RL atm_Qshortwave( 1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RL atm_Qshortwave( 1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
45 |
_RL atm_uVelGround( 1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RL atm_uVelGround( 1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
46 |
_RL atm_vVelGround( 1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RL atm_vVelGround( 1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
47 |
_RL atm_FWFlux( 1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RL atm_FWFlux( 1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
48 |
_RL atm_Hatm( 1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RL atm_Hatm( 1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
49 |
_RL ocn_SSTocn( 1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RL ocn_SSTocn( 1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
50 |
_RL ocn_Hocn( 1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RL ocn_Hocn( 1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
51 |
C == Local variables == |
C == Local variables == |
52 |
INTEGER iLoop |
INTEGER iLoop |
53 |
INTEGER myThid |
INTEGER myThid |
54 |
INTEGER myCurrentIter |
INTEGER myCurrentIter |
55 |
_RL myCurrentTime |
_RL myCurrentTime |
56 |
|
INTEGER I,J,K,bi,bj |
57 |
CEOP |
CEOP |
58 |
|
|
59 |
CALL COMP_OCN_PUSH_CONTEXT(1) |
CALL COMP_OCN_PUSH_CONTEXT(1) |
60 |
|
|
61 |
|
myThid = 1 |
62 |
|
CALL EEBOOT() |
63 |
C Read the import fields |
C Read the import fields |
64 |
DO bj=myByLo(myThid),myByHi(myThid) |
DO bj=myByLo(myThid),myByHi(myThid) |
65 |
DO bi=myBxLo(myThid),myBxHi(myThid) |
DO bi=myBxLo(myThid),myBxHi(myThid) |
66 |
DO j=1-OLy,sNy+OLy |
DO j=1-OLy,sNy+OLy |
67 |
DO i=1-OLx,sNx+OLx |
DO i=1-OLx,sNx+OLx |
68 |
HeatFlux( i,j,bi,bj)=atm_HeatFlux( i,j,bi,bj) |
HeatFlux( i,j,bi,bj)=atm_HeatFlux( i,j) |
69 |
TauX( i,j,bi,bj)=atm_tauX( i,j,bi,bj) |
TauX( i,j,bi,bj)=atm_tauX( i,j) |
70 |
TauY( i,j,bi,bj)=atm_tauY( i,j,bi,bj) |
TauY( i,j,bi,bj)=atm_tauY( i,j) |
71 |
Qlatent( i,j,bi,bj)=atm_qLatent( i,j,bi,bj) |
Qlatent( i,j,bi,bj)=atm_qLatent( i,j) |
72 |
Qsensible( i,j,bi,bj)=atm_qSensible( i,j,bi,bj) |
Qsensible( i,j,bi,bj)=atm_qSensible( i,j) |
73 |
Qlongwave( i,j,bi,bj)=atm_qLongwave( i,j,bi,bj) |
Qlongwave( i,j,bi,bj)=atm_qLongwave( i,j) |
74 |
Qshortwave( i,j,bi,bj)=atm_qShortwave(i,j,bi,bj) |
Qshortwave( i,j,bi,bj)=atm_qShortwave(i,j) |
75 |
uVelGround( i,j,bi,bj)=atm_uVelGround(i,j,bi,bj) |
uVelGround( i,j,bi,bj)=atm_uVelGround(i,j) |
76 |
vVelGround( i,j,bi,bj)=atm_vVelGround(i,j,bi,bj) |
vVelGround( i,j,bi,bj)=atm_vVelGround(i,j) |
77 |
FWFlux( i,j,bi,bj)=atm_FWFlux( i,j,bi,bj) |
FWFlux( i,j,bi,bj)=atm_FWFlux( i,j) |
78 |
Hatm( i,j,bi,bj)=atm_Hatm( i,j,bi,bj) |
Hatm( i,j,bi,bj)=atm_Hatm( i,j) |
79 |
ENDDO |
ENDDO |
80 |
ENDDO |
ENDDO |
81 |
ENDDO |
ENDDO |
82 |
ENDDO |
ENDDO |
83 |
|
|
|
CALL EEBOOT() |
|
|
myThid = 1 |
|
84 |
CALL INITIALISE_FIXED( myThid ) |
CALL INITIALISE_FIXED( myThid ) |
85 |
CALL INITIALISE_VARIA( myThid ) |
CALL INITIALISE_VARIA( myThid ) |
86 |
myCurrentTime = startTime |
myCurrentTime = startTime |
92 |
C Fill the export fields |
C Fill the export fields |
93 |
DO bj=myByLo(myThid),myByHi(myThid) |
DO bj=myByLo(myThid),myByHi(myThid) |
94 |
DO bi=myBxLo(myThid),myBxHi(myThid) |
DO bi=myBxLo(myThid),myBxHi(myThid) |
95 |
DO j=1-OLy,sNy+OLy |
DO j=1,sNy |
96 |
DO i=1-OLx,sNx+OLx |
DO i=1,sNx |
97 |
ocn_SSTocn(i,j,bi,bj) = SSTocn2cpl(i,j,bi,bj) |
ocn_SSTocn(i,j) = SSTocn2cpl(i,j,bi,bj) |
98 |
ocn_Hocn( i,j,bi,bj) = R_low( i,j,bi,bj) |
ocn_Hocn( i,j) = R_low( i,j,bi,bj) |
99 |
ENDDO |
ENDDO |
100 |
ENDDO |
ENDDO |
101 |
ENDDO |
ENDDO |
102 |
ENDDO |
ENDDO |
103 |
|
|
104 |
CLOSE(errorMessageUnit) |
WRITE( 6,*) 'MAXVAL SSTocn2cpl = ', MAXVAL(SSTocn2cpl(1:snx,1:sny,:,:)) |
105 |
CLOSE(standardMessageUnit) |
WRITE( 6,*) 'MINVAL SSTocn2cpl = ', MINVAL(SSTocn2cpl(1:snx,1:sny,:,:)) |
106 |
|
WRITE( 6,*) 'MAXVAL ocn_SSTocn = ', MAXVAL(ocn_SSTocn(1:snx,1:sny)) |
107 |
|
WRITE( 6,*) 'MINVAL ocn_SSTocn = ', MINVAL(ocn_SSTocn(1:snx,1:sny)) |
108 |
|
|
109 |
|
! CLOSE(errorMessageUnit) |
110 |
|
! CLOSE(standardMessageUnit) |
111 |
CALL COMP_OCN_POP_CONTEXT(1) |
CALL COMP_OCN_POP_CONTEXT(1) |
112 |
|
|
113 |
RETURN |
RETURN |