/[MITgcm]/MITgcm/pkg/atm_compon_interf/cpl_exch_configs.F
ViewVC logotype

Diff of /MITgcm/pkg/atm_compon_interf/cpl_exch_configs.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.1 by jmc, Mon Dec 15 02:44:48 2003 UTC revision 1.2 by jmc, Fri May 21 20:00:48 2004 UTC
# Line 30  C     *================================= Line 30  C     *=================================
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    
# Line 58  C     o Get ocean model configuration Line 59  C     o Get ocean model configuration
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.
# Line 72  C       atmospheres full depth is 10^5 P Line 73  C       atmospheres full depth is 10^5 P
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
# Line 91  C       atmospheres full depth is 10^5 P Line 93  C       atmospheres full depth is 10^5 P
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
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22