/[MITgcm]/MITgcm_contrib/jscott/pkg_atm2d/get_ocnvars.F
ViewVC logotype

Diff of /MITgcm_contrib/jscott/pkg_atm2d/get_ocnvars.F

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

revision 1.1 by jscott, Fri Aug 11 18:55:50 2006 UTC revision 1.3 by jscott, Tue May 1 19:38:47 2007 UTC
# Line 4  Line 4 
4  C     !INTERFACE:  C     !INTERFACE:
5        SUBROUTINE GET_OCNVARS( myTime, myIter, myThid )        SUBROUTINE GET_OCNVARS( myTime, myIter, myThid )
6  C     *==========================================================*  C     *==========================================================*
7  C     | |  C     | Gets the ocean data from either the ocean common block,  |
8  c     | |  C     | or, in MPI mode, from the passed common block. This data |
9    C     | is used by the atmosphere and sea ice models.            |
10  C     *==========================================================*  C     *==========================================================*
11          IMPLICIT NONE          IMPLICIT NONE
12    
# Line 16  C     === Global Atmosphere Variables == Line 17  C     === Global Atmosphere Variables ==
17  #include "THSICE_VARS.h"  #include "THSICE_VARS.h"
18  #include "ATM2D_VARS.h"  #include "ATM2D_VARS.h"
19    
20  #ifdef MPI_ON  #ifdef ATM2D_MPI_ON
21  #  include "OCNSIZE.h"  #  include "OCNSIZE.h"
22  #  include "OCNVARS.h"  #  include "OCNVARS.h"
23  #else  #else
# Line 45  C presuming ocnsize.h same as size.h Line 46  C presuming ocnsize.h same as size.h
46          DO i=1,sNx          DO i=1,sNx
47            sstFromOcn(i,j)= SST_ocn(i,j)            sstFromOcn(i,j)= SST_ocn(i,j)
48            sssFromOcn(i,j)= SSS_ocn(i,j)            sssFromOcn(i,j)= SSS_ocn(i,j)
49            mlDepth(i,j) = OcMx1D_ocn(i,j)            mlDepth(i,j) = OcMxlD_ocn(i,j)
50  C         fluxCO2(i,j) = new variable passed  C         fluxCO2(i,j) = new variable passed
51    
52  C Put in seaice common block  C Put in seaice common block
53            tOceMxl(i,j,1,1)= sst(i,j)            tOceMxl(i,j,1,1)= SST_ocn(i,j)
54            sOceMxl(i,j,1,1)= sss(i,j)            sOceMxl(i,j,1,1)= SSS_ocn(i,j)
55            hOceMxl(i,j,1,1) = mlDepth(i,j)            hOceMxl(i,j,1,1) = OcMxlD_ocn(i,j)
56            v2ocMxL(i,j,1,1) =  vSq_ocn(i,j)            v2ocMxL(i,j,1,1) =  vSq_ocn(i,j)
57          ENDDO          ENDDO
58        ENDDO        ENDDO
59    
60  #else  #else
 CJRS ask J-M about this next variable  
 C     IF ( cpl_earlyExpImpCall ) THEN  
61         DO j=1,sNy         DO j=1,sNy
62          DO i=1,sNx          DO i=1,sNx
63            sstFromOcn(i,j) = theta(i,j,1,1,1)            sstFromOcn(i,j) = theta(i,j,1,1,1)
64            sssFromOcn(i,j) = salt (i,j,1,1,1)            sssFromOcn(i,j) = salt (i,j,1,1,1)
65          ENDDO          ENDDO
66         ENDDO         ENDDO
67  C      ELSE  
 C       DO j=1,sNy  
 C        DO i=1,sNx  
 C          sstFromOcn(i,j) = gT(i,j,1,1,1)  
 C          sssFromOcn(i,j) = gS(i,j,1,1,1)  
 C        ENDDO  
 C       ENDDO  
 C      ENDIF  
