/[MITgcm]/MITgcm/pkg/aim_v23/aim_initialise.F
ViewVC logotype

Diff of /MITgcm/pkg/aim_v23/aim_initialise.F

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

revision 1.7 by jmc, Thu Jun 24 23:43:11 2004 UTC revision 1.8 by jmc, Sun Nov 14 19:54:01 2004 UTC
# Line 10  C     *================================= Line 10  C     *=================================
10  C     | Initialisation of AIM atmospheric physics package :  C     | Initialisation of AIM atmospheric physics package :
11  C     | 1) call iniphys (=> set parameters to default value)  C     | 1) call iniphys (=> set parameters to default value)
12  C     | 2) read AIM parameters  C     | 2) read AIM parameters
 C     | 3) initialisation of AIM arrays (time-ave)  
13  C     *==================================================================*  C     *==================================================================*
14        IMPLICIT NONE        IMPLICIT NONE
15    
# Line 113  C      of surface reference pressure Ro_ Line 112  C      of surface reference pressure Ro_
112         ENDDO         ENDDO
113        ENDDO        ENDDO
114                
115  C--   Initialise surface forcing fields (in AIM_FFIELDS.h):  C--   Initialise Land Fraction (in AIM_FFIELDS.h):
 c     DO bj = myByLo(myThid), myByHi(myThid)  
 c      DO bi = myBxLo(myThid), myBxHi(myThid)  
 c       DO j=1-Oly,sNy+Oly  
 c        DO i=1-Olx,sNx+Olx    
 c         aim_surfTemp(i,j,bi,bj) = 300.  
 c         aim_soilWater(i,j,bi,bj) = 0.  
 c         aim_albedo   (i,j,bi,bj) = 0.  
 c         aim_landFr   (i,j,bi,bj) = 0.  
 c        ENDDO  
 c       ENDDO  
 c      ENDDO  
 c     ENDDO  
   
 #ifdef ALLOW_AIM_TAVE  
 C     Initialise diagnostic counters (these are cleared on model start  
 C      i.e. not loaded from history file for now ).  
116        DO bj = myByLo(myThid), myByHi(myThid)        DO bj = myByLo(myThid), myByHi(myThid)
117         DO bi = myBxLo(myThid), myBxHi(myThid)         DO bi = myBxLo(myThid), myBxHi(myThid)
118          CALL TIMEAVE_RESET(USTRtave, 1, bi, bj, myThid)          DO j=1-Oly,sNy+Oly
119          CALL TIMEAVE_RESET(VSTRtave, 1, bi, bj, myThid)           DO i=1-Olx,sNx+Olx  
120          CALL TIMEAVE_RESET(TSRtave, 1, bi, bj, myThid)            aim_landFr   (i,j,bi,bj) = 0.
121          CALL TIMEAVE_RESET(OLRtave, 1, bi, bj, myThid)           ENDDO
         CALL TIMEAVE_RESET(SSRtave, 1, bi, bj, myThid)  
         CALL TIMEAVE_RESET(SLRtave, 1, bi, bj, myThid)  
         CALL TIMEAVE_RESET(SHFtave, 1, bi, bj, myThid)  
         CALL TIMEAVE_RESET(EVAPtave, 1, bi, bj, myThid)  
         CALL TIMEAVE_RESET(PRECNVtave,1, bi, bj, myThid)  
         CALL TIMEAVE_RESET(PRECLStave,1, bi, bj, myThid)  
         CALL TIMEAVE_RESET(CLOUDCtave,1, bi, bj, myThid)  
         CALL TIMEAVE_RESET(CLTOPtave, 1, bi, bj, myThid)  
         CALL TIMEAVE_RESET(CBMFtave,  1, bi, bj, myThid)  
         CALL TIMEAVE_RESET(DRAGtave,  1, bi, bj, myThid)  
         CALL TIMEAVE_RESET(aimV0tave, 1, bi, bj, myThid)  
         CALL TIMEAVE_RESET(aimT0tave, 1, bi, bj, myThid)  
         CALL TIMEAVE_RESET(aimQ0tave, 1, bi, bj, myThid)  
         CALL TIMEAVE_RESET(EnFxPrtave,1, bi, bj, myThid)  
         CALL TIMEAVE_RESET(albedotave,1, bi, bj, myThid)  
         CALL TIMEAVE_RESET(dTsurftave,1, bi, bj, myThid)  
         CALL TIMEAVE_RESET(aimRHtave,Nr, bi, bj, myThid)  
         DO k=1,Nr  
           aim_timeAve(k,bi,bj) = 0.  
122          ENDDO          ENDDO
123         ENDDO         ENDDO
124        ENDDO        ENDDO
125  #endif /* ALLOW_AIM_TAVE */  
126          IF ( aim_LandFile .NE. ' '  ) THEN
127             CALL READ_REC_XY_RS(aim_LandFile,aim_landFr,1,nIter0,myThid)
128          ENDIF
129    
130  #endif /* ALLOW_AIM */  #endif /* ALLOW_AIM */
131    
Legend:
Removed from v.1.7  
changed lines
  Added in v.1.8
  ViewVC Help
Powered by ViewVC 1.1.22