Calling UCVM With Fortran - Revision history

Davidgil: Undo revision 16775 by Davidgil (talk)

2016-06-22T23:46:55Z

Undo revision 16775 by Davidgil (talk)

Davidgil: /* Example Code */

2016-06-22T23:46:02Z

Example Code

Davidgil: /* Calling UCVM With Fortran */

2013-09-10T15:45:57Z

Calling UCVM With Fortran

Davidgil: /* Example Code */

2013-07-02T06:01:43Z

Example Code

Davidgil: /* Example Code */

2013-07-02T05:59:33Z

Example Code

Davidgil: /* Calling UCVM From Fortran */

2013-06-20T22:51:34Z

Calling UCVM From Fortran

Davidgil: /* Calling UCVM */

2013-06-20T22:39:59Z

Calling UCVM

Davidgil: /* Example Code */

2013-06-20T22:39:28Z

Example Code

Davidgil: /* Example Code */

2013-06-20T22:39:18Z

Example Code

Davidgil: /* Example Code */

2013-06-20T22:39:05Z

Example Code

@@ Line 14: / Line 14: @@
 == Example Code ==
-        program example
-−
-c       UCVM Configuration Location
+                program example
-        CHARACTER(LEN=80) ucvmconf
+        &nbsp;
-c       Model Name
+        c       UCVM Configuration Location
-       	CHARACTER(LEN=8) model
+                CHARACTER(LEN=80) ucvmconf
-c       Number of points we're passing to ucvm_query
+        c       Model Name
-        INTEGER pts
+                CHARACTER(LEN=8) model
+        c       Number of points we're passing to ucvm_query
-c       The UCVM point data structure.
+                INTEGER pts
-c       coord(1) is longitude
+        &nbsp;
-c       coord(2) is latitutde
+        c       The UCVM point data structure.
-c       coord(3) is depth
+        c       coord(1) is longitude
-        TYPE :: ucvm_point_t
+        c       coord(2) is latitutde
-		REAL*8 coord(3)
+        c       coord(3) is depth
-        END TYPE ucvm_point_t
+                TYPE :: ucvm_point_t
+                       REAL*8 coord(3)
-c       Generic property structure
+                END TYPE ucvm_point_t
-c       Source is where it comes from
+        &nbsp;
-c       vp is P-wave velocity in m/s
+        c       Generic property structure
-c       vs is S-wave velocity in m/s
+        c       Source is where it comes from
-c       rho is density in kg/m^3
+        c       vp is P-wave velocity in m/s
-        TYPE :: ucvm_prop_t
+        c       vs is S-wave velocity in m/s
-		INTEGER source
+        c       rho is density in kg/m^3
-                REAL*8 vp
+                TYPE :: ucvm_prop_t
-                REAL*8 vs
+                        INTEGER source
-                REAL*8 rho
+                        REAL*8 vp
-        END TYPE ucvm_prop_t
+                        REAL*8 vs
+                        REAL*8 rho
-c       Returned data structure
+                END TYPE ucvm_prop_t
-        TYPE :: ucvm_data_t
+        &nbsp;
-                REAL*8 surf
+        c       Returned data structure
-                REAL*8 vs30
+                TYPE :: ucvm_data_t
-                REAL*8 depth
+                        REAL*8 surf
-                INTEGER domain
+                        REAL*8 vs30
-                REAL*8 shift_cr
+                        REAL*8 depth
-                REAL*8 shift_gtl
+                        INTEGER domain
-                type(ucvm_prop_t) crust
+                        REAL*8 shift_cr
-                type(ucvm_prop_t) gtl
+                        REAL*8 shift_gtl
-                type(ucvm_prop_t) cmb
+                        type(ucvm_prop_t) crust
-        END TYPE ucvm_data_t
+                        type(ucvm_prop_t) gtl
+                        type(ucvm_prop_t) cmb
-c       For our example we'll query five points
+                END TYPE ucvm_data_t
-        type(ucvm_point_t) point(5)
+        &nbsp;
-c       And we'll get back five sets of material properties
+        c       For our example we'll query five points
-        type(ucvm_data_t) returnData(5)
+                type(ucvm_point_t) point(5)
+        c       And we'll get back five sets of material properties
-c       Number of points is 5.
+                type(ucvm_data_t) returnData(5)
-        pts = 5
+         &nbsp;
+        c       Number of points is 5.
-c       We'll start at -118, 34 at 0 depth and go down by 1000m
+                pts = 5
-c       each step
+        &nbsp;
-        do 10 i = 1, 5
+        c       We'll start at -118, 34 at 0 depth and go down by 1000m
-	        point(i)%coord(1) = -118
+        c       each step
-                point(i)%coord(2) = 34
+                do 10 i = 1, 5
-                point(i)%coord(3) = (i - 1) * 1000
+                        point(i)%coord(1) = -118
-      continue
+                        point(i)%coord(2) = 34
+                        point(i)%coord(3) = (i - 1) * 1000
-c       Where is our configuration file?
+      continue
-        ucvmconf = "/Users/davidgil/ucvm.conf" // CHAR(0)
+        &nbsp;
+        c       Where is our configuration file?
-c       What model are we querying?
+                ucvmconf = "/home/scec-01/davidgil/ucvm.conf" // CHAR(0)
-        model = "cvmsi" // CHAR(0)
+        &nbsp;
+        c       What model are we querying?
-c       Initialize UCVM
+                model = "cvmsi" // CHAR(0)
-        call ucvm_init(ucvmconf)
+        &nbsp;
+        c       Initialize UCVM
-c       Add the model to UCVM
+                call ucvm_init(ucvmconf)
-        call ucvm_add_model(model)
+        &nbsp;
+        c       Add the model to UCVM
-c       Query the model. Note that the number of points is passed
+                call ucvm_add_model(model)
-c       by value, not reference.
+        &nbsp;
-        call ucvm_query(%VAL(pts), point, returnData)
+        c       Query the model. Note that the number of points is passed
+        c       by value, not reference.
-        print *, model, " results for lon -118, lat 34"
+                call ucvm_query(%VAL(pts), point, returnData)
+        &nbsp;
-c       Print out the results.
+                print *, model, " results for lon -118, lat 34"
-        do 20 i = 1, 5
+        &nbsp;
- 		print *, "Depth ", (i - 1) * 1000
+        c       Print out the results.
-                print *, "Vs ", returnData(i)%crust%vs
+                do 20 i = 1, 5
-                print *, "Vp ", returnData(i)%crust%vp
+                        print *, "Depth ", (i - 1) * 1000
-                print *, "Rho ", returnData(i)%crust%rho
+                        print *, "Vs ", returnData(i)%crust%vs
-      continue
+                        print *, "Vp ", returnData(i)%crust%vp
+                        print *, "Rho ", returnData(i)%crust%rho
-c       Close UCVM now that we've queried the points
+      continue
-        call ucvm_finalize()
+        &nbsp;
+        c       Close UCVM now that we've queried the points
-        end
+                call ucvm_finalize()

