Changeset 5105 for LMDZ6/branches/Amaury_dev/libf/misc/juldate.f90

Timestamp:

Jul 23, 2024, 7:14:34 PM (8 weeks ago)

Author:

abarral

Message:

Replace 1DUTILS.h by module lmdz_1dutils.f90
Replace 1DConv.h by module lmdz_old_1dconv.f90 (it's only used by old_* files)
Convert *.F to *.f90
Fix gradsdef.h formatting
Remove unnecessary "RETURN" at the end of functions/subroutines

File:

• LMDZ6/branches/Amaury_dev/libf/misc/juldate.f90 (moved) (moved from LMDZ6/branches/Amaury_dev/libf/misc/juldate.F) (1 diff)

LMDZ6/branches/Amaury_dev/libf/misc/juldate.f90

@@ -2 +2 @@
 ! $Id$
 
-        SUBROUTINE juldate(ian,imoi,ijou,oh,om,os,tjd,tjdsec)
-c       Sous-routine de changement de date:
-c       gregorien>>>date julienne
-c       En entree:an,mois,jour,heure,min.,sec.
-c       En sortie:tjd
-        IMPLICIT NONE
-        INTEGER,INTENT(IN) :: ian,imoi,ijou,oh,om,os
-        REAL,INTENT(OUT) :: tjd,tjdsec
-        
-        REAL frac,year,rmon,cf,a,b
-        INTEGER ojou
-        
-        frac=((os/60.+om)/60.+oh)/24.
-        ojou=dble(ijou)+frac
-            year=dble(ian)
-            rmon=dble(imoi)
-        if (imoi <= 2) then
-            year=year-1.
-            rmon=rmon+12.
-        endif
-        cf=year+(rmon/100.)+(ojou/10000.)
-        if (cf >= 1582.1015) then
-            a=int(year/100)
-            b=2-a+int(a/4)
-        else
-            b=0
-        endif
-        tjd=int(365.25*year)+int(30.6001*(rmon+1))+int(ojou)
-     +   +1720994.5+b
-        tjdsec=(ojou-int(ojou))+(tjd-int(tjd))
-        tjd=int(tjd)+int(tjdsec)
-        tjdsec=tjdsec-int(tjdsec)
-        return
-        end
+  SUBROUTINE juldate(ian,imoi,ijou,oh,om,os,tjd,tjdsec)
+  ! Sous-routine de changement de date:
+  ! gregorien>>>date julienne
+  ! En entree:an,mois,jour,heure,min.,sec.
+  ! En sortie:tjd
+    IMPLICIT NONE
+    INTEGER,INTENT(IN) :: ian,imoi,ijou,oh,om,os
+    REAL,INTENT(OUT) :: tjd,tjdsec
+
+    REAL :: frac,year,rmon,cf,a,b
+    INTEGER :: ojou
+
+    frac=((os/60.+om)/60.+oh)/24.
+    ojou=dble(ijou)+frac
+        year=dble(ian)
+        rmon=dble(imoi)
+    if (imoi <= 2) then
+        year=year-1.
+        rmon=rmon+12.
+    endif
+    cf=year+(rmon/100.)+(ojou/10000.)
+    if (cf >= 1582.1015) then
+        a=int(year/100)
+        b=2-a+int(a/4)
+    else
+        b=0
+    endif
+    tjd=int(365.25*year)+int(30.6001*(rmon+1))+int(ojou) &
+          +1720994.5+b
+    tjdsec=(ojou-int(ojou))+(tjd-int(tjd))
+    tjd=int(tjd)+int(tjdsec)
+    tjdsec=tjdsec-int(tjdsec)
+
+end subroutine juldate