/[MITgcm]/MITgcm/pkg/flt/flt_runga2.F
ViewVC logotype

Diff of /MITgcm/pkg/flt/flt_runga2.F

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

revision 1.1 by adcroft, Thu Sep 13 17:43:56 2001 UTC revision 1.2 by molod, Mon Jul 19 15:13:07 2004 UTC
# Line 51  c     == local variables == Line 51  c     == local variables ==
51        character*(max_len_mbuf) msgbuf        character*(max_len_mbuf) msgbuf
52        _RL npart_dist        _RL npart_dist
53        Real*8 PORT_RAND        Real*8 PORT_RAND
54    #undef _USE_INTEGERS
55    #ifdef _USE_INTEGERS
56          integer seed
57    #else
58          Real*8 seed
59    #endif
60    
61  c     == end of interface ==  c     == end of interface ==
62    
# Line 107  c Line 113  c
113  #endif  #endif
114    
115              if (iup(ip,bi,bj).ne.-2.) then              if (iup(ip,bi,bj).ne.-2.) then
116                 uu = uu + uu*(PORT_RAND()-0.5)*flt_noise                 uu = uu + uu*(PORT_RAND(seed)-0.5)*flt_noise
117                 vv = vv + vv*(PORT_RAND()-0.5)*flt_noise                 vv = vv + vv*(PORT_RAND(seed)-0.5)*flt_noise
118              endif              endif
119    
120  c xx and xt are in indices. Therefore it is necessary to multiply  c xx and xt are in indices. Therefore it is necessary to multiply
# Line 135  c Line 141  c
141  #endif  #endif
142    
143              if (iup(ip,bi,bj).ne.-2.) then              if (iup(ip,bi,bj).ne.-2.) then
144                 u1 = u1 + u1*(PORT_RAND()-0.5)*flt_noise                 u1 = u1 + u1*(PORT_RAND(seed)-0.5)*flt_noise
145                 v1 = v1 + v1*(PORT_RAND()-0.5)*flt_noise                 v1 = v1 + v1*(PORT_RAND(seed)-0.5)*flt_noise
146              endif              endif
147    
148  c xpart is in coordinates. Therefore it is necessary to multiply  c xpart is in coordinates. Therefore it is necessary to multiply
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22