97 |
ENDIF |
ENDIF |
98 |
IF ( surfQFile .NE. ' ' ) THEN |
IF ( surfQFile .NE. ' ' ) THEN |
99 |
CALL READ_FLD_XY_RS( surfQFile, ' ', Qnet, 0, myThid ) |
CALL READ_FLD_XY_RS( surfQFile, ' ', Qnet, 0, myThid ) |
100 |
|
ELSEIF ( surfQnetFile .NE. ' ' ) THEN |
101 |
|
CALL READ_FLD_XY_RS( surfQnetFile, ' ', Qnet, 0, myThid ) |
102 |
ENDIF |
ENDIF |
103 |
IF ( EmPmRfile .NE. ' ' ) THEN |
IF ( EmPmRfile .NE. ' ' ) THEN |
104 |
CALL READ_FLD_XY_RS( EmPmRfile, ' ', EmPmR, 0, myThid ) |
CALL READ_FLD_XY_RS( EmPmRfile, ' ', EmPmR, 0, myThid ) |
112 |
#ifdef SHORTWAVE_HEATING |
#ifdef SHORTWAVE_HEATING |
113 |
IF ( surfQswFile .NE. ' ' ) THEN |
IF ( surfQswFile .NE. ' ' ) THEN |
114 |
CALL READ_FLD_XY_RS( surfQswFile, ' ', Qsw, 0, myThid ) |
CALL READ_FLD_XY_RS( surfQswFile, ' ', Qsw, 0, myThid ) |
115 |
|
IF ( surfQFile .NE. ' ' ) THEN |
116 |
|
C- Qnet is now (after c54) the net Heat Flux (including SW) |
117 |
|
DO j=1-OLy,sNy+OLy |
118 |
|
DO i=1-OLx,sNx+OLx |
119 |
|
Qnet(i,j,bi,bj) = Qnet(i,j,bi,bj) + Qsw(i,j,bi,bj) |
120 |
|
ENDDO |
121 |
|
ENDDO |
122 |
|
ENDIF |
123 |
ENDIF |
ENDIF |
124 |
#endif |
#endif |
125 |
#ifdef ATMOSPHERIC_LOADING |
#ifdef ATMOSPHERIC_LOADING |