/[MITgcm]/MITgcm/pkg/fizhi/fizhi_alarms.F
ViewVC logotype

Annotation of /MITgcm/pkg/fizhi/fizhi_alarms.F

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

Revision 1.13 - (hide annotations) (download)
Wed Apr 1 19:55:07 2009 UTC (16 years, 8 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint64y, checkpoint64x, checkpoint64z, checkpoint64q, checkpoint64p, checkpoint64s, checkpoint64r, checkpoint64u, checkpoint64t, checkpoint64w, checkpoint64v, checkpoint64i, checkpoint64h, checkpoint64k, checkpoint64j, checkpoint64m, checkpoint64l, checkpoint64o, checkpoint64n, checkpoint64a, checkpoint64c, checkpoint64b, checkpoint64e, checkpoint64d, checkpoint64g, checkpoint64f, checkpoint63p, checkpoint63q, checkpoint63r, checkpoint63s, checkpoint63l, checkpoint63m, checkpoint63n, checkpoint63o, checkpoint63h, checkpoint63i, checkpoint63j, checkpoint63k, checkpoint63d, checkpoint63e, checkpoint63f, checkpoint63g, checkpoint63a, checkpoint63b, checkpoint63c, checkpoint64, checkpoint65, checkpoint62, checkpoint63, checkpoint66g, checkpoint66f, checkpoint66e, checkpoint66d, checkpoint66c, checkpoint66b, checkpoint66a, checkpoint66o, checkpoint66n, checkpoint66m, checkpoint66l, checkpoint66k, checkpoint66j, checkpoint66i, checkpoint66h, checkpoint65z, checkpoint65x, checkpoint65y, checkpoint65r, checkpoint65s, checkpoint65p, checkpoint65q, checkpoint65v, checkpoint65w, checkpoint65t, checkpoint65u, checkpoint65j, checkpoint65k, checkpoint65h, checkpoint65i, checkpoint65n, checkpoint65o, checkpoint65l, checkpoint65m, checkpoint65b, checkpoint65c, checkpoint65a, checkpoint65f, checkpoint65g, checkpoint65d, checkpoint65e, checkpoint62c, checkpoint62b, checkpoint62a, checkpoint62g, checkpoint62f, checkpoint62e, checkpoint62d, checkpoint62k, checkpoint62j, checkpoint62i, checkpoint62h, checkpoint62o, checkpoint62n, checkpoint62m, checkpoint62l, checkpoint62s, checkpoint62r, checkpoint62q, checkpoint62p, checkpoint62w, checkpoint62v, checkpoint62u, checkpoint62t, checkpoint62z, checkpoint62y, checkpoint62x, checkpoint61n, checkpoint61o, checkpoint61l, checkpoint61m, checkpoint61v, checkpoint61w, checkpoint61t, checkpoint61u, checkpoint61r, checkpoint61s, checkpoint61p, checkpoint61q, checkpoint61z, checkpoint61x, checkpoint61y, HEAD
Changes since 1.12: +17 -27 lines 
remove unnecessary MPI stuff
1 jmc 1.13 C $Header: /u/gcmpack/MITgcm/pkg/fizhi/fizhi_alarms.F,v 1.12 2005/03/08 19:54:26 molod Exp $
2 molod 1.2 C $Name: $
3    
4 molod 1.6 #include "FIZHI_OPTIONS.h"
5 molod 1.8 subroutine fizhi_alarms (nymdb,nhmsb,deltaT)
6 molod 1.1
7 jmc 1.13 C***********************************************************************
8     C Purpose
9     C -------
10 molod 1.1 C Driver to Set Internal Model Alarms
11     C
12 jmc 1.13 C Argument Description
13     C --------------------
14     C nymdb ..... Begining Date to Start Alarm Clock
15     C nhmsb ..... Begining Time to Start Alarm Clock
16 molod 1.1 C ndt ....... Model Timestep
17     C
18     C Note: For ndt = 0, all Intermal Model Alarms are set to 1 second
19     C
20 jmc 1.13 C***********************************************************************
21 molod 1.1 implicit none
22    
23     c Input Parameters
24     c ----------------
25 molod 1.8 integer nymdb,nhmsb
26     _RL deltaT
27 molod 1.1
28 molod 1.5 integer myid
29 molod 1.1
30     c Physics Frequencies
31     c -------------------
32 molod 1.8 integer ndt
33 molod 1.1 integer nhmsf, m,n, nd
34     integer ndlw, ndsw, ndturb, ndmoist
35 jmc 1.13
36 molod 1.3 myid = 1
37 molod 1.8
38     ndt = deltaT
39    
40 molod 1.1 if( ndt.ne.0 ) then
41    
42 molod 1.4 C Longwave Radiation Frequency
43 jmc 1.13 ndlw = 003000
44 molod 1.4 C Shortwave Radiation Frequency
45 jmc 1.13 ndsw = 003000
46 molod 1.1
47     c Turbulence
48     c ----------
49     ndturb = 2*int( max(1.0,(300./iabs(ndt)+0.5)/2) ) * iabs(ndt)
50    
51     n = 1
52     m = -1
53     nd = ndturb
54     do while ( mod(3*3600,nd).ne.0 )
55     nd = ndturb + 2*n*m*ndt
56     m = -m
57     n = n+1
58     if( nd.lt.ndt ) then
59 molod 1.3 if( myid.eq.1 ) then
60 molod 1.1 print *
61     print *, 'Turbulence cannot occur with Model Timestep = ',ndt
62     print *
63     endif
64     call my_finalize
65     call my_exit (101)
66     endif
67     enddo
68 jmc 1.13 ndturb = nhmsf( nd )
69 molod 1.1
70     c Moist Processes
71     c ---------------
72     ndmoist = int( 600./iabs(ndt) + 0.5 ) * iabs(ndt)
73    
74     n = 1
75     m = -1
76     nd = ndmoist
77     do while ( mod(3*3600,nd).ne.0 )
78     nd = ndmoist + n*m*ndt
79     m = -m
80     n = n+1
81     if( nd.lt.ndt ) then
82 molod 1.3 if( myid.eq.1 ) then
83 molod 1.1 print *
84     print *, 'Convection cannot occur with Model Timestep = ',ndt
85     print *
86     endif
87     call my_finalize
88     call my_exit (101)
89     endif
90     enddo
91     ndmoist = nhmsf( nd )
92    
93     else
94    
95     c For NDT = 0, all Alarms set to 1 second
96     c ---------------------------------------
97 molod 1.4 C Longwave Radiation Frequency
98 jmc 1.13 ndlw = 000001
99 molod 1.4 C Shortwave Radiation Frequency
100 jmc 1.13 ndsw = 000001
101 molod 1.4 C Turbulence Frequency
102 jmc 1.13 ndturb = 000001
103 molod 1.4 C Moist Processes Frequency
104 jmc 1.13 ndmoist = 000001
105 molod 1.1
106     endif
107    
108     c Set Alarm Clocks
109     c ----------------
110 molod 1.9
111 molod 1.1 call set_alarm ( 'radsw',nymdb,nhmsb,ndsw )
112     call set_alarm ( 'radlw',nymdb,nhmsb,ndlw )
113     call set_alarm ( 'turb' ,nymdb,nhmsb,ndturb )
114     call set_alarm ( 'moist',nymdb,nhmsb,ndmoist )
115 molod 1.7 call set_alarm ( 'pnt',nymdb,nhmsb,0 )
116 molod 1.1
117     return
118     end
  ViewVC Help
Powered by ViewVC 1.1.22