Back to home page

MITgcm

 
 

    


File indexing completed on 2018-03-02 18:42:44 UTC

view on githubraw file Latest commit add29e06 on 2018-01-31 20:35:05 UTC
42c525bfb4 Alis*0001 #include "OBCS_OPTIONS.h"
                0002 
                0003       SUBROUTINE ORLANSKI_INIT(bi, bj, myThid )
                0004 C     /==========================================================\
                0005 C     | SUBROUTINE ORLANSKI_INIT                                 |
                0006 C     | o Initialise storage and phase speed arrays              |
                0007 C     |   for Orlanski radiation condition                       |
                0008 C     |==========================================================|
                0009 C     |                                                          |
                0010 C     \==========================================================/
                0011       IMPLICIT NONE
                0012 
                0013 C     === Global variables ===
                0014 #include "SIZE.h"
                0015 #include "EEPARAMS.h"
                0016 #include "PARAMS.h"
                0017 #ifdef ALLOW_ORLANSKI
                0018 #include "ORLANSKI.h"
                0019 #endif /* ALLOW_ORLANSKI */
a24437ad80 Jean*0020 
42c525bfb4 Alis*0021 C     == Routine arguments ==
                0022       INTEGER bi, bj, myThid
                0023 
                0024 #ifdef ALLOW_ORLANSKI
                0025 
                0026 C     == Local variables ==
                0027       INTEGER I, J, K
                0028 
                0029       DO K=1,Nr
                0030         DO I=1-Olx,sNx+Olx
                0031           UN_STORE_1(I,K,bi,bj)=0.
a24437ad80 Jean*0032           VN_STORE_1(I,K,bi,bj)=0.
42c525bfb4 Alis*0033           TN_STORE_1(I,K,bi,bj)=0.
                0034           SN_STORE_1(I,K,bi,bj)=0.
                0035           UN_STORE_2(I,K,bi,bj)=0.
a24437ad80 Jean*0036           VN_STORE_2(I,K,bi,bj)=0.
42c525bfb4 Alis*0037           TN_STORE_2(I,K,bi,bj)=0.
                0038           SN_STORE_2(I,K,bi,bj)=0.
                0039           UN_STORE_3(I,K,bi,bj)=0.
a24437ad80 Jean*0040           VN_STORE_3(I,K,bi,bj)=0.
42c525bfb4 Alis*0041           TN_STORE_3(I,K,bi,bj)=0.
                0042           SN_STORE_3(I,K,bi,bj)=0.
                0043           UN_STORE_4(I,K,bi,bj)=0.
a24437ad80 Jean*0044           VN_STORE_4(I,K,bi,bj)=0.
                0045           TN_STORE_4(I,K,bi,bj)=0.
42c525bfb4 Alis*0046           SN_STORE_4(I,K,bi,bj)=0.
                0047           US_STORE_1(I,K,bi,bj)=0.
a24437ad80 Jean*0048           VS_STORE_1(I,K,bi,bj)=0.
42c525bfb4 Alis*0049           TS_STORE_1(I,K,bi,bj)=0.
                0050           SS_STORE_1(I,K,bi,bj)=0.
                0051           US_STORE_2(I,K,bi,bj)=0.
a24437ad80 Jean*0052           VS_STORE_2(I,K,bi,bj)=0.
                0053           TS_STORE_2(I,K,bi,bj)=0.
                0054           SS_STORE_2(I,K,bi,bj)=0.
42c525bfb4 Alis*0055           US_STORE_3(I,K,bi,bj)=0.
a24437ad80 Jean*0056           VS_STORE_3(I,K,bi,bj)=0.
                0057           TS_STORE_3(I,K,bi,bj)=0.
                0058           SS_STORE_3(I,K,bi,bj)=0.
42c525bfb4 Alis*0059           US_STORE_4(I,K,bi,bj)=0.
a24437ad80 Jean*0060           VS_STORE_4(I,K,bi,bj)=0.
                0061           TS_STORE_4(I,K,bi,bj)=0.
                0062           SS_STORE_4(I,K,bi,bj)=0.
42c525bfb4 Alis*0063           CVEL_UN(I,K,bi,bj)=0.
a24437ad80 Jean*0064           CVEL_VN(I,K,bi,bj)=0.
42c525bfb4 Alis*0065           CVEL_TN(I,K,bi,bj)=0.
                0066           CVEL_SN(I,K,bi,bj)=0.
                0067           CVEL_US(I,K,bi,bj)=0.
