30 |
C == Global variables == |
C == Global variables == |
31 |
#include "SIZE.h" |
#include "SIZE.h" |
32 |
#include "EEPARAMS.h" |
#include "EEPARAMS.h" |
33 |
|
#include "PARAMS.h" |
34 |
#include "GRID.h" |
#include "GRID.h" |
35 |
#include "ATMCPL.h" |
#include "ATMCPL.h" |
36 |
|
|
59 |
|
|
60 |
C Summarise fields that were imported. |
C Summarise fields that were imported. |
61 |
C o Plot ocean depths |
C o Plot ocean depths |
62 |
CALL PLOT_FIELD_XYRL( Hocn, |
CALL PLOT_FIELD_XYRL( ocMxlD, |
63 |
& 'Ocean depths on atmos grid', |
& 'Ocean mixed-layer depth on atmos grid', |
64 |
& 1, myThid ) |
& 1, myThid ) |
65 |
|
|
66 |
C Do consistency checks on imported fields. |
C Do consistency checks on imported fields. |
73 |
biDiff = 0 |
biDiff = 0 |
74 |
bjDiff = 0 |
bjDiff = 0 |
75 |
seaSurface_ocn = 0. _d 0 |
seaSurface_ocn = 0. _d 0 |
76 |
groundLevel_atm = 1. _d 5 |
c groundLevel_atm = 1. _d 5 |
77 |
|
groundLevel_atm = Ro_SeaLevel |
78 |
DO bj=myByLo(myTHid),myByHi(myThid) |
DO bj=myByLo(myTHid),myByHi(myThid) |
79 |
DO bi=myBxLo(myThid),myBxHi(myThid) |
DO bi=myBxLo(myThid),myBxHi(myThid) |
80 |
DO J=1,sNy |
DO J=1,sNy |
81 |
DO I=1,sNx |
DO I=1,sNx |
82 |
IF ( Hocn(I,J,bi,bj) .NE. seaSurface_ocn .AND. |
IF ( ocMxlD(I,J,bi,bj) .NE. seaSurface_ocn .AND. |
83 |
& Ro_surf(I,J,bi,bj) .NE. groundLevel_atm ) THEN |
& Ro_surf(I,J,bi,bj) .NE. groundLevel_atm ) THEN |
84 |
iDiff = I |
iDiff = I |
85 |
jDiff = J |
jDiff = J |
93 |
& 'H (atmosphere) ==',Ro_surf(I,J,bi,bj) |
& 'H (atmosphere) ==',Ro_surf(I,J,bi,bj) |
94 |
CALL PRINT_ERROR( msgBuf, 1 ) |
CALL PRINT_ERROR( msgBuf, 1 ) |
95 |
WRITE(msgBuf,'(A,E30.15)') |
WRITE(msgBuf,'(A,E30.15)') |
96 |
& 'H (ocean) ==',Hocn(I,J,bi,bj) |
& 'H (ocean) ==',ocMxlD(I,J,bi,bj) |
97 |
CALL PRINT_ERROR( msgBuf, 1 ) |
CALL PRINT_ERROR( msgBuf, 1 ) |
98 |
ENDIF |
ENDIF |
99 |
ENDDO |
ENDDO |