1 |
#!/bin/csh -fe |
2 |
# RCSid $Id$ |
3 |
# |
4 |
# Test solar calculations |
5 |
# |
6 |
# Produces: |
7 |
# LAT LON SM TIME DAY MONTH SALT2 SAZI2 |
8 |
# all angles are degrees, as taken by sun2.cal |
9 |
cnt 1000 | rcalc -f sun2.cal -e YEAR=2020 \ |
10 |
-e 'LAT=160*rand(.359*recno-10)-80' \ |
11 |
-e 'LON=360*rand(3.561*recno+16.6)-180' \ |
12 |
-e 'SM=floor((LON+7.5)/15)*15' \ |
13 |
-e 'TIME=14*rand(-7.5858*recno-71)+5' \ |
14 |
-e 'MONTH=floor(1+11.99*rand(recno*.785+5.5))' \ |
15 |
-e 'ndays(m):select(m,31,28,31,30,31,30,31,31,30,31,30,31)' \ |
16 |
-e 'DAY=floor(1+.99*rand(recno*-71+9)*ndays(MONTH))' \ |
17 |
-e '$1=LAT;$2=LON;$3=SM;$4=TIME;$5=DAY;$6=MONTH;$7=SALT;$8=SAZI' \ |
18 |
> tests1.txt |
19 |
# Produces: |
20 |
# LAT LON SM TIME DAY MONTH SALT2 SAZI2 SALT SAZI SALT-SALT2 SAZI-SAZI2 |
21 |
rcalc -f sun.cal -e 'DEG:PI/180' -e 'RLAT=$1*DEG;RLON=$2*DEG;RSM=$3*DEG' \ |
22 |
-e 'TIME=$4;DAY=$5;MONTH=$6' \ |
23 |
-e '$1=$1;$2=$2;$3=$3;$4=$4;$5=$5;$6=$6;$7=$7;$8=$8' \ |
24 |
-e '$9=SALT;$10=SAZI;$11=SALT-$7;$12=SAZI-$8' tests1.txt \ |
25 |
> tests2.txt |