Changeset 3572 in Sophya for trunk/SophyaExt/XAstroPack/tstmollweide.cc
- Timestamp:
- Feb 7, 2009, 10:50:34 PM (17 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/SophyaExt/XAstroPack/tstmollweide.cc
r2743 r3572 20 20 a = 2.*tvrai+sin(2.*tvrai); 21 21 lat = asin(a/M_PI); 22 tcalc = MollWeide_XY(0.,lat, &x,&y);22 tcalc = MollWeide_XY(0.,lat,x,y); 23 23 if(fabs(tcalc-tvrai)>tll) {tll = fabs(tcalc-tvrai); txy=tvrai;} 24 24 if(file) fprintf(file,"%e %e %e %e\n",tvrai,tcalc,a,tcalc-tvrai); … … 28 28 29 29 for(lon=0.;lon<360.;lon+=2.5) for(lat=-90.;lat<=90.1;lat+=2.5) { 30 txy = MollWeide_XY(lon*M_PI/180.,lat*M_PI/180., &x,&y);31 tll = MollWeide_LL(x,y, &longitude,&latitude);30 txy = MollWeide_XY(lon*M_PI/180.,lat*M_PI/180.,x,y); 31 tll = MollWeide_LL(x,y,longitude,latitude); 32 32 longitude *= 180./M_PI; 33 33 latitude *= 180./M_PI;
Note:
See TracChangeset
for help on using the changeset viewer.