Changeset 2475 in lmdz_wrf

Timestamp:

Apr 26, 2019, 8:53:56 PM (6 years ago)

Author:

lfita

Message:

Adding:

• `GDATE': Subroutine to compute the gregorian calendar date (year,month,day) given the julian date (JD)
• `JD': Fucntion to compute the julian date (JD) given a gregorian calendar

File:

• trunk/tools/module_generic.f90 (modified) (3 diffs)

trunk/tools/module_generic.f90

@@ -12 +12 @@
 !   by a list of grid-point coordinates
 ! GetInNamelist: Subroutine to get a paramter from a namelistfile
+! GDATE: Subroutine to compute the gregorian calendar date (year,month,day) given the julian date (JD)
 ! get_xyconlimits: Subroutine for getting the limits of contiguous values from a given point in a 2D matrix
 ! index_list_coordsI: Function to provide the index of a given coordinate within a list of integer coordinates
@@ -20 +21 @@
 ! Index2DArrayR: Function to provide the first index of a given value inside a 2D real array
 ! Index2DArrayR_K: Function to provide the first index of a given value inside a 2D real(r_k) array
+! JD: Fucntion to compute the julian date (JD) given a gregorian calendar
 ! Nvalues_2DArrayI: Number of different values of a 2D integer array
 ! mat2DPosition: Function to provide the i, j indices of a given value inside a 2D matrix
@@ -1143 +1145 @@
   END SUBROUTINE from_coordlist_2DRKmatrix
 
+  INTEGER FUNCTION JD (year,month,day)
+! Fucntion to compute the julian date (JD) given a gregorian calendar
+!   FROM: https://aa.usno.navy.mil/faq/docs/JD_Formula.php
+
+  IMPLICIT NONE
+
+  INTEGER, INTENT(in)                                    :: year,month,day
+
+! Local
+  INTEGER                                                :: I,J,K
+
+  I= year
+  J= month
+  K= day
+
+  PRINT *,1461*(I+4800+(J-14)/12)/4, 1461*(I+4800+(J-14)/12)/4.
+
+  JD= K-32075+1461*(I+4800+(J-14)/12)/4+367*(J-2-(J-14)/12*12)/12-3*((I+4900+(J-14)/12)/100)/4
+
+  END FUNCTION JD
+
+  SUBROUTINE GDATE(jd,year,month,day)
+  ! Subroutine to compute the gregorian calendar date (year,month,day) given the julian date (JD)
+
+    IMPLICIT NONE
+
+    INTEGER, INTENT(in)                                  :: jd
+    INTEGER, INTENT(out)                                 :: year,month,day
+! Local
+    INTEGER                                              ::  I,J,K,L,N
+
+    L= JD+68569
+    N= 4*L/146097
+    L= L-(146097*N+3)/4
+    I= 4000*(L+1)/1461001
+    L= L-1461*I/4+31
+    J= 80*L/2447
+    K= L-2447*J/80
+    L= J/11
+    J= J+2-12*L
+    I= 100*(N-49)+I+L
+
+    year= I
+    month= J
+    day= K
+
+    RETURN
+
+  END SUBROUTINE GDATE
+
 END MODULE module_generic