68  #  ifdef NONLIN_FRSURF  #  ifdef NONLIN_FRSURF
69  Cjrs ask J-M about this  Cjrs ask J-M about next lines
70        IF ( nonlinFreeSurf.GT.0 .AND.        IF ( nonlinFreeSurf.GT.0 .AND.
71       &    (staggerTimeStep .OR. .NOT.cpl_earlyExpImpCall ) ) THEN       &    (staggerTimeStep .OR. .NOT.cpl_earlyExpImpCall ) ) THEN
72         IF ( select_rStar.GT.0 ) THEN         IF ( select_rStar.GT.0 ) THEN
# Line 118  Cjrs ask J-M about this Line 110  Cjrs ask J-M about this
110          DO i=1,sNx          DO i=1,sNx
111    
112  C         fluxCO2(i,j) = new variable passed  C         fluxCO2(i,j) = new variable passed
113    
114  C Put in seaice common block  C Put in seaice common block
115            tOceMxl(i,j,1,1)= sstFromOcn(i,j)            tOceMxl(i,j,1,1)= sstFromOcn(i,j)
116            sOceMxl(i,j,1,1)= sssFromOcn(i,j)            sOceMxl(i,j,1,1)= sssFromOcn(i,j)
# Line 128  C         v2ocMxL(i,j,1,1) =  velSq(i,j) Line 121  C         v2ocMxL(i,j,1,1) =  velSq(i,j)
121                    
122  #endif  #endif
123    
124        PRINT *,'After ocean step, sst,sss:', theta(JBUGI,JBUGJ,1,1,1),  C      PRINT *,'After ocean step, sst,sss:', theta(JBUGI,JBUGJ,1,1,1),
125       &         salt(JBUGI,JBUGJ,1,1,1)  C     &         salt(JBUGI,JBUGJ,1,1,1)
126        PRINT *,'***get_ocnvars, before ice thicken/extend'  C      PRINT *,'***get_ocnvars, before ice thicken/extend'
127        PRINT *,'Icemask: ',icemask(JBUGI,JBUGJ,1,1)  C      PRINT *,'Icemask: ',icemask(JBUGI,JBUGJ,1,1)
128        PRINT *,'Iceheight: ',iceheight(JBUGI,JBUGJ,1,1)  C      PRINT *,'Iceheight: ',iceheight(JBUGI,JBUGJ,1,1)
129        PRINT *,'pass_runoff:',pass_runoff(JBUGI,JBUGJ)  C      PRINT *,'pass_runoff:',pass_runoff(JBUGI,JBUGJ)
130        PRINT *,'pass_precip:',pass_precip(JBUGI,JBUGJ)  C      PRINT *,'pass_precip:',pass_precip(JBUGI,JBUGJ)
131        PRINT *,'pass_evap:', pass_evap(JBUGI,JBUGJ)  C      PRINT *,'pass_evap:', pass_evap(JBUGI,JBUGJ)
132        PRINT *,'pass_qnet:',pass_qnet(JBUGI,JBUGJ)  C      PRINT *,'pass_qnet:',pass_qnet(JBUGI,JBUGJ)
133        PRINT *,'icFrwAtm',icFrwAtm(JBUGI,JBUGJ,1,1)  C      PRINT *,'icFrwAtm',icFrwAtm(JBUGI,JBUGJ,1,1)
134        PRINT *,'icFlxSw:',icFlxSw(JBUGI,JBUGJ,1,1)  C      PRINT *,'icFlxSw:',icFlxSw(JBUGI,JBUGJ,1,1)
135        PRINT *,'pass_Prcatm',pass_prcAtm(JBUGI,JBUGJ)  C      PRINT *,'pass_Prcatm',pass_prcAtm(JBUGI,JBUGJ)
136        PRINT *,tOceMxl(JBUGI,JBUGJ,1,1),v2OcMxl(JBUGI,JBUGJ,1,1)  C      PRINT *,tOceMxl(JBUGI,JBUGJ,1,1),v2OcMxl(JBUGI,JBUGJ,1,1)
137    
138        RETURN        RETURN
139        END        END
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22