Changeset 6508 for TOOLS/WATER_BUDGET/lmdz.py

Timestamp:

06/13/23 12:58:38 (13 months ago)

Author:

omamce

Message:

WaterUtils?.pyVersion qui marche pour :
Water Budget
O.M.

Version qui marche pour :

• Grille LMDZ et routage SECHIBA
• Grille ICO avec sorties natives et routage SIMPLE : routage pas très précis.

Ne marche pas pour :

• Grille LMDZ et routage SIMPLE : pb sur runoff
• Grille ICO avec sorties interpolées :

# Erreurs relatives

### VALID-CM622-LR.01 - LMDZ : OK ? Sauf LIC

• LIC : 1.e-2
• SRF : 4.e-6
• SRF/ATM : 2e-10
• RUNOFF : 1.5e-4

### VALID-CM622-SIMPLE-ROUTING - LMDZ : Erreur RUNOFF, LIC

• LIC : 1.6
• SRF : 1e-6
• SRF/ATM : 1.e-10
• RUNOFF : -7

### TEST-CM72-SIMPLE-ROUTING.13 - ICO : Erreur SRF, RUNOFF, LIC

• LIC : 150
• SRF : 0.5
• SRF/ATM : 4e-2
• RUNOFF : 700

### CM71v420-LR-pd-02.ini - ICO interpolé : Erreur SRF, RUNOFF, LIC

• LIC : 15
• SRF : 0.7
• SRF/ATM : -5.e-2
• ROUTING : 3e3

### CM71v420-LR-pd-02.ini - ICO natif : Erreurs faibles RUNOFF, LIC

• LIC : 0.3
• SRF : 7.e-6
• SRF/ATM : -6.e-9
• ROUTING : 5.e-2

File:

• TOOLS/WATER_BUDGET/lmdz.py (modified) (5 diffs)

TOOLS/WATER_BUDGET/lmdz.py

@@ -238 +238 @@
     if math == xr :
         p2D = xr.DataArray (p2D)
-        p2D.attrs = p1D.attrs
+        for attr in p1D.attrs : 
+            p2D.attrs[attr] = p1D.attrs[attr]
         p2D = p2D.rename ( {p2D.dims[0]:p1D.dims[0], p2D.dims[-1]:'x', p2D.dims[-2]:'y'} )
         
     return p2D
 
-def geo2point ( p2D, cumulPoles=False, dim1D='points_physiques' ) : 
+def geo2point ( p2D, cumulPoles=False, dim1D='cell', debug=False ) : 
     '''
     From 2D (lat, lon) to 1D (points_phyiques)
@@ -259 +260 @@
 
     jpn = jpi*(jpj-2) + 2
+
+    if debug :
+        print ( f'{jpi=} {jpj=} {jpn=} {jpt=}' )
 
     if jpt == -1 :
@@ -277 +281 @@
     if math == xr :
         p1D       = xr.DataArray (p1D)
-        p1D.attrs = p2D.attrs
+        for attr in p2D.attrs :
+            p1D.attrs[attr] = p2D.attrs[attr]
         p1D       = p1D.rename ( {p1D.dims[0]:p2D.dims[0], p1D.dims[-1]:dim1D} )
 
@@ -286 +291 @@
     return p1D
 
-def geo3point ( p3D, cumulPoles=False, dim1D='points_physiques' ) : 
+def geo3point ( p3D, cumulPoles=False, dim1D='cell' ) : 
     '''
     From 3D (lev, lat, lon) to 2D (lev, points_phyiques)
@@ -314 +319 @@
     if math == xr :
         p2D       = xr.DataArray (p2D)
-        p2D.attrs = p3D.attrs
+        for attr in p2D.attrs : 
+            p2D.attrs[attr] = p3D.attrs[attr]
         p2D       = p2D.rename ( {p2D.dims[-1]:dim1D, p2D.dims[-2]:p3D.dims[-3]} )