/[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.2 by jscott, Tue Aug 22 20:21:39 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 57  C Put in seaice common block Line 58  C Put in seaice common block
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 130  C         v2ocMxL(i,j,1,1) =  velSq(i,j) Line 123  C         v2ocMxL(i,j,1,1) =  velSq(i,j)
123    
124  C      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  C     &         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.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22