COMBLOCK.GCM COMmon BLOCK definitions for run C070 2001/02/12 Contnets of PARAMETER statement IM = number of grid cells in east-west direction JM = number of grid cells in north-south direction LMA = number of atmospheric dynamical layers LMO = maximum number of ocean layers LMSI = number of thermal layers for sea ice or lake ice LMGI = number of thermal layers for glacial ice LMG = number of ground layers NMST = number of straits in ocean layer KAJ = number of J atmospheric diagnostics in AJ-FJ KAJL = number of JxL atmospheric diagnostics in AJL KAIJ = number of IxJ atmospheric diagnostics in AIJ KAIJL = number of IxJxL atmospheric diagnostics in AIJL KOIJL = number of IxJxL ocean diagnostics in OIJL KOLNST= number of LxN strait diagnostics in OLNST KCON = number of conservation quantities in CONSRV TWOPI = 2*pi = 6.283185307179586477 SDAY (s/day) = number of seconds per day = 86400 RDRY (J/C*kg) = gas constant for dry air = 287 RVAP (J/C*kg) = gas constant for water vapor = 461.5 RADIUS (m) = constant radius of the Earth = 6375000 GRAV (m/s^2) = Earth's mean effective gravitational acceleration = 9.81 SHCD (J/C*kg) = specific heat capacity of dry air = 1003.5 SHCV (J/C*kg) = specific heat capacity of water vapor = 1911 SHCW (J/C*kg) = specific heat capacity of liquid water = 4185 SHCI (J/C*kg) = specific heat capacity of ice = 2060 RKAP = exponent of exner function = RDRY/SHCD ELHE (J/kg) = latent heat of vaporization at 0 C = 2500000 ELHM (J/kg) = latent heat of melting at 0 C = 334000 ELHS (J/kg) = latent heat of sublimation at 0 C = 2834000 RRT3 = one over the square root of three = .57735 ACE1I (kg/m^2) = ice in first mass layer for all ice types = 182 ACE2GI (kg/m^2) = ice in second mass layer for glacial ice = 7280 Contents of COMMON /PARMCB/ DTS (s) = number of seconds per source time step DTA (s) = time step used in atmospheric momentum equation = DTS/NDYNA DTO (s) = time step used in ocean momentum equation = DTS/NDYNO MDRYA (kg/m^2) = global mean atmospheric mass = 10030.581 MTROP (kg/m^2) = unused MSTRAT (kg/m^2) = global atmospheric mass above dynamical top = 101.937 IDT = primary model time counter, DTS steps since IYEAR0 Jan 1, hr 0 IDT0 = primary model time counter at start of run IDAY = current model time, days since IYEAR0 Jan 1 IDAY0 = model time at start of run, days since IYEAR0 Jan 1 IHOUR = current model time, hours since IYEAR0 Jan 1, hr 0 IHOUR0 = model time at start of run, hours since IYEAR0 Jan 1, hr 0 IYEAR = year of current model time (A.D.) IYEAR0 = base year for model time (A.D.) JDATE = day of the month for current model time, 1 to 31 JDATE0 = day of the month for start of diagnostic accumulation period JDAY = Julian day of year for current model time = 1 + mod(IDAY,365) JDAY0 = Julian day of year for start of diagnostic accumulation period JHOUR = hour of the day for current model time = mod(IHOUR,24) JHOUR0 = hour of the day for start of diagnostic accumulation period JMON = character month of the year for current model time JMON0 = character month of the year for start of diagnostc accum period JYEAR = year of current model time = IYEAR JYEAR0 = year for start of diagnostic accumulation period LMCM = highest layer from which moist convection can start LSSM = highest layer from which large scale condenstion occurs LS1 = first atmospheric layer in the stratosphere LTM = highest atmospheric layer in the troposphere = LS1 - 1 MCNDS = computer time (.01 s) for condensation and precipitation MDIAG = computer time (.01 s) for diagnostic accumulating subroutines MDYNA = computer time (.01 s) for atmospheric dynamics MDYNO = computer time (.01 s) for ocean dynamics MELSE = computer time (.01 s) for initialization and writing output MRAD = computer time (.01 s) for radiation MSRCO = computer time (.01 s) for ocean source terms and straits MSURF = computer time (.01 s) for surface, ground, glacial ice, lakes NBERG = current number of ice bergs in the ocean NDAY = number of DTS steps per day = 86400/DTS NDIAGA = number of DTA steps between calling DIAGA NDIAGD = unused NDIAGE = unused NDIAGF = unused NDIAGO = number of DTS steps between accumulating ocean prognostic vars NDIAGS = number of surface steps between accumulating surface diagnostic NDISK = number of DTS steps between writing the ReStartFile NDYNA = number of atmo dynamic steps per source term call = DTS/DTA NDYNO = number of ocean dynamic steps per source term call = DTS/DTO NFILXA = number of DTS steps between applying zonal filter to UA,VA NFILXO = number of DTS steps between applying zonal filter to UO,VO NFILYO = number of DTS steps between applying merid filter to UO,VO NFILYO = number of DTS steps between applying merid filter to UO,VO NRAD = number of DTS steps between calling radiation routines NSLP = unused NSURF = number of surface steps per DTS step NTAPE = number of DTS steps between writing auxilliary tape QATMO = .true. for coupled or atmosphere only simulations QOCEAN = .true. for coupled or ocean only simulations QGRND = .false. for ocean only simulations without land subroutines QRAD = whether radiation should be called this source term step QPK = whether P^RKAP is consistent with atmospheric mass field QCHECK = whether CHECKT should be called to check prognostic variables Contents of COMMON /ATMOCB/ MA (kg/m^2) = atmospheric mass UA (m/s) = atmospheric eastward velocity component defined on C grid VA (m/s) = atmospheric northward velocity component defined on C grid H0M (J) = mean atmospheric potential enthalpy HXM (J) = eastward gradient of atmospheric potential enthalpy HYM (J) = northward gradient of atmospheric potential enthalpy HZM (J) = upward vertical gradient of atmospheric potential enthalpy Q0M (kg) = mean atmospheric water vapor QXM (kg) = eastward gradient of atmospheric water vapor QYM (kg) = northward gradient of atmospheric water vapor QZM (kg) = upward vertical gradient of atmospheric water vapor Contents of COMMON /OCENCB/ MO (kg/m^2) = ocean mass UO (m/s) = ocean eastward velocity component defined on C grid VO (m/s) = ocean northward velocity component defined on C grid G0M (J) = mean ocean potential enthalpy GXM (J) = eastward gradient of ocean potential enthalpy GYM (J) = northward gradient of ocean potential enthalpy GZM (J) = downward vertical gradient of ocean potential enthalpy S0M (kg) = mean ocean salt SXM (kg) = eastward gradient of ocean salt SYM (kg) = northward gradient of ocean salt SZM (kg) = downward vertical gradient of ocean salt Contents of COMMON /GICECB/ MGI (kg/m^2) = snow and glacial ice in mass layer 1 HGI (J/m^2) = heat content in glacial ice thermal layers Contents of COMMON /SICECB/ RSI = horizontal ratio of sea ice to ocean RSIX,RSIY = eastward and northward gradients of RSI MSI (kg/m^2) = mass in sea ice mass layers HSI (J/m^2) = heat content in sea ice thermal layers PSI (Pa) = internal sea ice pressure USI (m/s) = eastward sea ice velocity component on C grid VSI (m/s) = northward sea ice velocity component on C grid Contents of COMMON /FIXDCB/ (Read in or calculated in INPUT) FOCEAN = fraction of surface area covered by ocean (should be 0 or 1) FLAKE = fraction of surface area covered by lakes FGRND = fraction of surface area covered by ground FGICE = fraction of surface area covered by glacial or land ice FWATER = fraction of surface area covered by water = FOCEAN + FLAKE FLAND = fraction of surface area covered by land = FGRND + FGICE ZATMO (m) = mean topography that atmosphere sees ZSOLID (m) = mean topography beneath air, water and ice LMM = number of ocean layers for each primary horizontal grid cell LMU = lesser of ocean layers of grid cell and its eastern neighbor LMV = lesser of ocean layers of grid cell and its northern neighbor Contents of COMMON /GEOMCB/ (Calculated in subroutine INPUT) DXYP (m^2) = spherical area of primary grid cell DXP (m) = east-west distance of primary grid cell DYP (m) = north-south distance of primary grid cell DXV (m) = average distance at northern edge of primary grid cell DYV (m) = distance from grid cell centers of north-south neighbors COSV = average cosine of latitude at northern edge of grid cell COSP = average cosine of latitude of primary grid cell SINP = sine of latitude at center of primary grid cell DXYV (m^2) = spherical area of secondary grid cell DXYS (m^2) = half of area of primary grid cell except at poles DXYN (m^2) = half of area of primary grid cell except at poles RAMVS = ratio of areas = DXYS(J) / [DXYN(J-1)+DXYS(J)] RAMVN = ratio of areas = DXYN(J) / [DXYN(J)+DXYS(J+1)] RLAT = latitude in radians at center of primary grid cell SINI = sine of longitude of grid cell center measured from IDL COSI = cosine of longitude of grid cell center measured from IDL byDXYP (m^-2) = reciprocal of DXYP For both atmosphere and ocean models, horizontal velocity and momentum components are defined on the C grid. UA and UO are the average velocity between a grid cell and its eastern neighbor. VA and VO are the average velocity between a grid cell and its northern neighbor.