00001
00002
00003
00004
00005 !> @file
00006 !!
00007 !! Declarations of global variables for MAIC-2.
00008 !!
00009 !! @section Copyright
00010 !!
00011 !! Copyright 2010, 2011 Ralf Greve, Bjoern Grieger, Oliver J. Stenzel
00012 !!
00013 !! @section License
00014 !!
00015 !! This file is part of MAIC-2.
00016 !!
00017 !! MAIC-2 is free software: you can redistribute it and/or modify
00018 !! it under the terms of the GNU General Public License as published by
00019 !! the Free Software Foundation, either version 3 of the License, or
00020 !! (at your option) any later version.
00021 !!
00022 !! MAIC-2 is distributed in the hope that it will be useful,
00023 !! but WITHOUT ANY WARRANTY; without even the implied warranty of
00024 !! MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
00025 !! GNU General Public License for more details.
00026 !!
00027 !! You should have received a copy of the GNU General Public License
00028 !! along with MAIC-2. If not, see <http://www.gnu.org/licenses/>.
00029 !<
00030 !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
00031
00032
00033 !> Declarations of global variables for MAIC-2.
00034 !<------------------------------------------------------------------------------
00035 module maic2_variables
00036
00037 use maic2_types
00038
00039
00040
00041 !> phi_node(l): Latitude of grid point
00042 real(dp), dimension(0:LMAX) :: phi_node
00043
00044 real(dp), dimension(0:LMAX) :: phi_cb1
00045
00046 real(dp), dimension(0:LMAX) :: phi_cb2
00047
00048 real(dp), dimension(0:LMAX) :: cos_phi_cb1
00049
00050 real(dp), dimension(0:LMAX) :: cos_phi_cb2
00051
00052 real(dp), dimension(0:LMAX) :: sin_phi_cb1
00053
00054 real(dp), dimension(0:LMAX) :: sin_phi_cb2
00055
00056 real(dp), dimension(LMAX) :: dphi
00057
00058 real(dp), dimension(LMAX) :: dphi_inv
00059
00060 real(dp), dimension(0:LMAX) :: temp_surf
00061
00062 real(dp), dimension(0:LMAX) :: temp_surf_amp
00063
00064 real(dp), dimension(0:LMAX) :: temp_co2
00065
00066 real(dp), dimension(0:LMAX) :: p_surf
00067
00068 real(dp), dimension(0:LMAX) :: water
00069
00070 real(dp), dimension(0:LMAX) :: water_new
00071
00072 real(dp), dimension(0:LMAX) :: cond
00073
00074 real(dp), dimension(0:LMAX) :: evap
00075
00076 real(dp), dimension(0:LMAX) :: a_net
00077
00078 real(dp), dimension(0:LMAX) :: H
00079
00080 real(dp), dimension(0:LMAX) :: H_new
00081
00082
00083
00084 !> RHO: Density of ice-dust mixture
00085 real(dp) :: RHO
00086
00087 real(dp) :: RHO_I
00088
00089 real(dp) :: RHO_W
00090
00091 real(dp) :: G
00092
00093 real(dp) :: R
00094
00095 real(dp) :: rho_inv
00096
00097
00098
00099 !> insol_time_min: Minimum time of the data values for the insolation etc.
00100 integer(i4b) :: insol_time_min
00101
00102 integer(i4b) :: insol_time_stp
00103
00104 integer(i4b) :: insol_time_max
00105
00106 integer(i4b) :: itercount
00107
00108 !> insol_ma_90(n): Data for the mean-annual north- or south-polar insolation
00109 real(dp), dimension(0:100000) :: insol_ma_90
00110
00111 real(dp), dimension(0:100000) :: obl_data
00112
00113 real(dp), dimension(0:100000) :: ecc_data
00114
00115 !> (= 360 deg - Ls of perihelion )
00116 real(dp), dimension(0:100000) :: ave_data
00117
00118 !> = eccentricity
00119 !> *sin(Laskar's longitude of perihelion from moving equinox),
00120 !> ( where Laskar's longitude of perihelion from moving equinox
00121 !> = Ls of perihelion - 180 deg )
00122 real(dp), dimension(0:100000) :: cp_data
00123
00124 !> psi_tab(n): True anomalies (orbital positions with respect to perihelion)
00125 !> over a Martian year
00126 real(dp), dimension(0:NTIME) :: psi_tab
00127
00128 !> diff_aux(l): Auxiliary quantity needed for the diffusional transport
00129 real(dp), dimension(0:LMAX) :: diff_aux
00130
00131 !> n_output: Number of specified times for data output
00132 integer(i4b) :: n_output
00133
00134 real(dp) :: dtime_out
00135
00136 real(dp), dimension(100) :: time_output
00137
00138 integer(i4b) :: iter_out
00139
00140 integer(i4b), dimension(100) :: iter_output
00141
00142 !> pi: Mathematical constant
00143 real(dp), parameter :: pi = 3.141592653589793_dp
00144
00145 real(dp), parameter :: pi_inv = 1.0_dp/pi
00146
00147 real(dp), parameter :: pi_180 = pi/180.0_dp
00148
00149 real(dp), parameter :: pi_180_inv = 180.0_dp/pi
00150
00151 real(dp), parameter :: eps = 1.0e-05_dp
00152
00153 end module maic2_variables
00154