Changeset 1404 for trunk/UTIL/SPECTRA

Timestamp:

Mar 31, 2015, 6:03:30 PM (10 years ago)

Author:

milmd

Message:

Add the opportunity to create some harmonical wind fields for testing. Spherepack library can also be compiled with ifort.

Location:

trunk/UTIL/SPECTRA

Files:

• makefile (modified) (1 diff)
• readme (modified) (2 diffs)
• spectra_analysis.f90 (modified) (8 diffs)
• test_harmonic.f90 (added)

trunk/UTIL/SPECTRA/makefile

@@ -12 +12 @@
 TMP=$(SRCS:=~)
 
+all: spectra_analysis test_harmonic
+
 %.o: %.f90
         $(FC) $(FFLAGS) -c $<

trunk/UTIL/SPECTRA/readme

@@ -1 @@
->>>> install (64b) spherepack 3.2 library with pgf90
-(g95,ifort,gfortran... not  possible yet)
+>>>> install spectra tool and spherepack 3.2 library with pgf90 or ifort
+(g95,gfortran... not  possible yet)
 
 wget -c https://www2.cisl.ucar.edu/sites/default/files/spherepack3.2.tar
@@ -6 +6 @@
 cd spherepack3.2
 
-> install
-for pgf90, after discussion with Richard A. Valent (valent@ucar.edu), change line 12 of make.inc in spherepack3.2 folder:
-F90 := pgf90 -module ../lib -I../lib
+>>>> install spherepack 3.2 library
+>for pgf90, after discussion with Richard A. Valent (valent@ucar.edu), change line 12 of make.inc in spherepack3.2 folder:
+  F90 := pgf90 -module ../lib -I../lib
 and replace by
-F90 := pgf90 -r8 -module ../lib -I../lib
+  F90 := pgf90 -r8 -module ../lib -I../lib
 to run in double precision (better).
 
-make all > log_install 2>&1
+>for ifort:
+change line 10 to 20 of make.inc
+  ifeq ($(PGI),pgf90-Warning-No files to process)
+
+    F90 := pgf90 -module ../lib -I../lib
+    CPP := pgf90 -E
+
+  else
+
+    F90 := g95 -DG95 -g -fmod=../lib -I../lib 
+    CPP := g95 -E -DG95
+
+  endif
+and replace by
+  F90 := ifort -r8 -module ../lib -I../lib
+  CPP := ifort -fpp
+change line 16 of test/Makefile
+        $(F90) $< -o $@ -L../lib -l spherepack
+and replace by
+        $(F90) $< -o $@ -L../lib -lspherepack
+        
+
+> make all > log_install 2>&1
 
 > OK 'lib' ,... are in /planeto/milmd/library/spherepack/spherepack3.2_levan_pgf90
 
 
->>>> compile spectra_analysis program with pgf90 or g95
+>>>> compile spectra_analysis program with pgf90 or ifort
 
 cd ..
 
 change the first 3 lines of the makefile: path of spherepack library you have just installed, path of netcdf libray, compiler
-netcdfpath=/planeto/milmd/library/netcdf/netcdf-4.0.1_levan_pgf90
-spherepackpath=/planeto/milmd/library/spherepack/spherepack3.2_levan_pgf90
-FC=pgf90
+netcdfpath=/planeto/milmd/library/netcdf/netcdf-4.0.1_levan_pgf90 (on levan)
+netcdfpath=/smplocal/pub/NetCDF/4.1.3/seq (on ada)
+netcdfpath=/usr/local (on gnome)
+spherepackpath=/path/to/spherepack3.2
+FC=pgf90 or ifort
 
 make spectra_analysis
+
+>>>>> test on some harmonic cases
+
+> compile test program
+make test_analysis
+
+> create harmonic winds
+./test_analysis
+
+> compute kinetic energy spectrum for one harmonic
+./spectra_analysis harmonic_64x48_lmdz_22.nc -alt none -time none -o harmonic_64x48_lmdz_22_spectra
+

trunk/UTIL/SPECTRA/spectra_analysis.f90

@@ -132 +132 @@
       print*,'[-lat str]\t: name of latitude field in netcdf file'
       print*,'[-lon str]\t: name of longitude field in netcdf file'
-      print*,'[-alt str]\t: name of altitude field in netcdf file. \'none\' if no altitude axis.'
-      print*,'[-time str]\t: name of time field in netcdf file. \'none\' if no time axis.'
+      print*,'[-alt str]\t: name of altitude field in netcdf file. ''none'' if no altitude axis.'
+      print*,'[-time str]\t: name of time field in netcdf file. ''none'' if no time axis.'
       print*,'[-divrot]\t: total, divergence and vorticity spectra coefficient in output file'
       stop 'End help'
@@ -166 +166 @@
 else
   print*,'latitude dimension not found!'
-  stop ' must have this... use \'-lat latitudefieldname\' option.'
+  stop ' must have this... use ''-lat latitudefieldname'' option.'
 end if
 
@@ -180 +180 @@
 else
   print*,'longitude dimension not found!'
-  stop ' must have this... use \'-lon longitudefieldname\' option.'
+  stop ' must have this... use ''-lon longitudefieldname'' option.'
 end if
 
@@ -197 +197 @@
   else 
     print*,'altitude dimension not found!'
-    stop ' if no altitude axis, use \'-alt none\' option...'
+    stop ' if no altitude axis, use ''-alt none'' option...'
   end if
 end if
@@ -215 +215 @@
   else 
     print*,'time dimension not found!'
-    stop ' if no time axis, use \'-time none\' option...'
+    stop ' if no time axis, use ''-time none'' option...'
   end if
 end if
@@ -295 +295 @@
 else
   print*,'zonal wind variable not found!'
-  stop ' must have this... use \'-u zonalwindfieldname\' option.'
+  stop ' must have this... use ''-u zonalwindfieldname'' option.'
 end if
 
@@ -324 +324 @@
 else
   print*,'meridional wind variable not found!'
-  stop ' must have this... use \'-v meridionalwindfieldname\' option.'
+  stop ' must have this... use ''-v meridionalwindfieldname'' option.'
 end if
 
@@ -542 +542 @@
   number_config = (t-1)*n_indz+z
   if (is_div_rot) then
-    write(10,'(e12.6E2,a4,e12.6E2,a4,e12.6E2,a6)',advance='no') spectra_config(j,number_config),' ',&
+    write(10,'(e13.6E2,a4,e13.6E2,a4,e13.6E2,a6)',advance='no') spectra_config(j,number_config),' ',&
                         spectra_b_config(j,number_config),' ',spectra_c_config(j,number_config),' '
   else
-    write(10,'(e12.6E2,a38)',advance='no') spectra_config(j,number_config),' '
+    write(10,'(e13.6E2,a38)',advance='no') spectra_config(j,number_config),' '
   end if
   end do