Changeset 8712

Timestamp:

2024-11-10T11:00:15+01:00 (7 weeks ago)

Author:

jan.polcher

Message:

Brought the defitions of the grid type in agreement with the latest versions of WRF. That essentially exchanges the code between LatLon? and Cassini. This change impacts the few users who have forcing data on a Cassini projection. They are already informed of the change.

File:

• branches/ORCHIDEE_2_2/ORCHIDEE/src_global/grid.f90 (modified) (3 diffs)

branches/ORCHIDEE_2_2/ORCHIDEE/src_global/grid.f90

@@ -735 +735 @@
     REAL(r_std)                :: user_stand_lon, user_truelat1, user_truelat2, user_dxkm, user_dykm
     REAL(r_std)                :: user_dlat, user_dlon, user_known_x, user_known_y, user_known_lat, user_known_lon
+    REAL(r_std)                :: user_lon0, user_lat0, proj_known_lon, proj_known_lat
     REAL(r_std), DIMENSION(16) :: corner_lons, corner_lats
     !
@@ -768 +769 @@
     iret = NF90_GET_ATT(fid, NF90_GLOBAL, "DX", user_dxkm)
     iret = NF90_GET_ATT(fid, NF90_GLOBAL, "DY", user_dykm)
-    user_dlat = undef
-    user_dlon = undef
+    !
+    user_dlat = user_dykm/(R_Earth * pi * 2.0 / 360.0)
+    user_dlon = user_dxkm/(R_Earth * pi * 2.0 / 360.0)
+    !
+    iret = NF90_GET_ATT(fid, NF90_GLOBAL, "POLE_LON", user_lon0)
+    iret = NF90_GET_ATT(fid, NF90_GLOBAL, "POLE_LAT", user_lat0)
+
+    
     !
     IF ( current_proj == PROJ_LATLON ) THEN
@@ -851 +858 @@
        
     else if (current_proj == PROJ_CASSINI) then
-       call ipslerr(3,"grid_initproj",'Should not have PROJ_CASSINI as projection for', &
-            'source data in push_source_projection()', " ")
-       
+       CALL rotate_coords(user_known_lat, user_known_lon, proj_known_lat, &
+            proj_known_lon, user_lat0, user_lon0, user_stand_lon, -1)
+       CALL map_set(current_proj, proj_stack(idom), &
+            loninc=user_dlon, &
+            latinc=user_dlat, &
+            lat1=user_known_lat, &
+            lon1=user_known_lon, &
+            lat0=user_lat0 ,&
+            lon0=user_lon0 ,&
+            knowni=user_known_x, &
+            knownj=user_known_y, &
+            stdlon=user_stand_lon, &
+            r_earth=R_Earth)
+
     else if (current_proj == PROJ_LC) then
        call map_set(current_proj, proj_stack(idom), &