Changeset 856 for codes/icosagcm/devel/src/unstructured

Timestamp:

05/06/19 01:49:12 (5 years ago)

Author:

dubos

Message:

devel : towards Fortran driver for unstructured/LAM meshes

File:

• codes/icosagcm/devel/src/unstructured/init_unstructured.f90 (modified) (3 diffs)

codes/icosagcm/devel/src/unstructured/init_unstructured.f90

@@ -96 +96 @@
     USE ioipsl
     USE field_mod
+    USE geometry, ONLY : lon_i, lat_i, lon_e, lat_e, ep_e
     IMPLICIT NONE
 
     !!-------------Declare local variables-------------------
-    CHARACTER(LEN=*),PARAMETER :: f="mesh_information.nc"
-    INTEGER :: id_nc, ierr, status, descriptionLength
+    CHARACTER(LEN=*),PARAMETER :: f="input/mesh_information.nc"
+    INTEGER :: id_nc, ierr, status, descriptionLength, ij
     CHARACTER(LEN= 80) :: description
-
+    REAL(rstd), ALLOCATABLE :: angle_e(:)
+    REAL(rstd) :: clon, slon, clat, slat, & ! COS/SIN of lon/lat
+         elon(3), elat(3) ! lon/lat unit vectors
     print *,"------------------ READING FILE " , f, "----------------------- "
     !open and read the input file
@@ -167 +170 @@
     ALLOCATE(dual_vertex, source = Idata_read2)
 
+    CALL read_from_gridfile(id_nc, 'float', 'lon_i')
+    ALLOCATE(lon_i, source = Ddata_read1)
+    CALL read_from_gridfile(id_nc, 'float', 'lat_i')
+    ALLOCATE(lat_i, source = Ddata_read1)
+    CALL read_from_gridfile(id_nc, 'float', 'lon_e')
+    ALLOCATE(lon_e, source = Ddata_read1)
+    CALL read_from_gridfile(id_nc, 'float', 'lat_e')
+    ALLOCATE(lat_e, source = Ddata_read1)
+    CALL read_from_gridfile(id_nc, 'float', 'angle_e')
+    ALLOCATE(angle_e, source = Ddata_read1)
+
     CALL free_data_read ! free buffers after reading all data from grid file
 
@@ -176 +190 @@
     max_trisk_deg = SIZE(trisk,1)
 
+    ! now post-process some of the data we just read
+    ALLOCATE(ep_e(edge_num,3))
+    DO ij=1,edge_num
+       clon = COS(lon_e(ij))
+       slon = SIN(lon_e(ij))
+       clat = COS(lat_e(ij))
+       slat = SIN(lat_e(ij))
+       ! x,y,z = clat*clon, clat*slon, slat
+       elon = (/ -slon, clon, 0. /)
+       elat = (/ -slat*clon, -slat*slon, clat /)
+       ep_e(ij,:) = COS(angle_e(ij))*elon + SIN(angle_e(ij))*elat
+    END DO
+    
+    DEALLOCATE(angle_e)
   END SUBROUTINE read_dump_partition