/[MITgcm]/MITgcm/pkg/exf/exf_getffields.F
ViewVC logotype

Diff of /MITgcm/pkg/exf/exf_getffields.F

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

revision 1.58 by gforget, Sun Nov 2 21:23:30 2014 UTC revision 1.59 by gforget, Fri May 22 23:13:23 2015 UTC
# Line 391  C     Runoff temperature Line 391  C     Runoff temperature
391       &     myTime, myIter, myThid )       &     myTime, myIter, myThid )
392  #endif /* ALLOW_RUNOFTEMP */  #endif /* ALLOW_RUNOFTEMP */
393    
394    #ifdef ALLOW_ROTATE_UV_CONTROLS
395          if ( useCTRL ) then
396            DO bj = myByLo(myThid),myByHi(myThid)
397             DO bi = myBxLo(myThid),mybxhi(myThid)
398              DO j = 1-OLy,sNy+OLy
399               DO i = 1-OLx,sNx+OLx
400                 tmpUE(i,j,bi,bj) = 0. _d 0
401                 tmpVN(i,j,bi,bj) = 0. _d 0
402                 tmpUX(i,j,bi,bj) = 0. _d 0
403                 tmpVY(i,j,bi,bj) = 0. _d 0
404               ENDDO
405              ENDDO
406             ENDDO
407            ENDDO
408          endif
409    #endif
410    
411  # if (!defined (ALLOW_ECCO) || defined (ECCO_CTRL_DEPRECATED))  # if (!defined (ALLOW_ECCO) || defined (ECCO_CTRL_DEPRECATED))
412    
413  C-- Control variables for atmos. state  C-- Control variables for atmos. state
# Line 505  C-- Control variables for atmos. state Line 522  C-- Control variables for atmos. state
522  #else  #else
523    
524  #if defined(ALLOW_UWIND_CONTROL) && defined(ALLOW_VWIND_CONTROL)  #if defined(ALLOW_UWIND_CONTROL) && defined(ALLOW_VWIND_CONTROL)
         DO bj = myByLo(myThid),myByHi(myThid)  
          DO bi = myBxLo(myThid),mybxhi(myThid)  
           DO j = 1-OLy,sNy+OLy  
            DO i = 1-OLx,sNx+OLx  
              tmpUE(i,j,bi,bj) = 0. _d 0  
              tmpVN(i,j,bi,bj) = 0. _d 0  
              tmpUX(i,j,bi,bj) = 0. _d 0  
              tmpVY(i,j,bi,bj) = 0. _d 0  
            ENDDO  
           ENDDO  
          ENDDO  
         ENDDO  
   
525        CALL CTRL_GET_GEN (        CALL CTRL_GET_GEN (
526       &     xx_uwind_file, xx_uwindstartdate, xx_uwindperiod,       &     xx_uwind_file, xx_uwindstartdate, xx_uwindperiod,
527       &     maskc, tmpUE, xx_uwind0, xx_uwind1, xx_uwind_dummy,       &     maskc, tmpUE, xx_uwind0, xx_uwind1, xx_uwind_dummy,
# Line 613  cdm #endif Line 617  cdm #endif
617             if (xx_gentim2d_file(iarr).EQ.'xx_precip')             if (xx_gentim2d_file(iarr).EQ.'xx_precip')
618       &       precip(i,j,bi,bj)=precip(i,j,bi,bj)+       &       precip(i,j,bi,bj)=precip(i,j,bi,bj)+
619       &                         xx_gentim2d(i,j,bi,bj,iarr)       &                         xx_gentim2d(i,j,bi,bj,iarr)
620    #ifndef ALLOW_ROTATE_UV_CONTROLS
621               if (xx_gentim2d_file(iarr).EQ.'xx_uwind')
622         &       uwind(i,j,bi,bj)=uwind(i,j,bi,bj)+
623         &                         xx_gentim2d(i,j,bi,bj,iarr)
624               if (xx_gentim2d_file(iarr).EQ.'xx_vwind')
625         &       vwind(i,j,bi,bj)=vwind(i,j,bi,bj)+
626         &                         xx_gentim2d(i,j,bi,bj,iarr)
627    #else
628               if (xx_gentim2d_file(iarr).EQ.'xx_uwind')
629         &       uwind(i,j,bi,bj)=uwind(i,j,bi,bj)+
630         &                         xx_gentim2d(i,j,bi,bj,iarr)
631               if (xx_gentim2d_file(iarr).EQ.'xx_vwind')
632         &       vwind(i,j,bi,bj)=vwind(i,j,bi,bj)+
633         &                         xx_gentim2d(i,j,bi,bj,iarr)
634    
635    #endif
636             enddo             enddo
637           ENDDO           ENDDO
638          ENDDO          ENDDO
639         ENDDO         ENDDO
640        ENDDO        ENDDO
641    #ifdef ALLOW_ROTATE_UV_CONTROLS
642          CALL ROTATE_UV2EN_RL(tmpUX,tmpVY,tmpUE,tmpVN,
643         &     .FALSE.,.FALSE.,.TRUE.,1,myThid)
644    
645           DO bj = myByLo(myThid),myByHi(myThid)
646             DO bi = myBxLo(myThid),mybxhi(myThid)
647              DO j = 1,sNy
648               DO i = 1,sNx
649                 uwind(i,j,bi,bj)=uwind(i,j,bi,bj)+tmpUX(i,j,bi,bj)
650                 vwind(i,j,bi,bj)=vwind(i,j,bi,bj)+tmpVY(i,j,bi,bj)
651               ENDDO
652              ENDDO
653             ENDDO
654           ENDDO
655    #endif /* ALLOW_ROTATE_UV_CONTROLS */
656    
657        endif !if (ctrlUseGen) then        endif !if (ctrlUseGen) then
658  #endif  #endif
659    
Legend:
Removed from v.1.58  
changed lines
  Added in v.1.59
  ViewVC Help
Powered by ViewVC 1.1.22