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

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

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

Revision 1.4 - (show annotations) (download)
Tue Feb 18 05:33:54 2003 UTC (21 years, 4 months ago) by dimitri
Branch: MAIN
CVS Tags: checkpoint48f_post, checkpoint48i_post, checkpoint50e_post, checkpoint50c_post, checkpoint48h_post, checkpoint50c_pre, checkpoint50h_post, checkpoint50d_pre, checkpoint50b_pre, checkpoint51e_post, checkpoint51b_post, checkpoint51c_post, checkpoint49, checkpoint50i_post, checkpoint48g_post, checkpoint51, checkpoint50, checkpoint50d_post, checkpoint51b_pre, checkpoint50g_post, checkpoint50b_post, checkpoint50f_post, checkpoint50a_post, checkpoint50f_pre, checkpoint51d_post, checkpoint51a_post, checkpoint50e_pre
Changes since 1.3: +12 -15 lines 
Merging from release1_p12:
o Modifications for using pkg/exf with pkg/seaice
  - improved description of the various forcing configurations
  - added basic radiation bulk formulae to pkg/exf
  - units/sign fix for evap computation in exf_getffields.F
  - updated verification/global_with_exf/results/output.txt
o Added pkg/sbo for computing IERS Special Bureau for the Oceans
  (SBO) core products, including oceanic mass, center-of-mass,
  angular, and bottom pressure (see pkg/sbo/README.sbo).
o Lower bound for viscosity/diffusivity in pkg/kpp/kpp_routines.F
  to avoid negative values in shallow regions.
  - updated verification/natl_box/results/output.txt
  - updated verification/lab_sea/results/output.txt
o MPI gather, scatter: eesupp/src/gather_2d.F and scatter_2d.F
o Added useSingleCpuIO option (see PARAMS.h).
o Updated useSingleCpuIO option in mdsio_writefield.F to
  work with multi-field files, e.g., for single-file pickup.
o pkg/seaice:
  - bug fix in growth.F: QNET for no shortwave case
  - added HeffFile for specifying initial sea-ice thickness
  - changed SEAICE_EXTERNAL_FLUXES wind stress implementation
o Added missing /* */ to CPP comments in pkg/seaice, pkg/exf,
  kpp_transport_t.F, forward_step.F, and the_main_loop.F
o pkg/seaice:
  - adjoint-friendly modifications
  - added a SEAICE_WRITE_PICKUP at end of the_model_main.F
1 #include "EXF_CPPOPTIONS.h"
2
3 subroutine exf_init_runoff( mythid )
4
5 c ==================================================================
6 c SUBROUTINE exf_init_runoff
7 c ==================================================================
8 c
9 c o
10 c
11 c started: Ralf.Giering@FastOpt.de 25-Mai-2000
12 c mods for pkg/seaice: menemenlis@jpl.nasa.gov 20-Dec-2002
13 c
14 c ==================================================================
15 c SUBROUTINE exf_init_runoff
16 c ==================================================================
17
18 implicit none
19
20 c == global variables ==
21
22 #include "EEPARAMS.h"
23 #include "SIZE.h"
24 #include "GRID.h"
25
26 #include "exf_param.h"
27 #include "exf_fields.h"
28 #include "exf_constants.h"
29
30 c == routine arguments ==
31
32 integer mythid
33
34 #if (defined (ALLOW_RUNOFF) || defined (ALLOW_SEAICE))
35
36 c == local variables ==
37
38 integer bi, bj
39 integer i, j
40
41 c == end of interface ==
42
43 do bj = mybylo(mythid), mybyhi(mythid)
44 do bi = mybxlo(mythid), mybxhi(mythid)
45 do j = 1, sny
46 do i = 1, snx
47 runoff(i,j,bi,bj) = 0. _d 0
48 enddo
49 enddo
50 enddo
51 enddo
52
53 if ( runofffile .NE. ' ' ) then
54 call mdsreadfield( runofffile, exf_iprec, exf_yftype, 1
55 & , runoff, 1, mythid
56 & )
57 if (exf_yftype .eq. 'RL') then
58 call exf_filter_rl( runoff, runoffmask, mythid )
59 else
60 call exf_filter_rs( runoff, runoffmask, mythid )
61 end if
62
63 c-- scale runoff array
64 do bj = mybylo(mythid), mybyhi(mythid)
65 do bi = mybxlo(mythid), mybxhi(mythid)
66 do j = 1, sny
67 do i = 1, snx
68 if (runoff(i,j,bi,bj).gt.20.) runoff(i,j,bi,bj)=0.
69 runoff(i,j,bi,bj) = exf_inscal_runoff*runoff(i,j,bi,bj)
70 enddo
71 enddo
72 enddo
73 enddo
74 endif
75
76 #endif
77
78 end
  ViewVC Help
Powered by ViewVC 1.1.22