| 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 |
