Ignore:
Timestamp:
Jan 15, 2010, 11:35:59 AM (14 years ago)
Author:
Laurent Fairhead
Message:
  • sauvegarde de variables des thermiques dans le restart physique
  • modification de iostart pour pouvoir sauvegarder des champs en klevp1
  • nettoyage sur physiq.F pour gfortran
  • variables from the new version of the thermics are saved to the restart file
  • iostart is modified to allow the writing of fields dimensionned in klevp1
  • some clean-up work on physiq.F for the gfortran compilation
File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/phyetat0.F

    r1279 r1298  
    999999      PRINT*,'(ecart-type) wake_fip:', xmin, xmax
    10001000c
     1001c  thermiques
     1002c
     1003
     1004      CALL get_field("FM_THERM",fm_therm,found)
     1005      IF (.NOT. found) THEN
     1006         PRINT*, "phyetat0: Le champ <fm_therm> est absent"
     1007         PRINT*, "Depart legerement fausse. Mais je continue"
     1008         fm_therm=0.
     1009      ENDIF
     1010      xmin = 1.0E+20
     1011      xmax = -1.0E+20
     1012      xmin = MINval(fm_therm)
     1013      xmax = MAXval(fm_therm)
     1014      PRINT*,'(ecart-type) fm_therm:', xmin, xmax
     1015
     1016      CALL get_field("ENTR_THERM",entr_therm,found)
     1017      IF (.NOT. found) THEN
     1018         PRINT*, "phyetat0: Le champ <entr_therm> est absent"
     1019         PRINT*, "Depart legerement fausse. Mais je continue"
     1020         entr_therm=0.
     1021      ENDIF
     1022      xmin = 1.0E+20
     1023      xmax = -1.0E+20
     1024      xmin = MINval(entr_therm)
     1025      xmax = MAXval(entr_therm)
     1026      PRINT*,'(ecart-type) entr_therm:', xmin, xmax
     1027
     1028      CALL get_field("DETR_THERM",detr_therm,found)
     1029      IF (.NOT. found) THEN
     1030         PRINT*, "phyetat0: Le champ <detr_therm> est absent"
     1031         PRINT*, "Depart legerement fausse. Mais je continue"
     1032         detr_therm=0.
     1033      ENDIF
     1034      xmin = 1.0E+20
     1035      xmax = -1.0E+20
     1036      xmin = MINval(detr_therm)
     1037      xmax = MAXval(detr_therm)
     1038      PRINT*,'(ecart-type) detr_therm:', xmin, xmax
     1039
     1040
     1041
     1042c
    10011043c Read and send field trs to traclmdz
    10021044c
Note: See TracChangeset for help on using the changeset viewer.