Changeset 4039 for trunk/LMDZ.MARS/libf/dynphy_lonlat

Timestamp:

Feb 4, 2026, 5:10:55 PM (12 days ago)

Author:

emillour

Message:

Mars PCM:
Minor updates/fixes on "xvik" to comply with picky gfortran compiler.
EM

File:

• trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/xvik.F (modified) (4 diffs)

trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/xvik.F

@@ -846 +846 @@
         if (iv==1) then
          open(ifile(iv), file = 'ps_VL1_year'//trim(chr2)//'_diurnal')
-         open(ifile(iv)+20, file = 'ps_VL1_year'//trim(chr2)//'_harmonics')
+         open(ifile(iv)+20, file = 
+     &                      'ps_VL1_year'//trim(chr2)//'_harmonics')
         
          IF(Time_unit == 1) THEN          
@@ -864 +865 @@
         elseif (iv==2) then
          open(ifile(iv), file = 'ps_VL2_year'//trim(chr2)//'_diurnal')
-         open(ifile(iv)+20, file = 'ps_VL2_year'//trim(chr2)//'_harmonics')
+         open(ifile(iv)+20, file = 
+     &                      'ps_VL2_year'//trim(chr2)//'_harmonics')
          
          IF(Time_unit == 1) THEN          
@@ -882 +884 @@
         else 
          open(ifile(iv), file = 'ps_INS_year'//trim(chr2)//'_diurnal')
-         open(ifile(iv)+20, file = 'ps_INS_year'//trim(chr2)//'_harmonics')
+         open(ifile(iv)+20, file = 
+     &                      'ps_INS_year'//trim(chr2)//'_harmonics')
          
          IF(Time_unit == 1) THEN          
@@ -944 +947 @@
         enddo 
          
-         write(ifile(iv)+20,'(a)') 'Fourrier coefficients ak/bk 
+        write(ifile(iv)+20,'(a)') 'Fourrier coefficients ak/bk 
      &(interpolated altitude)'
-         write(ifile(iv)+20,'(a,7f)') '#',(a_tab(k),k=1,n_harmo+1)
-         write(ifile(iv)+20,'(a,7f)') '#',(b_tab(k),k=1,n_harmo+1)
-         write(ifile(iv)+20,'(a)') 'Fourrier coefficients ak/bk 
+        write(ifile(iv)+20,'(a,7E12.4)') '#',(a_tab(k),k=1,n_harmo+1)
+        write(ifile(iv)+20,'(a,7E12.4)') '#',(b_tab(k),k=1,n_harmo+1)
+        write(ifile(iv)+20,'(a)') 'Fourrier coefficients ak/bk 
      &(GCM altitude)'
-         write(ifile(iv)+20,'(a,7f)') '#',(a_tab_gcm(k),k=1,n_harmo+1)
-         write(ifile(iv)+20,'(a,7f)') '#',(b_tab_gcm(k),k=1,n_harmo+1)   
+        write(ifile(iv)+20,'(a,7E12.4)')'#',(a_tab_gcm(k),k=1,n_harmo+1)
+        write(ifile(iv)+20,'(a,7E12.4)')'#',(b_tab_gcm(k),k=1,n_harmo+1)