pkg/fizhi/fizhi_alarms.F

C $Header: /u/gcmpack/MITgcm/pkg/fizhi/fizhi_alarms.F,v 1.12 2005/03/08 19:54:26 molod Exp $
C $Name:  $

#include "FIZHI_OPTIONS.h"
      subroutine fizhi_alarms (nymdb,nhmsb,deltaT)

C***********************************************************************
C  Purpose
C  -------
C     Driver to Set Internal Model Alarms
C
C  Argument Description
C  --------------------
C     nymdb ..... Begining Date to Start Alarm Clock
C     nhmsb ..... Begining Time to Start Alarm Clock
C     ndt ....... Model Timestep
C
C     Note: For ndt = 0, all Intermal Model Alarms are set to 1 second
C
C***********************************************************************
      implicit none

c Input Parameters
c ----------------
      integer  nymdb,nhmsb
      _RL deltaT

      integer  myid

c Physics Frequencies
c -------------------
      integer ndt
      integer  nhmsf, m,n, nd
      integer  ndlw, ndsw, ndturb, ndmoist

      myid = 1

      ndt = deltaT

      if( ndt.ne.0 ) then

C Longwave  Radiation Frequency
      ndlw = 003000
C Shortwave Radiation Frequency
      ndsw = 003000

c Turbulence
c ----------
      ndturb = 2*int( max(1.0,(300./iabs(ndt)+0.5)/2) ) * iabs(ndt)

      n  =  1
      m  = -1
      nd = ndturb
      do while ( mod(3*3600,nd).ne.0 )
      nd = ndturb + 2*n*m*ndt
      m  = -m
      n  =  n+1
      if( nd.lt.ndt ) then
          if( myid.eq.1 ) then
          print *
          print *, 'Turbulence cannot occur with Model Timestep = ',ndt
          print *
          endif
      call my_finalize
      call my_exit (101)
      endif
      enddo
      ndturb  = nhmsf( nd )

c Moist Processes
c ---------------
      ndmoist = int( 600./iabs(ndt) + 0.5 ) * iabs(ndt)

      n  =  1
      m  = -1
      nd = ndmoist
      do while ( mod(3*3600,nd).ne.0 )
      nd = ndmoist + n*m*ndt
      m  = -m
      n  =  n+1
      if( nd.lt.ndt ) then
          if( myid.eq.1 ) then
          print *
          print *, 'Convection cannot occur with Model Timestep = ',ndt
          print *
          endif
      call my_finalize
      call my_exit (101)
      endif
      enddo
      ndmoist = nhmsf( nd )

      else

c For NDT = 0, all Alarms set to 1 second
c ---------------------------------------
C Longwave  Radiation Frequency
      ndlw    = 000001
C Shortwave Radiation Frequency
      ndsw    = 000001
C Turbulence          Frequency
      ndturb  = 000001
C Moist Processes     Frequency
      ndmoist = 000001

      endif

c Set Alarm Clocks
c ----------------

      call set_alarm ( 'radsw',nymdb,nhmsb,ndsw    )
      call set_alarm ( 'radlw',nymdb,nhmsb,ndlw    )
      call set_alarm ( 'turb' ,nymdb,nhmsb,ndturb  )
      call set_alarm ( 'moist',nymdb,nhmsb,ndmoist )
      call set_alarm ( 'pnt',nymdb,nhmsb,0 )

      return
      end
1	C $Header: /u/gcmpack/MITgcm/pkg/fizhi/fizhi_alarms.F,v 1.12 2005/03/08 19:54:26 molod Exp $
2	C $Name: $
3
4	#include "FIZHI_OPTIONS.h"
5	subroutine fizhi_alarms (nymdb,nhmsb,deltaT)
6
7	C***********************************************************************
8	C Purpose
9	C -------
10	C Driver to Set Internal Model Alarms
11	C
12	C Argument Description
13	C --------------------
14	C nymdb ..... Begining Date to Start Alarm Clock
15	C nhmsb ..... Begining Time to Start Alarm Clock
16	C ndt ....... Model Timestep
17	C
18	C Note: For ndt = 0, all Intermal Model Alarms are set to 1 second
19	C
20	C***********************************************************************
21	implicit none
22
23	c Input Parameters
24	c ----------------
25	integer nymdb,nhmsb
26	_RL deltaT
27
28	integer myid
29
30	c Physics Frequencies
31	c -------------------
32	integer ndt
33	integer nhmsf, m,n, nd
34	integer ndlw, ndsw, ndturb, ndmoist
35
36	myid = 1
37
38	ndt = deltaT
39
40	if( ndt.ne.0 ) then
41
42	C Longwave Radiation Frequency
43	ndlw = 003000
44	C Shortwave Radiation Frequency
45	ndsw = 003000
46
47	c Turbulence
48	c ----------
49	ndturb = 2int( 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 + 2nm*ndt
56	m = -m
57	n = n+1
58	if( nd.lt.ndt ) then
59	if( myid.eq.1 ) then
60	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	ndturb = nhmsf( nd )
69
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 + nmndt
79	m = -m
80	n = n+1
81	if( nd.lt.ndt ) then
82	if( myid.eq.1 ) then
83	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	C Longwave Radiation Frequency
98	ndlw = 000001
99	C Shortwave Radiation Frequency
100	ndsw = 000001
101	C Turbulence Frequency
102	ndturb = 000001
103	C Moist Processes Frequency
104	ndmoist = 000001
105
106	endif
107
108	c Set Alarm Clocks
109	c ----------------
110
111	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	call set_alarm ( 'pnt',nymdb,nhmsb,0 )
116
117	return
118	end