Changeset 1619 for LMDZ5/trunk/libf/phylmd/phyetat0.F

Timestamp:

Mar 27, 2012, 10:16:51 AM (13 years ago)

Author:

jghattas

Message:

• Changed variable ratqs, clwcon and rnebcon in restartphy file : before they were saved only at surface layer, now the whole variables are saved. Added specific option in iostart to be able to restart with

old restart files.

• Added variable uancien and vancien in restartphy.nc file. These variables are only needed for calculation of diagnostiques.

No change in result for standard physics (old or new).

File:

• LMDZ5/trunk/libf/phylmd/phyetat0.F (modified) (3 diffs)

LMDZ5/trunk/libf/phylmd/phyetat0.F

@@ -749 +749 @@
       ENDIF
 
-      u_ancien = 0.0   !AXC: We don't have u_ancien and v_ancien in the start
-      v_ancien = 0.0   !AXC: files, therefore they have to be initialized.
-c
+      CALL get_field("UANCIEN",u_ancien,found)
+      IF (.NOT. found) THEN
+         PRINT*, "phyetat0: Le champ <UANCIEN> est absent"
+         PRINT*, "Depart legerement fausse. Mais je continue"
+         ancien_ok = .FALSE.
+      ENDIF
+
+      CALL get_field("VANCIEN",v_ancien,found)
+      IF (.NOT. found) THEN
+         PRINT*, "phyetat0: Le champ <VANCIEN> est absent"
+         PRINT*, "Depart legerement fausse. Mais je continue"
+         ancien_ok = .FALSE.
+      ENDIF
 
       clwcon=0.
-      CALL get_field("CLWCON",clwcon(:,1),found)
+      CALL get_field("CLWCON",clwcon,found)
       IF (.NOT. found) THEN
          PRINT*, "phyetat0: Le champ CLWCON est absent"
@@ -766 +776 @@
 c
       rnebcon = 0.
-      CALL get_field("RNEBCON",rnebcon(:,1),found)
+      CALL get_field("RNEBCON",rnebcon,found)
       IF (.NOT. found) THEN
          PRINT*, "phyetat0: Le champ RNEBCON est absent"
@@ -781 +791 @@
 c
       ratqs=0.
-      CALL get_field("RATQS",ratqs(:,1),found)
+      CALL get_field("RATQS",ratqs,found)
       IF (.NOT. found) THEN
          PRINT*, "phyetat0: Le champ <RATQS> est absent"