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

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

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

revision 1.2 by mlosch, Tue May 30 22:47:08 2006 UTC revision 1.3 by mlosch, Sat Jun 3 21:36:15 2006 UTC
# Line 38  c     == routine arguments == Line 38  c     == routine arguments ==
38        integer myiter        integer myiter
39        _RL     mytime        _RL     mytime
40    
 #ifdef ALLOW_BULKFORMULAE  
   
41  c     == local variables ==  c     == local variables ==
42    
43        integer bi,bj        integer bi,bj
# Line 53  c     == external functions == Line 51  c     == external functions ==
51        integer  ilnblnk        integer  ilnblnk
52        external ilnblnk        external ilnblnk
53    
54  #ifndef ALLOW_ATM_WIND        _RL       tmp1
55        _RL       TMP1        _RL       tmp2
56        _RL       TMP2        _RL       tmp3
57        _RL       TMP3        _RL       tmp4
58        _RL       TMP4        _RL       tmp5
       _RL       TMP5  
 #endif  
59    
60  c     == end of interface ==  c     == end of interface ==
61    
# Line 93  c             Wind speed and direction. Line 89  c             Wind speed and direction.
89            endif            endif
90  #else  /* ifndef ALLOW_ATM_WIND */  #else  /* ifndef ALLOW_ATM_WIND */
91  c  c
 #ifdef ALLOW_ATM_TEMP  
 c  
92  c             The variables us, sh and rdn have to be computed from  c             The variables us, sh and rdn have to be computed from
93  c             given wind stresses inverting relationship for neutral  c             given wind stresses inverting relationship for neutral
94  c             drag coeff. cdn.  c             drag coeff. cdn.
# Line 127  c             cdn(umps); ustar can be di Line 121  c             cdn(umps); ustar can be di
121               us(i,j,bi,bj)   = (tmp4 + tmp5)**(1/3) +               us(i,j,bi,bj)   = (tmp4 + tmp5)**(1/3) +
122       &            tmp3**2 * (tmp4 + tmp5)**(-1/3) - tmp3       &            tmp3**2 * (tmp4 + tmp5)**(-1/3) - tmp3
123            endif            endif
124  c                    uwind(i,j,bi,bj) = us(i,j,bi,bj)*cw(i,j,bi,bj)
125  #endif /* ALLOW_ATM_TEMP */            vwind(i,j,bi,bj) = us(i,j,bi,bj)*sw(i,j,bi,bj)
126  c  c
127  #endif /* ifndef ALLOW_ATM_WIND */  #endif /* ifndef ALLOW_ATM_WIND */
128    
# Line 149  c--   otherwise fill wspeed array for la Line 143  c--   otherwise fill wspeed array for la
143         enddo         enddo
144        enddo        enddo
145    
146  #endif /* ALLOW_BULKFORMULAE */        return
   
147        end        end
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22