Back to home page

MITgcm

 
 

    


File indexing completed on 2018-03-02 18:39:17 UTC

view on githubraw file Latest commit add29e06 on 2018-01-31 20:35:05 UTC
d8206d87ee Patr*0001 #include "EBM_OPTIONS.h"
                0002 
d2780bc091 Jean*0003 CBOP 0
                0004 C !ROUTINE: EBM_INI_VARS
                0005 
                0006 C !INTERFACE:
d8206d87ee Patr*0007       SUBROUTINE EBM_INI_VARS( myThid )
                0008 
d2780bc091 Jean*0009 C     !DESCRIPTION:
                0010 C     *==========================================================*
                0011 C     | S/R INITIALISE VARIABLES FOR EBM
                0012 C     *==========================================================*
                0013 
                0014 C     !USES:
d8206d87ee Patr*0015       IMPLICIT NONE
                0016 
                0017 C     === Global variables ===
                0018 #include "SIZE.h"
                0019 #include "EEPARAMS.h"
                0020 #include "PARAMS.h"
d2780bc091 Jean*0021 c#include "DYNVARS.h"
                0022 c#include "GRID.h"
                0023 c#include "FFIELDS.h"
d8206d87ee Patr*0024 #ifdef ALLOW_EBM
                0025 # include "EBM.h"
                0026 #endif
                0027 
d2780bc091 Jean*0028 C     !INPUT PARAMETERS:
                0029 C     myThid  :: my Thread Id number
d8206d87ee Patr*0030       INTEGER myThid
d2780bc091 Jean*0031 CEOP
d8206d87ee Patr*0032 
                0033 #ifdef ALLOW_EBM
d2780bc091 Jean*0034 C     !LOCAL VARIABLES:
                0035 C     == Local variables ==
                0036 C     i, j, k :: Loop counters
                0037       INTEGER bi,bj,i,j
d8206d87ee Patr*0038 
                0039 C--   initially all forcing fields are zero
                0040       DO bj = myByLo(myThid), myByHi(myThid)
                0041        DO bi = myBxLo(myThid), myBxHi(myThid)
                0042         DO j=1-OLy,sNy+OLy
                0043          DO i=1-OLx,sNx+OLx
d2780bc091 Jean*0044           Run    (i,j,bi,bj) = 0.
d8206d87ee Patr*0045           winPert(i,j,bi,bj) = 0.
                0046          ENDDO
                0047         ENDDO
                0048        ENDDO
                0049       ENDDO
                0050 
                0051       CALL EBM_LOAD_CLIMATOLOGY( myThid )
                0052 
d2780bc091 Jean*0053 c     CALL EBM_DRIVER( startTime, niter0, myThid )
d8206d87ee Patr*0054 
d2780bc091 Jean*0055 c     print *, TmlS - 273.15, TtS - 273.15
                0056 c     print *, TmlN - 273.15, TtN - 273.15
d8206d87ee Patr*0057 
                0058 #endif /* ALLOW_EBM */
                0059 
                0060       RETURN
                0061       END