Changeset 3874 for trunk/LMDZ.COMMON/libf/dyn3d_common/disvert_noterre.F

Timestamp:

Aug 1, 2025, 4:03:00 PM (5 days ago)

Author:

afalco

Message:

check size of z2sig.def when reading it.
AF

File:

• trunk/LMDZ.COMMON/libf/dyn3d_common/disvert_noterre.F (modified) (5 diffs)

trunk/LMDZ.COMMON/libf/dyn3d_common/disvert_noterre.F

@@ -24 +24 @@
 c
 c=======================================================================
-c    Discretisation verticale en coordonnée hybride (ou sigma)
+c    Discretisation verticale en coordonn�e hybride (ou sigma)
 c
 c=======================================================================
@@ -138 +138 @@
          READ(99,*) scaleheight
          do l=1,llm
-            read(99,*) zsig(l)
+            read(99,*,iostat=ierr1) zsig(l)
+            if (ierr1 .ne. 0) then
+               write(lunout,*) 'ERROR with line ', l,' of z2sig.def'
+               call abort_gcm(modname, "z2sig.def too short?", 1)
+            endif
          end do
+         ! Check if there are more lines to read in the file
+         read(99,*,iostat=ierr1)
+         if (ierr1 .eq. 0) then
+            write(lunout,*) 'ERROR: Expected nb of levels:', llm
+            call abort_gcm(modname, 
+     &        'z2sig.def has more lines than expected', 1)
+         endif
          CLOSE(99)
 
@@ -208 +219 @@
 c     Calcul au milieu des couches :
 c     WARNING : le choix de placer le milieu des couches au niveau de
-c     pression intermédiaire est arbitraire et pourrait etre modifié.
+c     pression interm�diaire est arbitraire et pourrait etre modifi�.
 c     Le calcul du niveau pour la derniere couche 
 c     (on met la meme distance (en log pression)  entre P(llm)
 c     et P(llm -1) qu'entre P(llm-1) et P(llm-2) ) est
-c     Specifique.  Ce choix est spécifié ici ET dans exner_milieu.F
+c     Specifique.  Ce choix est sp�cifi� ici ET dans exner_milieu.F
 
       DO l = 1, llm-1
@@ -284 +295 @@
 c     L'objectif est de calculer newsig telle que
 c       (1 -pa/preff)*exp(1-1./newsig**2)+(pa/preff)*newsig = sig
-c     Cela ne se résoud pas analytiquement: 
-c     => on résoud par iterration bourrine 
+c     Cela ne se r�soud pas analytiquement: 
+c     => on r�soud par iterration bourrine 
 c     ----------------------------------------------
 c     Information  : where exp(1-1./x**2) become << x
@@ -320 +331 @@
               newsig=(X2+newsig)*0.5
           end if
-c         Test : on arete lorsque on approxime sig à moins de 0.01 m près 
+c         Test : on arete lorsque on approxime sig � moins de 0.01 m pr�s 
 c         (en pseudo altitude) :
           IF(abs(10.*log(F)).LT.1.E-5) goto 999