@@ Line 14: / Line 14: @@
 == Example Code ==
+        program example
-−
-                program example
+c       UCVM Configuration Location
-        &nbsp;
+        CHARACTER(LEN=80) ucvmconf
-        c       UCVM Configuration Location
+c       Model Name
-                CHARACTER(LEN=80) ucvmconf
+       	CHARACTER(LEN=8) model
-        c       Model Name
+c       Number of points we're passing to ucvm_query
-                CHARACTER(LEN=4) model
+        INTEGER pts
-        c       Number of points we're passing to ucvm_query
-                INTEGER pts
+c       The UCVM point data structure.
-        &nbsp;
+c       coord(1) is longitude
-        c       The UCVM point data structure.
+c       coord(2) is latitutde
-        c       coord(1) is longitude
+c       coord(3) is depth
-        c       coord(2) is latitutde
+        TYPE :: ucvm_point_t
-        c       coord(3) is depth
+		REAL*8 coord(3)
-                TYPE :: ucvm_point_t
+        END TYPE ucvm_point_t
-                       REAL*8 coord(3)
-                END TYPE ucvm_point_t
+c       Generic property structure
-        &nbsp;
+c       Source is where it comes from
-        c       Generic property structure
+c       vp is P-wave velocity in m/s
-        c       Source is where it comes from
+c       vs is S-wave velocity in m/s
-        c       vp is P-wave velocity in m/s
+c       rho is density in kg/m^3
-        c       vs is S-wave velocity in m/s
+        TYPE :: ucvm_prop_t
-        c       rho is density in kg/m^3
+		INTEGER source
-                TYPE :: ucvm_prop_t
+                REAL*8 vp
-                        INTEGER source
+                REAL*8 vs
-                        REAL*8 vp
+                REAL*8 rho
-                        REAL*8 vs
+        END TYPE ucvm_prop_t
-                        REAL*8 rho
-                END TYPE ucvm_prop_t
+c       Returned data structure
-        &nbsp;
+        TYPE :: ucvm_data_t
-        c       Returned data structure
+                REAL*8 surf
-                TYPE :: ucvm_data_t
+                REAL*8 vs30
-                        REAL*8 surf
+                REAL*8 depth
-                        REAL*8 vs30
+                INTEGER domain
-                        REAL*8 depth
+                REAL*8 shift_cr
-                        INTEGER domain
+                REAL*8 shift_gtl
-                        REAL*8 shift_cr
+                type(ucvm_prop_t) crust
-                        REAL*8 shift_gtl
+                type(ucvm_prop_t) gtl
-                        type(ucvm_prop_t) crust
+                type(ucvm_prop_t) cmb
-                        type(ucvm_prop_t) gtl
+        END TYPE ucvm_data_t
-                        type(ucvm_prop_t) cmb
-                END TYPE ucvm_data_t
+c       For our example we'll query five points
-        &nbsp;
+        type(ucvm_point_t) point(5)
-        c       For our example we'll query five points
+c       And we'll get back five sets of material properties
-                type(ucvm_point_t) point(5)
+        type(ucvm_data_t) returnData(5)
-        c       And we'll get back five sets of material properties
-                type(ucvm_data_t) returnData(5)
+c       Number of points is 5.
-         &nbsp;
+        pts = 5
-        c       Number of points is 5.
-                pts = 5
+c       We'll start at -118, 34 at 0 depth and go down by 1000m
-        &nbsp;
+c       each step
-        c       We'll start at -118, 34 at 0 depth and go down by 1000m
+        do 10 i = 1, 5
-        c       each step
+	        point(i)%coord(1) = -118
-                do 10 i = 1, 5
+                point(i)%coord(2) = 34
-                        point(i)%coord(1) = -118
+                point(i)%coord(3) = (i - 1) * 1000
-                        point(i)%coord(2) = 34
+      continue
-                        point(i)%coord(3) = (i - 1) * 1000
-      continue
+c       Where is our configuration file?
-        &nbsp;
+        ucvmconf = "/Users/davidgil/ucvm.conf" // CHAR(0)
-        c       Where is our configuration file?
-                ucvmconf = "/home/scec-01/davidgil/ucvm.conf" // CHAR(0)
+c       What model are we querying?
-        &nbsp;
+        model = "cvmsi" // CHAR(0)
-        c       What model are we querying?
-                model = "cvms"
+c       Initialize UCVM
-        &nbsp;
+        call ucvm_init(ucvmconf)
-        c       Initialize UCVM
-                call ucvm_init(ucvmconf)
+c       Add the model to UCVM
-        &nbsp;
+        call ucvm_add_model(model)
-        c       Add the model to UCVM
-                call ucvm_add_model(model)
+c       Query the model. Note that the number of points is passed
-        &nbsp;
+c       by value, not reference.
-        c       Query the model. Note that the number of points is passed
+        call ucvm_query(%VAL(pts), point, returnData)
-        c       by value, not reference.
-                call ucvm_query(%VAL(pts), point, returnData)
+         print *, model, " results for lon -118, lat 34"
-         &nbsp;
-                print *, model, " results for lon -118, lat 34"
+c       Print out the results.
-        &nbsp;
+        do 20 i = 1, 5
-        c       Print out the results.
+ 		print *, "Depth ", (i - 1) * 1000
-                do 20 i = 1, 5
+                print *, "Vs ", returnData(i)%crust%vs
-                        print *, "Depth ", (i - 1) * 1000
+                print *, "Vp ", returnData(i)%crust%vp
-                        print *, "Vs ", returnData(i)%crust%vs
+                print *, "Rho ", returnData(i)%crust%rho
-                        print *, "Vp ", returnData(i)%crust%vp
+      continue
-                        print *, "Rho ", returnData(i)%crust%rho
-      continue
+c       Close UCVM now that we've queried the points
-        &nbsp;
+        call ucvm_finalize()
-        c       Close UCVM now that we've queried the points
-                call ucvm_finalize()
+         end

