pkg/atm_phys/atm_phys_init_varia.F

b2ea1d2979 Jean*0001 #include "ATM_PHYS_OPTIONS.h"
                0002
                0003 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
                0004 CBOP 0
                0005 C !ROUTINE: ATM_PHYS_INIT_VARIA
                0006
                0007 C !INTERFACE:
                0008       SUBROUTINE ATM_PHYS_INIT_VARIA( myThid )
                0009
                0010 C     !DESCRIPTION:
                0011 C     Initialize Atm_Phys variables
                0012
                0013 C     !USES:
                0014       IMPLICIT NONE
                0015 #include "SIZE.h"
                0016 #include "EEPARAMS.h"
                0017 #include "PARAMS.h"
                0018 #include "GRID.h"
                0019 #include "ATM_PHYS_PARAMS.h"
                0020 #include "ATM_PHYS_VARS.h"
                0021
                0022 C     !INPUT/OUTPUT PARAMETERS:
                0023 C     myThid ::  my Thread Id number
                0024       INTEGER myThid
                0025 CEOP
                0026
                0027 #ifdef ALLOW_ATM_PHYS
                0028
                0029 C     !LOCAL VARIABLES:
                0030 C     === Local variables ===
                0031       INTEGER i,j,k,bi,bj
                0032
                0033       DO bj = myByLo(myThid), myByHi(myThid)
                0034        DO bi = myBxLo(myThid), myBxHi(myThid)
                0035
8a8a10bf89 Jean*0036 C--   Initialize arrays in common blocks :
b2ea1d2979 Jean*0037 C     3-dim.
                0038         DO k=1,Nr
                0039          DO j=1-OLy,sNy+OLy
                0040           DO i=1-OLx,sNx+OLx
                0041            atmPhys_dT(i,j,k,bi,bj) = 0. _d 0
                0042            atmPhys_dQ(i,j,k,bi,bj) = 0. _d 0
                0043            atmPhys_dU(i,j,k,bi,bj) = 0. _d 0
                0044            atmPhys_dV(i,j,k,bi,bj) = 0. _d 0
                0045           ENDDO
                0046          ENDDO
                0047         ENDDO
                0048
                0049 C     2-dim.
                0050         DO j=1-OLy,sNy+OLy
                0051          DO i=1-OLx,sNx+OLx
                0052           atmPhys_SST(i,j,bi,bj) = 0. _d 0
                0053          ENDDO
                0054         ENDDO
                0055
                0056 C- end bi,bj loops
                0057        ENDDO
                0058       ENDDO
                0059
8a8a10bf89 Jean*0060 C--   Initialize Boundary conditions:
                0061       IF ( atmPhys_SSTFile .NE. ' ' ) THEN
b2ea1d2979 Jean*0062          CALL READ_FLD_XY_RL( atmPhys_SSTFile,' ',
                0063      &                        atmPhys_SST, 0, myThid )
                0064          CALL EXCH_XY_RL( atmPhys_SST, myThid )
8a8a10bf89 Jean*0065       ENDIF
                0066
                0067 C--   Initialize ATM_PHYS State Variables:
                0068       IF ( startTime.EQ.baseTime .AND. nIter0.EQ.0
                0069      &                           .AND. pickupSuff.EQ.' ' ) THEN
                0070 C-    start new run from iter-zero:
b2ea1d2979 Jean*0071
                0072       ELSE
8a8a10bf89 Jean*0073 C-    restart from a pickup:
b2ea1d2979 Jean*0074          CALL ATM_PHYS_READ_PICKUP( nIter0, myThid )
                0075
8a8a10bf89 Jean*0076 C-    end start-from-iter-zero if/else block
b2ea1d2979 Jean*0077       ENDIF
                0078
                0079 #endif /* ALLOW_ATM_PHYS */
                0080
                0081       RETURN
                0082       END