a24437ad80 Jean*0068           CVEL_VS(I,K,bi,bj)=0.
42c525bfb4 Alis*0069           CVEL_TS(I,K,bi,bj)=0.
                0070           CVEL_SS(I,K,bi,bj)=0.
                0071 #ifdef ALLOW_NONHYDROSTATIC
                0072           WN_STORE_1(I,K,bi,bj)=0.
                0073           WN_STORE_2(I,K,bi,bj)=0.
                0074           WN_STORE_3(I,K,bi,bj)=0.
                0075           WN_STORE_4(I,K,bi,bj)=0.
                0076           WS_STORE_1(I,K,bi,bj)=0.
                0077           WS_STORE_2(I,K,bi,bj)=0.
                0078           WS_STORE_3(I,K,bi,bj)=0.
                0079           WS_STORE_4(I,K,bi,bj)=0.
                0080           CVEL_WN(I,K,bi,bj)=0.
                0081           CVEL_WS(I,K,bi,bj)=0.
                0082 #endif
                0083         ENDDO
                0084         DO J=1-Oly,sNy+Oly
                0085           UE_STORE_1(J,K,bi,bj)=0.
a24437ad80 Jean*0086           VE_STORE_1(J,K,bi,bj)=0.
42c525bfb4 Alis*0087           TE_STORE_1(J,K,bi,bj)=0.
                0088           SE_STORE_1(J,K,bi,bj)=0.
                0089           UE_STORE_2(J,K,bi,bj)=0.
a24437ad80 Jean*0090           VE_STORE_2(J,K,bi,bj)=0.
                0091           TE_STORE_2(J,K,bi,bj)=0.
                0092           SE_STORE_2(J,K,bi,bj)=0.
42c525bfb4 Alis*0093           UE_STORE_3(J,K,bi,bj)=0.
a24437ad80 Jean*0094           VE_STORE_3(J,K,bi,bj)=0.
                0095           TE_STORE_3(J,K,bi,bj)=0.
                0096           SE_STORE_3(J,K,bi,bj)=0.
42c525bfb4 Alis*0097           UE_STORE_4(J,K,bi,bj)=0.
a24437ad80 Jean*0098           VE_STORE_4(J,K,bi,bj)=0.
42c525bfb4 Alis*0099           TE_STORE_4(J,K,bi,bj)=0.
                0100           SE_STORE_4(J,K,bi,bj)=0.
                0101           UW_STORE_1(J,K,bi,bj)=0.
a24437ad80 Jean*0102           VW_STORE_1(J,K,bi,bj)=0.
42c525bfb4 Alis*0103           TW_STORE_1(J,K,bi,bj)=0.
                0104           SW_STORE_1(J,K,bi,bj)=0.
                0105           UW_STORE_2(J,K,bi,bj)=0.
a24437ad80 Jean*0106           VW_STORE_2(J,K,bi,bj)=0.
                0107           TW_STORE_2(J,K,bi,bj)=0.
                0108           SW_STORE_2(J,K,bi,bj)=0.
42c525bfb4 Alis*0109           UW_STORE_3(J,K,bi,bj)=0.
a24437ad80 Jean*0110           VW_STORE_3(J,K,bi,bj)=0.
                0111           TW_STORE_3(J,K,bi,bj)=0.
                0112           SW_STORE_3(J,K,bi,bj)=0.
42c525bfb4 Alis*0113           UW_STORE_4(J,K,bi,bj)=0.
a24437ad80 Jean*0114           VW_STORE_4(J,K,bi,bj)=0.
42c525bfb4 Alis*0115           TW_STORE_4(J,K,bi,bj)=0.
                0116           SW_STORE_4(J,K,bi,bj)=0.
                0117           CVEL_UE(J,K,bi,bj)=0.
a24437ad80 Jean*0118           CVEL_VE(J,K,bi,bj)=0.
42c525bfb4 Alis*0119           CVEL_TE(J,K,bi,bj)=0.
                0120           CVEL_SE(J,K,bi,bj)=0.
                0121           CVEL_UW(J,K,bi,bj)=0.
a24437ad80 Jean*0122           CVEL_VW(J,K,bi,bj)=0.
42c525bfb4 Alis*0123           CVEL_TW(J,K,bi,bj)=0.
                0124           CVEL_SW(J,K,bi,bj)=0.
                0125 #ifdef ALLOW_NONHYDROSTATIC
                0126           WE_STORE_1(J,K,bi,bj)=0.
                0127           WE_STORE_2(J,K,bi,bj)=0.
                0128           WE_STORE_3(J,K,bi,bj)=0.
                0129           WE_STORE_4(J,K,bi,bj)=0.
                0130           WW_STORE_1(J,K,bi,bj)=0.
                0131           WW_STORE_2(J,K,bi,bj)=0.
                0132           WW_STORE_3(J,K,bi,bj)=0.
                0133           WW_STORE_4(J,K,bi,bj)=0.
                0134           CVEL_WE(J,K,bi,bj)=0.
                0135           CVEL_WW(J,K,bi,bj)=0.
                0136 #endif
                0137         ENDDO
                0138       ENDDO
                0139 
                0140 #endif /* ALLOW_ORLANSKI */
                0141       RETURN
                0142       END