Changeset 1048 for trunk/LMDZ.TITAN/libf/phytitan/interface_surf.F90

Timestamp:

Sep 23, 2013, 7:02:50 PM (11 years ago)

Author:

slebonnois

Message:

SL: update pour divers details titan + quelques modifs arch et makelmdz

File:

• trunk/LMDZ.TITAN/libf/phytitan/interface_surf.F90 (modified) (5 diffs)

trunk/LMDZ.TITAN/libf/phytitan/interface_surf.F90

@@ -45 +45 @@
       & tsol_rad, tsurf_new, alb_new)
 
+      use write_field_phy
+      use cpdet_mod, only: cpdet
+
+      IMPLICIT none
 
 ! Cette routine sert d'aiguillage entre l'atmosphere et la surface en general 
@@ -141 +145 @@
 
   if (check) write(*,*) 'Entree ', modname
-  
+
 ! Initialisations diverses
 !
@@ -149 +153 @@
 
 ! ADAPTATION GCM POUR CP(T)
-       zcp(1:klon) = cpdet(tsurf(1:klon))
+       do ii=1,klon
+         zcp(ii)=cpdet(tsurf(ii))
+       enddo
 
        IF (soil_model) THEN 
            CALL soil(dtime, knon, tsurf, tsoil,soilcap, soilflux)
            cal(1:knon) = zcp(1:knon) / soilcap(1:knon)
+! for tests:
+!  call writefield_phy('interfsurf_hq_zcp',zcp,1)
+!  call writefield_phy('interfsurf_hq_cal',cal,1)
+!  call writefield_phy('interfsurf_hq_soilcap',soilcap,1)
 !       print*,"DIAGNOSTIC SOIL"
 !       print*,"soilcap=",soilcap
@@ -186 +196 @@
      & tsurf_new, fluxsens, dflux_s)
 
+  use write_field_phy
+  use cpdet_mod, only: t2tpot, tpot2t
+
+  IMPLICIT none
+
 ! Cette routine calcule les fluxs en h a l'interface et eventuellement
 ! une temperature de surface (au cas ou ok_veget = false)
@@ -317 +332 @@
   ENDDO
 
+! for tests: write output fields...
+!  call writefield_phy('calcul_fluxs_d_ts',d_ts,1)
+!  call writefield_phy('calcul_fluxs_fluxsens',fluxsens,1)
+!  call writefield_phy('calcul_fluxs_dflux_s',dflux_s,1)
+
   END SUBROUTINE calcul_fluxs
 !