| 1 | //MODULE phycst |
|---|
| 2 | // !!====================================================================== |
|---|
| 3 | // !! *** MODULE phycst *** |
|---|
| 4 | // !! Definition of of both ocean and ice parameters used in the code |
|---|
| 5 | // !!===================================================================== |
|---|
| 6 | double rpi = 3.141592653589793;// !: pi |
|---|
| 7 | double rad = 3.141592653589793/180.;//_wp !: conversion from degre into radian |
|---|
| 8 | //double rsmall = 0.5 * EPSILON( 1. );// !: smallest real computer value |
|---|
| 9 | double |
|---|
| 10 | rday = 24.*60.*60. ,// & !: day (s) |
|---|
| 11 | rsiyea ,// & !: sideral year (s) |
|---|
| 12 | rsiday ,// & !: sideral day (s) |
|---|
| 13 | raajj = 365. ,// & !: number of days in one year |
|---|
| 14 | raamo = 12. ,// & !: number of months in one year |
|---|
| 15 | rjjhh = 24. ,// & !: number of hours in one day |
|---|
| 16 | rhhmm = 60. ,// & !: number of minutes in one hour |
|---|
| 17 | rmmss = 60. ,// & !: number of seconds in one minute |
|---|
| 18 | raass ,// & !: number of seconds in one year |
|---|
| 19 | rmoss ,// & !: number of seconds in one month |
|---|
| 20 | rjjss ,// & !: number of seconds in one day |
|---|
| 21 | omega ,// & !: earth rotation parameter |
|---|
| 22 | ra = 6371229. ,// & !: earth radius (meter) |
|---|
| 23 | grav = 9.80665; // !: gravity (m/s2) |
|---|
| 24 | |
|---|
| 25 | double |
|---|
| 26 | rtt = 273.16 ,// & !: triple point of temperature (Kelvin) |
|---|
| 27 | rt0 = 273.15 ,// & !: freezing point of water (Kelvin) |
|---|
| 28 | rt0_snow = 273.15 ,// & !: melting point of snow (Kelvin) |
|---|
| 29 | rt0_ice = 273.05 ,// & !: melting point of ice (Kelvin) |
|---|
| 30 | rau0 = 1020. ,// & !: volumic mass of reference (kg/m3) |
|---|
| 31 | rauw = 1000. ,// & !: volumic mass of pure water (kg/m3) |
|---|
| 32 | rcp = 4.e+3,// & !: ocean specific heat |
|---|
| 33 | ro0cpr; // !: = 1. / ( rau0 * rcp ) |
|---|
| 34 | |
|---|
| 35 | double |
|---|
| 36 | rcdsn = 0.22 ,// & !: conductivity of the snow |
|---|
| 37 | rcdic = 2.034396 ,// & !: conductivity of the ice |
|---|
| 38 | rcpsn = 6.9069e+5,// & !: density times specific heat for snow |
|---|
| 39 | rcpic = 1.8837e+6,// & !: volumetric latent heat fusion of sea ice |
|---|
| 40 | xlsn = 110.121e+6 ,// & !: volumetric latent heat fusion of snow |
|---|
| 41 | xlic = 300.33e+6 ,// & !: volumetric latent heat fusion of ice |
|---|
| 42 | xsn = 2.8e+6 ,// & !: latent heat of sublimation of snow |
|---|
| 43 | rhoic = 900. ,// & !: volumic mass of sea ice (kg/m3) |
|---|
| 44 | rhosn = 330. ,// & !: volumic mass of snow (kg/m3) |
|---|
| 45 | emic = 0.97 ,// & !: emissivity of snow or ice |
|---|
| 46 | sice = 6.0 ,// & !: salinity of ice (psu) |
|---|
| 47 | soce = 34.7 ,// & !: salinity of sea (psu) |
|---|
| 48 | cevap = 2.5e+6 ,// & !: latent heat of evaporation (water) |
|---|
| 49 | srgamma = 0.9 ,// & !: correction factor for solar radiation (Oberhuber, 1974) |
|---|
| 50 | vkarmn = 0.4 ,// & !: von Karman constant |
|---|
| 51 | stefan = 5.67e-8 ;// !: Stefan-Boltzmann constant |
|---|
| 52 | |
|---|
| 53 | void phy_cst(); |
|---|
| 54 | void phy_cst(){ |
|---|
| 55 | rsiyea = 365.25 * rday * 2. * rpi / 6.283076; |
|---|
| 56 | rsiday = rday / ( 1. + rday / rsiyea ); |
|---|
| 57 | omega = 2. * rpi / rsiday ; |
|---|
| 58 | rjjss = rjjhh * rhhmm * rmmss; |
|---|
| 59 | ro0cpr = 1. / ( rau0 * rcp ); |
|---|
| 60 | raajj = 30*raamo; // a revoir |
|---|
| 61 | raass = raajj*rjjss;// a revoir |
|---|
| 62 | rmoss = 30*rjjss;// a revoir |
|---|
| 63 | } |
|---|
| 64 | // !!====================================================================== |
|---|
![(please configure the [header_logo] section in trac.ini)](/nemoyao/chrome/site/your_project_logo.png){kind=logo}
