Changeset 5095 for LMDZ6/branches/Amaury_dev/libf/phylmd/cosp2/optics_lib.F90

Timestamp:

Jul 22, 2024, 9:46:57 AM (4 months ago)

Author:

abarral

Message:

Revert cosp*/ from the trunk, as it's external code
Add missing bits from FCM2 source

File:

• LMDZ6/branches/Amaury_dev/libf/phylmd/cosp2/optics_lib.F90 (modified) (8 diffs)

LMDZ6/branches/Amaury_dev/libf/phylmd/cosp2/optics_lib.F90

@@ -539 +539 @@
     if (alam < cutice) then
        ! Region from 0.045 microns to 167.0 microns - no temperature depend
+       do i=2,nwl
+          if(alam < wl(i)) continue
+       enddo
        x1  = log(wl(i-1))
        x2  = log(wl(i))
@@ -555 +558 @@
        if(tk < temref(4)) tk=temref(4)
        do i=2,4
-          if(tk>=temref(i)) go to 12
+          if(tk.ge.temref(i)) go to 12
        enddo
 12     lt1 = i
        lt2 = i-1
        do i=2,nwlt
-          if(alam<=wlt(i)) go to 14
+          if(alam.le.wlt(i)) go to 14
        enddo
 14     x1  = log(wlt(i-1))
@@ -649 +652 @@
     Complex(wp) :: A1
     
-    If ((Dx>Imaxx) .Or. (InP>ImaxNP)) Then
+    If ((Dx.Gt.Imaxx) .Or. (InP.Gt.ImaxNP)) Then
        Error = 1
        Return
@@ -656 +659 @@
     Ir = 1 / Cm
     Y =  Dx * Cm
-    If (Dx<0.02) Then
+    If (Dx.Lt.0.02) Then
        NStop = 2
     Else
-       If (Dx<=8.0) Then
+       If (Dx.Le.8.0) Then
           NStop = Dx + 4.00*Dx**(1./3.) + 2.0
        Else
-          If (Dx< 4200.0) Then
+          If (Dx.Lt. 4200.0) Then
              NStop = Dx + 4.05*Dx**(1./3.) + 2.0
           Else
@@ -670 +673 @@
     End If
     NmX = Max(Real(NStop),Real(Abs(Y))) + 15.
-    If (Nmx > Itermax) then
+    If (Nmx .gt. Itermax) then
        Error = 1
        Return
@@ -723 +726 @@
 !ds       Dqxt = Tnp1 *      Dble(A + B)          + Dqxt
        Dqsc = Tnp1 * (A*Conjg(A) + B*Conjg(B)) + Dqsc
-       If (N>1) then
+       If (N.Gt.1) then
           Dg = Dg + (dN*dN - 1) * (ANM1*Conjg(A) + BNM1 * Conjg(B)) / dN + TNM1 *(ANM1*Conjg(BNM1)) / (dN*dN - dN)
 !ds          Dg = Dg + (dN*dN - 1) * Dble(ANM1*Conjg(A) + BNM1 * Conjg(B)) / dN + TNM1 * Dble(ANM1*Conjg(BNM1)) / (dN*dN - dN)
@@ -732 +735 @@
        AMB = A2 * (A - B)
        Do I = 1,Inp2
-          If (I>Inp) Then
+          If (I.GT.Inp) Then
              S(I) = -Pi1(I)
           Else
@@ -753 +756 @@
     End Do
 
-    If (Dg >0) Dg = 2 * Dg / Dqsc
+    If (Dg .GT.0) Dg = 2 * Dg / Dqsc
     Dqsc =  2 * Dqsc / Dx**2
     Dqxt =  2 * Dqxt / Dx**2