@@ Line 9: / Line 9: @@
 As an example, suppose we have a Fortran file, example.f, that calls UCVM. We have compiled UCVM with CVM-S and CVM-H. The code to compile example.f would be as follows:
-gfortran example.f -o ./example -L./ucvm-12.2.0/inst/lib -L./cvms/lib -L./cvmh-11.9.1/inst/lib -L./proj-4.7.0/inst/lib -L./euclid3-1.3/libsrc -lucvm -lcvms -lvxapi -lgeo -lproj -letree -fno-underscoring
+gfortran example.f -o ./example -L./ucvm-13.9.0/inst/lib -L./cvms/lib -L./cvmh-11.9.1/inst/lib -L./proj-4.8.0/inst/lib -L./euclid3-1.3/libsrc -lucvm -lcvms -lvxapi -lgeo -lproj -letree -fno-underscoring
 The basic structure of how to call UCVM within Fortran is outlined in the example below.

@@ Line 19: / Line 19: @@
          &nbsp;
          c       UCVM Configuration Location
-                 CHARACTER(LEN=32) ucvmconf
+                 CHARACTER(LEN=80) ucvmconf
          c       Model Name
                  CHARACTER(LEN=4) model

@@ Line 75: / Line 75: @@
          &nbsp;
          c       Where is our configuration file?
-                 ucvmconf = "/home/scec-01/davidgil/ucvm.conf"
+                 ucvmconf = "/home/scec-01/davidgil/ucvm.conf" // CHAR(0)
          &nbsp;
          c       What model are we querying?

@@ Line 1: / Line 1: @@
-== Calling UCVM From Fortran==
+== Calling UCVM With Fortran==
-Calling
+Calling UCVM using Fortran is a relatively trivial procedure. In order to do so, you must install [http://scec.usc.edu/scecpedia/UCVM_User_Guide UCVM] as per the guide, including the etree and Proj.4 libraries. It is also necessary to download and install any community velocity models that you will be calling. UCVM requires those velocity models to be added to a configuration file so that they can be referenced within the application.
 == Example Code ==

← Older revision		Revision as of 22:39, 20 June 2013
Line 1:		Line 1:
−	== Calling UCVM ==	+	== Calling UCVM From Fortran==
		+
		+	Calling

	== Example Code ==		== Example Code ==