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

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

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

revision 1.3 by dimitri, Sat Dec 28 10:11:11 2002 UTC revision 1.4 by dimitri, Tue Feb 18 05:33:54 2003 UTC
# Line 9  c Line 9  c
9  c     o  c     o
10  c  c
11  c     started: Ralf.Giering@FastOpt.de 25-Mai-2000  c     started: Ralf.Giering@FastOpt.de 25-Mai-2000
12  c     changed: menemenlis@jpl.nasa.gov 20-Dec-2002  c     mods for pkg/seaice: menemenlis@jpl.nasa.gov 20-Dec-2002
13  c  c
14  c     ==================================================================  c     ==================================================================
15  c     SUBROUTINE exf_init_runoff  c     SUBROUTINE exf_init_runoff
# Line 31  c     == routine arguments == Line 31  c     == routine arguments ==
31    
32        integer mythid        integer mythid
33    
34  #ifdef ALLOW_RUNOFF  #if (defined (ALLOW_RUNOFF) || defined (ALLOW_SEAICE))
35    
36  c     == local variables ==  c     == local variables ==
37    
38        integer bi, bj        integer bi, bj
39        integer i, j        integer i, j
40    
       _RL fac  
   
41  c     == end of interface ==  c     == end of interface ==
42    
43        do bj = mybylo(mythid), mybyhi(mythid)        do bj = mybylo(mythid), mybyhi(mythid)
# Line 62  c     == end of interface == Line 60  c     == end of interface ==
60              call exf_filter_rs( runoff, runoffmask, mythid )              call exf_filter_rs( runoff, runoffmask, mythid )
61           end if           end if
62                    
63  cds   conversion from m/yr to m/s and changing sign.  c--   scale runoff array
          fac = -1.0/(86400.*365.0)  
64           do bj = mybylo(mythid), mybyhi(mythid)           do bj = mybylo(mythid), mybyhi(mythid)
65              do bi = mybxlo(mythid), mybxhi(mythid)             do bi = mybxlo(mythid), mybxhi(mythid)
66                 do j = 1, sny               do j = 1, sny
67                    do i = 1, snx                 do i = 1, snx
68                       if (runoff(i,j,bi,bj).gt.20.) runoff(i,j,bi,bj)=0.                   if (runoff(i,j,bi,bj).gt.20.) runoff(i,j,bi,bj)=0.
69                       runoff(i,j,bi,bj) = fac*runoff(i,j,bi,bj)                   runoff(i,j,bi,bj) = exf_inscal_runoff*runoff(i,j,bi,bj)
70                    enddo                enddo
                enddo  
71              enddo              enddo
72           enddo            enddo
73            enddo
74        endif        endif
75    
76  #endif ALLOW_RUNOFF  #endif
77    
78        end        end
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22