/[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.48 by heimbach, Tue Jul 31 16:08:16 2012 UTC revision 1.49 by gforget, Tue Aug 28 19:17:46 2012 UTC
# Line 60  c     == end of interface == Line 60  c     == end of interface ==
60  c--   read forcing fields from files and temporal interpolation  c--   read forcing fields from files and temporal interpolation
61    
62  c     Zonal and meridional wind stress.  c     Zonal and meridional wind stress.
63          IF ( .NOT.useAtmWind ) THEN
64        call exf_set_uv(        call exf_set_uv(
65       I     ustressfile, ustressstartdate, ustressperiod,       I     ustressfile, ustressstartdate, ustressperiod,
66       I     exf_inscal_ustress,       I     exf_inscal_ustress,
# Line 77  c     Zonal and meridional wind stress. Line 78  c     Zonal and meridional wind stress.
78       I     uvInterp_stress,       I     uvInterp_stress,
79  #endif /* USE_EXF_INTERPOLATION */  #endif /* USE_EXF_INTERPOLATION */
80       I     mytime, myiter, mythid )       I     mytime, myiter, mythid )
81          ELSE
82           do bj = mybylo(mythid),mybyhi(mythid)
83            do bi = mybxlo(mythid),mybxhi(mythid)
84             do j=1-Oly,sNy+Oly
85              do i=1-Olx,sNx+Olx
86               ustress(i,j,bi,bj) = 0. _d 0
87               vstress(i,j,bi,bj) = 0. _d 0
88              enddo
89             enddo
90            enddo
91           enddo
92          ENDIF
93    
94  c--   wind speed  c--   wind speed
95        call exf_set_gen(        call exf_set_gen(
# Line 91  c--   wind speed Line 104  c--   wind speed
104  #endif  #endif
105       &     mytime, myiter, mythid )       &     mytime, myiter, mythid )
106    
 #ifdef ALLOW_ATM_WIND  
   
107  c     Zonal and meridional wind.  c     Zonal and meridional wind.
108          IF ( useAtmWind ) THEN
109        call exf_set_uv(        call exf_set_uv(
110       I     uwindfile, uwindstartdate, uwindperiod,       I     uwindfile, uwindstartdate, uwindperiod,
111       I     exf_inscal_uwind,       I     exf_inscal_uwind,
# Line 127  C     Subtract UVEL and VVEL from UWIND Line 139  C     Subtract UVEL and VVEL from UWIND
139         enddo         enddo
140        endif        endif
141    
142  #endif /* ALLOW_ATM_WIND */        ELSE
143           do bj = mybylo(mythid),mybyhi(mythid)
144            do bi = mybxlo(mythid),mybxhi(mythid)
145             do j=1-Oly,sNy+Oly
146              do i=1-Olx,sNx+Olx
147               uwind(i,j,bi,bj) = 0. _d 0
148               vwind(i,j,bi,bj) = 0. _d 0
149              enddo
150             enddo
151            enddo
152           enddo
153          ENDIF
154    
155  c     Atmospheric heat flux.  c     Atmospheric heat flux.
156        call exf_set_gen  (        call exf_set_gen  (
# Line 427  c-- Control variables for atmos. state Line 450  c-- Control variables for atmos. state
450       &     mytime, myiter, mythid )       &     mytime, myiter, mythid )
451  #endif  #endif
452    
453          IF ( useAtmWind ) THEN
454  #ifndef ALLOW_ROTATE_UV_CONTROLS  #ifndef ALLOW_ROTATE_UV_CONTROLS
455    
456  #ifdef ALLOW_UWIND_CONTROL  #ifdef ALLOW_UWIND_CONTROL
# Line 489  c-- Control variables for atmos. state Line 513  c-- Control variables for atmos. state
513  #endif  #endif
514    
515  #endif /* ALLOW_ROTATE_UV_CONTROLS */  #endif /* ALLOW_ROTATE_UV_CONTROLS */
516          ENDIF
517    
518  #ifdef ALLOW_ATM_MEAN_CONTROL  #ifdef ALLOW_ATM_MEAN_CONTROL
519        do bj = mybylo(mythid),mybyhi(mythid)        do bj = mybylo(mythid),mybyhi(mythid)
Legend:
Removed from v.1.48  
changed lines
  Added in v.1.49
  ViewVC Help
Powered by ViewVC 1.1.22